Nature, position, and frequency of mutations made in a single cycle of HIV-1 replication

HIV-1 Reverse Transcriptase Error Rates and Transcriptional Thresholds Based on Single-strand Consensus Sequencing of Target RNA Derived From In Vitro-transcription and HIV-infected Cells

Nucleotide incorporation and lacZ-based forward mutation assays have been widely used to determine the accuracy of reverse transcriptases (RTs) in RNA-dependent DNA polymerization reactions. However, they involve quite complex and laborious procedures, and cannot provide accurate error rates. Recently, NGS-based methods using barcodes opened the possibility of detecting all errors introduced by the RT, although their widespread use is limited by cost, due to the large size of libraries to be sequenced. In this study, we describe a novel and relatively simple NGS assay based on single-strand consensus sequencing that provides robust results with a relatively small number of raw sequences (around 60 Mb). The method has been validated by determining the error rate of HIV-1 (BH10 strain) RT using the HIV-1 protease-coding sequence as target. HIV-1 reverse transcription error rates in standard conditions (37 °C/3 mM Mg2+) using an in vitro-transcribed RNA were around 7.3 × 10-5. In agreement with previous reports, an 8-fold increase in RT's accuracy was observed after reducing Mg2+ concentration to 0.5 mM. The fidelity of HIV-1 RT was also higher at 50 °C than at 37 °C (error rate 1.5 × 10-5). Interestingly, error rates obtained with HIV-1 RNA from infected cells as template of the reverse transcription at 3 mM Mg2+ (7.4 × 10-5) were similar to those determined with the in vitro-transcribed RNA, and were reduced to 1.8 × 10-5 in the presence of 0.5 mM Mg2+. Values obtained at low magnesium concentrations were modestly higher than the transcription error rates calculated for human cells, thereby suggesting a realistic transcriptional threshold for our NGS-based error rate determinations.

Biological Barriers for Drug Delivery and Development of Innovative Therapeutic Approaches in HIV, Pancreatic Cancer, and Hemophilia A/B

Biological barriers remain a major obstacle for the development of innovative therapeutics. Depending on a disease's pathophysiology, the involved tissues, cell populations, and cellular components, drugs often have to overcome several biological barriers to reach their target cells and become effective in a specific cellular compartment. Human biological barriers are incredibly diverse and include multiple layers of protection and obstruction. Importantly, biological barriers are not only found at the organ/tissue level, but also include cellular structures such as the outer plasma membrane, the endolysosomal machinery, and the nuclear envelope. Nowadays, clinicians have access to a broad arsenal of therapeutics ranging from chemically synthesized small molecules, biologicals including recombinant proteins (such as monoclonal antibodies and hormones), nucleic-acid-based therapeutics, and antibody-drug conjugates (ADCs), to modern viral-vector-mediated gene therapy. In the past decade, the therapeutic landscape has been changing rapidly, giving rise to a multitude of innovative therapy approaches. In 2018, the FDA approval of patisiran paved the way for small interfering RNAs (siRNAs) to become a novel class of nucleic-acid-based therapeutics, which-upon effective drug delivery to their target cells-allow to elegantly regulate the post-transcriptional gene expression. The recent approvals of valoctocogene roxaparvovec and etranacogene dezaparvovec for the treatment of hemophilia A and B, respectively, mark the breakthrough of viral-vector-based gene therapy as a new tool to cure disease. A multitude of highly innovative medicines and drug delivery methods including mRNA-based cancer vaccines and exosome-targeted therapy is on the verge of entering the market and changing the treatment landscape for a broad range of conditions. In this review, we provide insights into three different disease entities, which are clinically, scientifically, and socioeconomically impactful and have given rise to many technological advancements: acquired immunodeficiency syndrome (AIDS) as a predominant infectious disease, pancreatic carcinoma as one of the most lethal solid cancers, and hemophilia A/B as a hereditary genetic disorder. Our primary objective is to highlight the overarching principles of biological barriers that can be identified across different disease areas. Our second goal is to showcase which therapeutic approaches designed to cross disease-specific biological barriers have been promising in effectively treating disease. In this context, we will exemplify how the right selection of the drug category and delivery vehicle, mode of administration, and therapeutic target(s) can help overcome various biological barriers to prevent, treat, and cure disease.

Combating antimicrobial resistance in malaria, HIV and tuberculosis

Antimicrobial resistance poses a significant threat to the sustainability of effective treatments against the three most prevalent infectious diseases: malaria, human immunodeficiency virus (HIV) infection and tuberculosis. Therefore, there is an urgent need to develop novel drugs and treatment protocols capable of reducing the emergence of resistance and combating it when it does occur. In this Review, we present an overview of the status and underlying molecular mechanisms of drug resistance in these three diseases. We also discuss current strategies to address resistance during the research and development of next-generation therapies. These strategies vary depending on the infectious agent and the array of resistance mechanisms involved. Furthermore, we explore the potential for cross-fertilization of knowledge and technology among these diseases to create innovative approaches for minimizing drug resistance and advancing the discovery and development of new anti-infective treatments. In conclusion, we advocate for the implementation of well-defined strategies to effectively mitigate and manage resistance in all interventions against infectious diseases.

Chronic HIV Transcription, Translation, and Persistent Inflammation

People with HIV exhibit persistent inflammation that correlates with HIV-associated comorbidities including accelerated aging, increased risk of cardiovascular disease, and neuroinflammation. Mechanisms that perpetuate chronic inflammation in people with HIV undergoing antiretroviral treatments are poorly understood. One hypothesis is that the persistent low-level expression of HIV proviruses, including RNAs generated from defective proviral genomes, drives the immune dysfunction that is responsible for chronic HIV pathogenesis. We explore factors during HIV infection that contribute to the generation of a pool of defective proviruses as well as how HIV-1 mRNA and proteins alter immune function in people living with HIV.

The Question of HIV Vaccine: Why Is a Solution Not Yet Available?

Ever since its discovery, human immunodeficiency virus type 1 (HIV-1) infection has remained a significant public health concern. The number of HIV-1 seropositive individuals currently stands at 40.1 million, yet definitive treatment for the virus is still unavailable on the market. Vaccination has proven to be a potent tool in combating infectious diseases, as evidenced by its success against other pathogens. However, despite ongoing efforts and research, the unique viral characteristics have prevented the development of an effective anti-HIV-1 vaccine. In this review, we aim to provide an historical overview of the various approaches attempted to create an effective anti-HIV-1 vaccine. Our objective is to explore the reasons why specific methods have failed to induce a protective immune response and to analyze the different modalities of immunogen presentation. This trial is registered with NCT05414786, NCT05471076, NCT04224701, and NCT01937455.

Modeling resistance to the broadly neutralizing antibody PGT121 in people living with HIV-1

PGT121 is a broadly neutralizing antibody in clinical development for the treatment and prevention of HIV-1 infection via passive administration. PGT121 targets the HIV-1 V3-glycan and demonstrated potent antiviral activity in a phase I clinical trial. Resistance to PGT121 monotherapy rapidly occurred in the majority of participants in this trial with the sampled rebound viruses being entirely resistant to PGT121 mediated neutralization. However, two individuals experienced long-term ART-free viral suppression following antibody infusion and retained sensitivity to PGT121 upon viral rebound. Here, we develop mathematical models of the HIV-1 dynamics during this phase I clinical trial. We utilize these models to understand the dynamics leading to PGT121 resistance and to identify the mechanisms driving the observed long-term viral control. Our modeling highlights the importance of the relative fitness difference between PGT121 sensitive and resistant subpopulations prior to treatment. Specifically, by fitting our models to data, we identify the treatment-induced competitive advantage of previously existing or newly generated resistant population as a primary driver of resistance. Finally, our modeling emphasizes the high neutralization ability of PGT121 in both participants who exhibited long-term viral control.

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.

From viral democratic genomes to viral wild bunch of quasispecies

The tremendous majority of RNA genomes from pathogenic viruses analyzed and deposited in databases are consensus or "democratic" genomes. They represent the genomes most frequently found in the clinical samples of patients but do not account for the huge genetic diversity of coexisting genomes, which is better described as quasispecies. A viral quasispecies is defined as the dynamic distribution of nonidentical but closely related mutants, variants, recombinant, or reassortant viral genomes. Viral quasispecies have collective behavior and dynamics and are the subject of internal interactions that comprise interference, complementation, or cooperation. In the setting of SARS-CoV-2 infection, intrahost SARS-CoV-2 genetic diversity was recently notably reported for immunocompromised, chronically infected patients, for patients treated with monoclonal antibodies targeting the viral spike protein, and for different body compartments of a single patient. A question that deserves attention is whether such diversity is generated postinfection from a clonal genome in response to selection pressure or is already present at the time of infection as a quasispecies. In the present review, we summarize the data supporting that hosts are infected by a "wild bunch" of viruses rather than by multiple virions sharing the same genome. Each virion in the "wild bunch" may have different virulence and tissue tropisms. As the number of viruses replicated during host infections is huge, a viral quasispecies at any time of infection is wide and is also influenced by host-specific selection pressure after infection, which accounts for the difficulty in deciphering and predicting the appearance of more fit variants and the evolution of epidemics of novel RNA viruses.

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.

Challenges and Progress in Designing Broad-Spectrum Vaccines Against Rapidly Mutating Viruses

Viruses evolve to evade prior immunity, causing significant disease burden. Vaccine effectiveness deteriorates as pathogens mutate, requiring redesign. This is a problem that has grown worse due to population increase, global travel, and farming practices. Thus, there is significant interest in developing broad-spectrum vaccines that mitigate disease severity and ideally inhibit disease transmission without requiring frequent updates. Even in cases where vaccines against rapidly mutating pathogens have been somewhat effective, such as seasonal influenza and SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2), designing vaccines that provide broad-spectrum immunity against routinely observed viral variation remains a desirable but not yet achieved goal. This review highlights the key theoretical advances in understanding the interplay between polymorphism and vaccine efficacy, challenges in designing broad-spectrum vaccines, and technology advances and possible avenues forward. We also discuss data-driven approaches for monitoring vaccine efficacy and predicting viral escape from vaccine-induced protection. In each case, we consider illustrative examples in vaccine development from influenza, SARS-CoV-2, and HIV (human immunodeficiency virus)-three examples of highly prevalent rapidly mutating viruses with distinct phylogenetics and unique histories of vaccine technology development.

Antigen-specific and cross-reactive T cells in protection and disease

Human T cells have a diverse T-cell receptor (TCR) repertoire that endows them with the ability to identify and defend against a broad spectrum of antigens. The universe of possible antigens that T cells may encounter, however, is even larger. To effectively surveil such a vast universe, the T-cell repertoire must adopt a high degree of cross-reactivity. Likewise, antigen-specific and cross-reactive T-cell responses play pivotal roles in both protective and pathological immune responses in numerous diseases. In this review, we explore the implications of these antigen-driven T-cell responses, with a particular focus on CD8+ T cells, using infection, neurodegeneration, and cancer as examples. We also summarize recent technological advances that facilitate high-throughput profiling of antigen-specific and cross-reactive T-cell responses experimentally, as well as computational biology approaches that predict these interactions.

HIV and Drug-Resistant Subtypes

Acquired Immunodeficiency Syndrome (AIDS) is a human viral infectious disease caused by the positive-sense single-stranded (ss) RNA Human Immunodeficiency Virus (HIV) (Retroviridae family, Ortervirales order). HIV-1 can be distinguished into various worldwide spread groups and subtypes. HIV-2 also causes human immunodeficiency, which develops slowly and tends to be less aggressive. HIV-2 only partially homologates to HIV-1 despite the similar derivation. Antiretroviral therapy (ART) is the treatment approved to control HIV infection, based on multiple antiretroviral drugs that belong to different classes: (i) NNRTIs, (ii) NRTIs, (iii) PIs, (iv) INSTIs, and (v) entry inhibitors. These drugs, acting on different stages of the HIV life cycle, decrease the patient's total burden of HIV, maintain the function of the immune system, and prevent opportunistic infections. The appearance of several strains resistant to these drugs, however, represents a problem today that needs to be addressed as best as we can. New outbreaks of strains show a widespread geographic distribution and a highly variable mortality rate, even affecting treated patients significantly. Therefore, novel treatment approaches should be explored. The present review discusses updated information on HIV-1- and HIV-2-resistant strains, including details on different mutations responsible for drug resistance.

Deep-sequencing of viral genomes from a large and diverse cohort of treatment-naive HIV-infected persons shows associations between intrahost genetic diversity and viral load

BACKGROUND

Infection with human immunodeficiency virus type 1 (HIV) typically results from transmission of a small and genetically uniform viral population. Following transmission, the virus population becomes more diverse because of recombination and acquired mutations through genetic drift and selection. Viral intrahost genetic diversity remains a major obstacle to the cure of HIV; however, the association between intrahost diversity and disease progression markers has not been investigated in large and diverse cohorts for which the majority of the genome has been deep-sequenced. Viral load (VL) is a key progression marker and understanding of its relationship to viral intrahost genetic diversity could help design future strategies for HIV monitoring and treatment.

METHODS

We analysed deep-sequenced viral genomes from 2,650 treatment-naive HIV-infected persons to measure the intrahost genetic diversity of 2,447 genomic codon positions as calculated by Shannon entropy. We tested for associations between VL and amino acid (AA) entropy accounting for sex, age, race, duration of infection, and HIV population structure.

RESULTS

We confirmed that the intrahost genetic diversity is highest in the env gene. Furthermore, we showed that mean Shannon entropy is significantly associated with VL, especially in infections of >24 months duration. We identified 16 significant associations between VL (p-value<2.0x10-5) and Shannon entropy at AA positions which in our association analysis explained 13% of the variance in VL. Finally, equivalent analysis based on variation in HIV consensus sequences explained only 2% of VL variance.

CONCLUSIONS

Our results elucidate that viral intrahost genetic diversity is associated with VL and could be used as a better disease progression marker than HIV consensus sequence variants, especially in infections of longer duration. We emphasize that viral intrahost diversity should be considered when studying viral genomes and infection outcomes.

TRIAL REGISTRATION

Samples included in this study were derived from participants who consented in the clinical trial, START (NCT00867048) (23), run by the International Network for Strategic Initiatives in Global HIV Trials (INSIGHT). All the participant sites are listed here: http://www.insight-trials.org/start/my_phpscript/participating.php?by=site.

Antiretroviral therapy resistance mutations among HIV infected people in Kazakhstan

In Kazakhstan, the number of people living with HIV (PLHIV) has increased steadily by 39% since 2010. Development of antiretroviral therapy (ART) resistance mutations (ARTRM) is a major hurdle in achieving effective treatment and prevention against HIV. Using HIV pol sequences from 602 PLHIV from Kazakhstan, we analyzed ARTRMs for their association with factors that may promote development of ARTRMs. 56% PLHIV were infected with HIV subtype A6 and 42% with CRF02_AG. The ARTRM Q174K was associated with increased viral load and decreased CD4+ cell count, while infection with CRF02_AG was associated with a lower likelihood of Q174K. Interestingly, CRF02_AG was positively associated with the ARTRM L10V that, in turn, was observed frequently with darunavir administration. Infection with CRF02_AG was positively associated with the ARTRM S162A that, in turn, was frequently observed with the administration of nevirapine, also associated with lower CD4 counts. Zidovudine or Nevirapine receipt was associated with the development of the ARTRM E138A, that, in turn, was associated with lower CD4 counts. Determination of a patient's HIV variant can help guide ART choice in Kazakhstan. For example, PLHIV infected with CRF02_AG will benefit less from darunavir and nevirapine, and emtricitabine should replace zidovudine.

[The prospects for the use of drugs based on the phenomenon of RNA interference against HIV infection]

The human immunodeficiency virus (HIV) is currently one of the most pressing global health problems. Since its discovery in 1978, HIV has claimed the lives of more than 35 million people, and the number of people infected today reaches 37 million. In the absence of highly active antiretroviral therapy (HAART), HIV infection is characterized by a steady decrease in the number of CD4+ T-lymphocytes, but its manifestations can affect the central nervous, cardiovascular, digestive, endocrine and genitourinary systems. At the same time, complications induced by representatives of pathogenic and opportunistic microflora, which can lead to the development of bacterial, fungal and viral concomitant infections, are of particular danger. It should be borne in mind that an important problem is the emergence of viruses resistant to standard therapy, as well as the toxicity of the drugs themselves for the body. In the context of this review, of particular interest is the assessment of the prospects for the creation and clinical use of drugs based on small interfering RNAs aimed at suppressing the reproduction of HIV, taking into account the experience of similar studies conducted earlier. RNA interference is a cascade of regulatory reactions in eukaryotic cells, which results in the degradation of foreign messenger RNA. The development of drugs based on the mechanism of RNA interference will overcome the problem of viral resistance. Along with this, this technology makes it possible to quickly respond to outbreaks of new viral diseases.

Mechanisms of Coronavirus Genome Stability As Potential Targets for Antiviral Drugs

The COVID-19 pandemic has made it necessary to create antivirals active against the SARS-CoV-2 coronavirus. One of the widely used strategies to fight off viral infections is the use of modified nucleoside analogues that inhibit viral replication by incorporating DNA or RNA into the growing chain, thus stopping its synthesis. The difficulty of using this method of treatment in the case of SARS-CoV-2 is that coronaviruses have an effective mechanism for maintaining genome stability. Its central element is the nsp14 protein, which is characterized by exonuclease activity, due to which incorrectly included and noncanonical nucleotides are removed from the 3' end of the growing RNA chain. Inhibitors of nsp14 exonuclease and nucleoside analogues resistant to its action are viewed as potential targets for anticoronavirus therapy.

Mechanistic Interplay between HIV-1 Reverse Transcriptase Enzyme Kinetics and Host SAMHD1 Protein: Viral Myeloid-Cell Tropism and Genomic Mutagenesis

Human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) has been the primary interest among studies on antiviral discovery, viral replication kinetics, drug resistance, and viral evolution. Following infection and entry into target cells, the HIV-1 core disassembles, and the viral RT concomitantly converts the viral RNA into double-stranded proviral DNA, which is integrated into the host genome. The successful completion of the viral life cycle highly depends on the enzymatic DNA polymerase activity of RT. Furthermore, HIV-1 RT has long been known as an error-prone DNA polymerase due to its lack of proofreading exonuclease properties. Indeed, the low fidelity of HIV-1 RT has been considered as one of the key factors in the uniquely high rate of mutagenesis of HIV-1, which leads to efficient viral escape from immune and therapeutic antiviral selective pressures. Interestingly, a series of studies on the replication kinetics of HIV-1 in non-dividing myeloid cells and myeloid specific host restriction factor, SAM domain, and HD domain-containing protein, SAMHD1, suggest that the myeloid cell tropism and high rate of mutagenesis of HIV-1 are mechanistically connected. Here, we review not only HIV-1 RT as a key antiviral target, but also potential evolutionary and mechanistic crosstalk among the unique enzymatic features of HIV-1 RT, the replication kinetics of HIV-1, cell tropism, viral genetic mutation, and host SAMHD1 protein.

Molecular Biology and Diversification of Human Retroviruses

Studies of retroviruses have led to many extraordinary discoveries that have advanced our understanding of not only human diseases, but also molecular biology as a whole. The most recognizable human retrovirus, human immunodeficiency virus type 1 (HIV-1), is the causative agent of the global AIDS epidemic and has been extensively studied. Other human retroviruses, such as human immunodeficiency virus type 2 (HIV-2) and human T-cell leukemia virus type 1 (HTLV-1), have received less attention, and many of the assumptions about the replication and biology of these viruses are based on knowledge of HIV-1. Existing comparative studies on human retroviruses, however, have revealed that key differences between these viruses exist that affect evolution, diversification, and potentially pathogenicity. In this review, we examine current insights on disparities in the replication of pathogenic human retroviruses, with a particular focus on the determinants of structural and genetic diversity amongst HIVs and HTLV.

Accelerated Evolution by Diversity-Generating Retroelements

Diversity-generating retroelements (DGRs) create vast amounts of targeted, functional diversity by facilitating the rapid evolution of ligand-binding protein domains. Thousands of DGRs have been identified in bacteria, archaea, and their respective viruses. They are broadly distributed throughout the microbial world, with enrichment observed in certain taxa and environments. The diversification machinery works through a novel mechanism termed mutagenic retrohoming, whereby nucleotide sequence information is copied from an invariant DNA template repeat (TR) into an RNA intermediate, selectively mutagenized at TR adenines during cDNA synthesis by a DGR-encoded reverse transcriptase, and transferred to a variable repeat (VR) region within a variable-protein gene (54). This unidirectional flow of information leaves TR-DNA sequences unmodified, allowing for repeated rounds of mutagenic retrohoming to optimize variable-protein function. DGR target genes are often modular and can encode one or more of a wide variety of discrete functional domains appended to a diversifiable ligand-binding motif. Bacterial variable proteins often localize to cell surfaces, although a subset appear to be cytoplasmic, while phage-encoded DGRs commonly diversify tail fiber-associated receptor-binding proteins. Here, we provide a comprehensive review of the mechanism and consequences of accelerated protein evolution by these unique and beneficial genetic elements. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Plasma-Derived HIV-1 Virions Contain Considerable Levels of Defective Genomes

Genetically-characterizing full-length HIV-1 RNA is critical for identifying genetically-intact genomes and for comparing these RNA genomes to proviral DNA. We have developed a method for sequencing plasma-derived RNA using long-range sequencing (PRLS assay; ∼8.3 kb from gag to the 3′ end or ∼5 kb from integrase to the 3′ end). We employed the gag-3′ PRLS assay to sequence HIV-1 RNA genomes from ART-naive participants during acute/early infection (n = 6) or chronic infection (n = 2). On average, only 65% of plasma-derived genomes were genetically-intact. Defects were found in all genomic regions but were concentrated in env and pol. We compared these genomes to near-full-length proviral sequences from paired peripheral blood mononuclear cell (PBMC) samples for the acute/early group and found that near-identical (>99.98% identical) sequences were identified only during acute infection. For three participants who initiated therapy during acute infection, we used the int-3′ PRLS assay to sequence plasma-derived genomes from an analytical treatment interruption and identified 100% identical genomes between pretherapy and rebound time points. The PRLS assay provides a new level of sensitivity for understanding the genetic composition of plasma-derived HIV-1 RNA from viremic individuals either pretherapy or after treatment interruption, which will be invaluable in assessing possible HIV-1 curative strategies.

IMPORTANCE We developed novel plasma-derived RNA using long-range sequencing assays (PRLS assay; 8.3 kb, gag-3′, and 5.0 kb, int-3′). Employing the gag-3′ PRLS assay, we found that 26% to 51% of plasma-derived genomes are genetically-defective, largely as a result of frameshift mutations and deletions. These genetic defects were concentrated in the env region compared to gag and pol, likely a reflection of viral immune escape in env during untreated HIV-1 infection. Employing the int-3′ PRLS assay, we found that analytical treatment interruption (ATI) plasma-derived sequences were identical and genetically-intact. Several sequences from the ATI plasma samples were identical to viral sequences from pretherapy plasma and PBMC samples, indicating that HIV-1 reservoirs established prior to therapy contribute to viral rebound during an ATI. Therefore, near-full-length sequencing of HIV-1 particles is required to gain an accurate picture of the genetic landscape of plasma HIV-1 virions in studies of HIV-1 replication and persistence.

Engineered Zinc Finger Protein Targeting 2LTR Inhibits HIV Integration in Hematopoietic Stem and Progenitor Cell-Derived Macrophages: In Vitro Study

Human hematopoietic stem/progenitor cell (HSPC)-based gene therapy is a promising direction for curing HIV-1-infected individuals. The zinc finger protein (2LTRZFP) designed to target the 2-LTR-circle junction of HIV-1 cDNA was previously reported as an intracellular antiviral molecular scaffold that prevents HIV integration. Here, we elucidate the efficacy and safety of using 2LTRZFP in human CD34⁺ HSPCs. We transduced 2LTRZFP which has the mCherry tag (2LTRZFPmCherry) into human CD34⁺ HSPCs using a lentiviral vector. The 2LTRZFPmCherry-transduced HSPCs were subsequently differentiated into macrophages. The expression levels of pro-apoptotic proteins of the 2LTRZFPmCherry-transduced HSPCs showed no significant difference from those of the non-transduced control. Furthermore, the 2LTRZFPmCherry-transduced HSPCs were successfully differentiated into mature macrophages, which had normal phagocytic function. The cytokine secretion assay demonstrated that 2LTRZFPmCherry-transduced CD34⁺ derived macrophages promoted the polarization towards classically activated (M1) subtypes. More importantly, the 2LTRZFPmCherry transduced cells significantly exhibited resistance to HIV-1 integration in vitro. Our findings demonstrate that the 2LTRZFPmCherry-transduced macrophages were found to be functionally and phenotypically normal, with no adverse effects of the anti-HIV-1 scaffold. Our data suggest that the anti-HIV-1 integrase scaffold is a promising antiviral molecule that could be applied to human CD34⁺ HSPC-based gene therapy for AIDS patients.

New Approaches to Multi-Parametric HIV-1 Genetics Using Multiple Displacement Amplification: Determining the What, How, and Where of the HIV-1 Reservoir

Development of potential HIV-1 curative interventions requires accurate characterization of the proviral reservoir, defined as host-integrated viral DNA genomes that drive rebound of viremia upon halting ART (antiretroviral therapy). Evaluation of such interventions necessitates methods capable of pinpointing the rare, genetically intact, replication-competent proviruses within a background of defective proviruses. This evaluation can be achieved by identifying the distinct integration sites of intact proviruses within host genomes and monitoring the dynamics of these proviruses and host cell lineages over longitudinal sampling. Until recently, molecular genetic approaches at the single proviral level have been generally limited to one of a few metrics, such as proviral genome sequence/intactness, host-proviral integration site, or replication competency. New approaches, taking advantage of MDA (multiple displacement amplification) for WGA (whole genome amplification), have enabled multiparametric proviral characterization at the single-genome level, including proviral genome sequence, host-proviral integration site, and phenotypic characterization of the host cell lineage, such as CD4 memory subset and antigen specificity. In this review, we will examine the workflow of MDA-augmented molecular genetic approaches to study the HIV-1 reservoir, highlighting technical advantages and flexibility. We focus on a collection of recent studies in which investigators have used these approaches to comprehensively characterize intact and defective proviruses from donors on ART, investigate mechanisms of elite control, and define cell lineage identity and antigen specificity of infected CD4+ T cell clones. The highlighted studies exemplify how these approaches and their future iterations will be key in defining the targets and evaluating the impacts of HIV curative interventions.

Physiological magnesium concentrations increase fidelity of diverse reverse transcriptases from HIV-1, HIV-2, and foamy virus, but not MuLV or AMV

Reverse transcriptases (RTs) are typically assayed using optimized Mg2+ concentrations (~5-10 mM) several-fold higher than physiological cellular free Mg2+ (~0.5 mM). Recent analyses demonstrated that HIV-1, but not Moloney murine leukaemia (MuLV) or avain myeloblastosis (AMV) virus RTs has higher fidelity in low Mg2+. In the current report, lacZα-based α-complementation assays were used to measure the fidelity of several RTs including HIV-1 (subtype B and A/E), several drug-resistant HIV-1 derivatives, HIV-2, and prototype foamy virus (PFV), all which showed higher fidelity using physiological Mg2+, while MuLV and AMV RTs demonstrated equivalent fidelity in low and high Mg2+. In 0.5 mM Mg2+, all RTs demonstrated approximately equal fidelity, except for PFV which showed higher fidelity. A Next Generation Sequencing (NGS) approach that used barcoding to determine mutation profiles was used to examine the types of mutations made by HIV-1 RT (type B) in low (0.5 mM) and high (6 mM) Mg2+ on a lacZα template. Unlike α-complementation assays which are dependent on LacZα activity, the NGS assay scores mutations at all positions and of every type. Consistent with α-complementation assays, a ~four-fold increase in mutations was observed in high Mg2+. These findings help explain why HIV-1 RT displays lower fidelity in vitro (with high Mg2+ concentrations) than other RTs (e.g. MuLV and AMV), yet cellular fidelity for these viruses is comparable. Establishing in vitro conditions that accurately represent RT's activity in cells is pivotal to determining the contribution of RT and other factors to the mutation profile observed with HIV-1.

virusMED: an atlas of hotspots of viral proteins

Metal binding sites, antigen epitopes and drug binding sites are the hotspots in viral proteins that control how viruses interact with their hosts. virusMED (virus Metal binding sites, Epitopes and Drug binding sites) is a rich internet application based on a database of atomic interactions around hotspots in 7041 experimentally determined viral protein structures. 25306 hotspots from 805 virus strains from 75 virus families were characterized, including influenza, HIV-1 and SARS-CoV-2 viruses. Just as Google Maps organizes and annotates points of interest, virusMED presents the positions of individual hotspots on each viral protein and creates an atlas upon which newly characterized functional sites can be placed as they are being discovered. virusMED contains an extensive set of annotation tags about the virus species and strains, viral hosts, viral proteins, metal ions, specific antibodies and FDA-approved drugs, which permits rapid screening of hotspots on viral proteins tailored to a particular research problem. The virusMED portal (https://virusmed.biocloud.top) can serve as a window to a valuable resource for many areas of virus research and play a critical role in the rational design of new preventative and therapeutic agents targeting viral infections.

Understanding patterns of HIV multi-drug resistance through models of temporal and spatial drug heterogeneity

Triple-drug therapies have transformed HIV from a fatal condition to a chronic one. These therapies should prevent HIV drug resistance evolution, because one or more drugs suppress any partially resistant viruses. In practice, such therapies drastically reduced, but did not eliminate, resistance evolution. In this article, we reanalyze published data from an evolutionary perspective and demonstrate several intriguing patterns about HIV resistance evolution - resistance evolves (1) even after years on successful therapy, (2) sequentially, often via one mutation at a time and (3) in a partially predictable order. We describe how these observations might emerge under two models of HIV drugs varying in space or time. Despite decades of work in this area, much opportunity remains to create models with realistic parameters for three drugs, and to match model outcomes to resistance rates and genetic patterns from individuals on triple-drug therapy. Further, lessons from HIV may inform other systems.

Development of a User-Friendly Pipeline for Mutational Analyses of HIV Using Ultra-Accurate Maximum-Depth Sequencing

Human immunodeficiency virus type 2 (HIV-2) accumulates fewer mutations during replication than HIV type 1 (HIV-1). Advanced studies of HIV-2 mutagenesis, however, have historically been confounded by high background error rates in traditional next-generation sequencing techniques. In this study, we describe the adaptation of the previously described maximum-depth sequencing (MDS) technique to studies of both HIV-1 and HIV-2 for the ultra-accurate characterization of viral mutagenesis. We also present the development of a user-friendly Galaxy workflow for the bioinformatic analyses of sequencing data generated using the MDS technique, designed to improve replicability and accessibility to molecular virologists. This adapted MDS technique and analysis pipeline were validated by comparisons with previously published analyses of the frequency and spectra of mutations in HIV-1 and HIV-2 and is readily expandable to studies of viral mutation across the genomes of both viruses. Using this novel sequencing pipeline, we observed that the background error rate was reduced 100-fold over standard Illumina error rates, and 10-fold over traditional unique molecular identifier (UMI)-based sequencing. This technical advancement will allow for the exploration of novel and previously unrecognized sources of viral mutagenesis in both HIV-1 and HIV-2, which will expand our understanding of retroviral diversity and evolution.

HIV-1 proviral landscape characterization varies by pipeline analysis

INTRODUCTION

HIV rebounds after cessation of antiretroviral therapy, representing a barrier to cure. To better understand the virus reservoir, analysis pipelines have been developed that categorize proviral sequences as intact or defective, and further determine the precise nature of the sequence defects that may be present. We investigated the effects that different analysis pipelines had on the characterization of HIV-1 proviral sequences.

METHODS

We used single genome amplification to generate near full-length (NFL) HIV-1 proviral DNA sequences, defined as amplicons greater than 8000 base pairs in length, isolated from peripheral blood mononuclear cells (PBMC) of treated suppressed participants with HIV-1. Amplicons underwent direct next-generation single genome sequencing and were analysed using four HIV-1 proviral characterization pipelines. Sequences were characterized as intact or defective; defective sequences were assessed for the number and types of defects present. To confirm and extend our findings, 691 proviruses from the Proviral Sequence Database (PSD) were analysed and the ProSeq-IT tool of the PSD was used to characterize both the participant and PSD proviruses.

RESULTS AND DISCUSSION

Virus sequences derived from thirteen ART-treated virologically suppressed participants with HIV were studied. A total of 693 HIV-1 proviral sequences were generated, 282 of which were NFL. An average of 53 sequences per participant was analysed. We found that proviruses often harbour multiple sequence defect types (mean 2.7, 95% confidence interval [CI] 2.5, 3.0); the elimination order used by each pipeline affected the percentage of proviruses allotted into each defect category. These differences varied between participants, depending on the number of defect categories present in a given provirus sequence. Pipeline-specific differences in characterizing the HIV-1 5' untranslated region (5' UTR) led to an overestimation of the number of intact NFL proviral sequences, a finding corroborated in the independent PSD analysis. A comparison of the four published pipelines to ProSeq-IT found that ProSeq IT was more likely to characterize proviruses as intact.

CONCLUSIONS

The choice of pipeline used for HIV-1 provirus landscape analysis may bias the classification of defective sequences. To improve the comparison of provirus characterizations across research groups, the development of a consensus elimination pipeline should be prioritized.

TRAL 2.0: Tandem Repeat Detection With Circular Profile Hidden Markov Models and Evolutionary Aligner

The Tandem Repeat Annotation Library (TRAL) focuses on analyzing tandem repeat units in genomic sequences. TRAL can integrate and harmonize tandem repeat annotations from a large number of external tools, and provides a statistical model for evaluating and filtering the detected repeats. TRAL version 2.0 includes new features such as a module for identifying repeats from circular profile hidden Markov models, a new repeat alignment method based on the progressive Poisson Indel Process, an improved installation procedure and a docker container. TRAL is an open-source Python 3 library and is available, together with documentation and tutorials viavital-it.ch/software/tral.

Identification of a Novel Cis-Acting Regulator of HIV-1 Genome Packaging

Human immunodeficiency virus type 1 (HIV-1) uptakes homo-dimerized viral RNA genome into its own particle. A cis-acting viral RNA segment responsible for this event, termed packaging signal (psi), is located at the 5′-end of the viral genome. Although the psi segment exhibits nucleotide variation in nature, its effects on the psi function largely remain unknown. Here we show that a psi sequence from an HIV-1 regional variant, subtype D, has a lower packaging ability compared with that from another regional variant, HIV-1 subtype B, despite maintaining similar genome dimerization activities. A series of molecular genetic investigations narrowed down the responsible element of the selective attenuation to the two sequential nucleotides at positions 226 and 227 in the psi segment. Molecular dynamics simulations predicted that the dinucleotide substitution alters structural dynamics, fold, and hydrogen-bond networks primarily of the psi-SL2 element that contains the binding interface of viral nucleocapsid protein for the genome packaging. In contrast, such structural changes were minimal within the SL1 element involved in genome dimerization. These results suggest that the psi 226/227 dinucleotide pair functions as a cis-acting regulator to control the psi structure to selectively tune the efficiency of packaging, but not dimerization of highly variable HIV-1 genomes.

The Role of Tissue Resident Memory CD4 T Cells in Herpes Simplex Viral and HIV Infection

Tissue-resident memory T cells (TRM) were first described in 2009. While initially the major focus was on CD8+ TRM, there has recently been increased interest in defining the phenotype and the role of CD4+ TRM in diseases. Circulating CD4+ T cells seed CD4+ TRM, but there also appears to be an equilibrium between CD4+ TRM and blood CD4+ T cells. CD4+ TRM are more mobile than CD8+ TRM, usually localized deeper within the dermis/lamina propria and yet may exhibit synergy with CD8+ TRM in disease control. This has been demonstrated in herpes simplex infections in mice. In human recurrent herpes infections, both CD4+ and CD8+ TRM persisting between lesions may control asymptomatic shedding through interferon-gamma secretion, although this has been more clearly shown for CD8+ T cells. The exact role of the CD4+/CD8+ TRM axis in the trigeminal ganglia and/or cornea in controlling recurrent herpetic keratitis is unknown. In HIV, CD4+ TRM have now been shown to be a major target for productive and latent infection in the cervix. In HSV and HIV co-infections, CD4+ TRM persisting in the dermis support HIV replication. Further understanding of the role of CD4+ TRM and their induction by vaccines may help control sexual transmission by both viruses.

Accelerating phylogeny-aware alignment with indel evolution using short time Fourier transform

Recently we presented a frequentist dynamic programming (DP) approach for multiple sequence alignment based on the explicit model of indel evolution Poisson Indel Process (PIP). This phylogeny-aware approach produces evolutionary meaningful gap patterns and is robust to the ‘over-alignment’ bias. Despite linear time complexity for the computation of marginal likelihoods, the overall method’s complexity is cubic in sequence length. Inspired by the popular aligner MAFFT, we propose a new technique to accelerate the evolutionary indel based alignment. Amino acid sequences are converted to sequences representing their physicochemical properties, and homologous blocks are identified by multi-scale short-time Fourier transform. Three three-dimensional DP matrices are then created under PIP, with homologous blocks defining sparse structures where most cells are excluded from the calculations. The homologous blocks are connected through intermediate ‘linking blocks’. The homologous and linking blocks are aligned under PIP as independent DP sub-matrices and their tracebacks merged to yield the final alignment. The new algorithm can largely profit from parallel computing, yielding a theoretical speed-up estimated to be proportional to the cubic power of the number of sub-blocks in the DP matrices. We compare the new method to the original PIP approach and demonstrate it on real data.

Avoiding Drug Resistance in HIV Reverse Transcriptase

HIV reverse transcriptase (RT) is an enzyme that plays a major role in the replication cycle of HIV and has been a key target of anti-HIV drug development efforts. Because of the high genetic diversity of the virus, mutations in RT can impart resistance to various RT inhibitors. As the prevalence of drug resistance mutations is on the rise, it is necessary to design strategies that will lead to drugs less susceptible to resistance. Here we provide an in-depth review of HIV reverse transcriptase, current RT inhibitors, novel RT inhibitors, and mechanisms of drug resistance. We also present novel strategies that can be useful to overcome RT's ability to escape therapies through drug resistance. While resistance may not be completely avoidable, designing drugs based on the strategies and principles discussed in this review could decrease the prevalence of drug resistance.

Survey of Pretreatment HIV Drug Resistance and Genetic Transmission Network Analysis Among HIV Patients in a High Drug-Use Area of Southwest China

Background: Pretreatment drug resistance (PDR) poses an increasing threat to the success of antiretroviral treatment (ART) programs in China. We aimed to conduct a survey of PDR among HIV patients in an area in Southwest China with extensive drug trafficking.

Methods: Consecutive cross-sectional surveys were conducted in Liangshan Prefecture of Sichuan Province from 2009 to 2018 based on the WHO-recommended method. PDR was identified by testing pol region sequences with the Stanford HIVdb algorithm (version 7.0). PDR prevalence and related factors were assessed by multivariable logistic regression. The transmission of HIV drug resistance was analyzed using a genetic transmission network.

Results: HIV-1 pol genes from 1889 patients were successfully amplified. The distribution of HIV-1 genotypes was as follows: CRF07_BC (94.0%), CRF08_BC (2.3%), CRF01_AE (2.0%) and others (1.4%). Of the participants, 6.9% (95% CI: 4.1-8.1%) had pretreatment resistance to 12 antiretroviral drugs recommended by the WHO, and nucleoside reverse transcriptase inhibitor (NRTI), non-nucleoside reverse transcriptase inhibitor (NNRTI) and protease inhibitors (PI) resistance were identified among 1.4% (95% CI: 0.7-3.4%), 5.8% (95% CI: 1.2-8.7%) and 0.4% (95% CI: 0.1-3.0%) of the patients, respectively. In the multivariate logistic model, the prevalence of PDR was 1.52-fold higher among intravenous drug users (IDUs) than among patients infected by heterosexual transmission (95% CI: 1.07-2.38; P=0.049), and the prevalence of PDR among patients diagnosed from 2017-2018 was 2.03-fold higher than that among patients diagnosed from 2009-2016 (95% CI: 1.18-5.76; P=0.018). A total of 26 clusters containing PDR and a rapidly growing drug resistance-related cluster containing the E138Q and V179D mutations were identified by genetic transmission network analysis.

Conclusion: The results show a moderate overall level of PDR prevalence and rapidly growing drug resistance over time. Preventive intervention should be focused on controlling the HIV epidemic among drug users, and surveillance is urgently needed to monitor the trend of PDR.

The clarifying role of time series data in the population genetics of HIV

HIV can evolve remarkably quickly in response to antiretroviral therapies and the immune system. This evolution stymies treatment effectiveness and prevents the development of an HIV vaccine. Consequently, there has been a great interest in using population genetics to disentangle the forces that govern the HIV adaptive landscape (selection, drift, mutation, and recombination). Traditional population genetics approaches look at the current state of genetic variation and infer the processes that can generate it. However, because HIV evolves rapidly, we can also sample populations repeatedly over time and watch evolution in action. In this paper, we demonstrate how time series data can bound evolutionary parameters in a way that complements and informs traditional population genetic approaches. Specifically, we focus on our recent paper (Feder et al., 2016, eLife), in which we show that, as improved HIV drugs have led to fewer patients failing therapy due to resistance evolution, less genetic diversity has been maintained following the fixation of drug resistance mutations. Because soft sweeps of multiple drug resistance mutations spreading simultaneously have been previously documented in response to the less effective HIV therapies used early in the epidemic, we interpret the maintenance of post-sweep diversity in response to poor therapies as further evidence of soft sweeps and therefore a high population mutation rate (θ) in these intra-patient HIV populations. Because improved drugs resulted in rarer resistance evolution accompanied by lower post-sweep diversity, we suggest that both observations can be explained by decreased population mutation rates and a resultant transition to hard selective sweeps. A recent paper (Harris et al., 2018, PLOS Genetics) proposed an alternative interpretation: Diversity maintenance following drug resistance evolution in response to poor therapies may have been driven by recombination during slow, hard selective sweeps of single mutations. Then, if better drugs have led to faster hard selective sweeps of resistance, recombination will have less time to rescue diversity during the sweep, recapitulating the decrease in post-sweep diversity as drugs have improved. In this paper, we use time series data to show that drug resistance evolution during ineffective treatment is very fast, providing new evidence that soft sweeps drove early HIV treatment failure.

Spontaneous Mutations in HIV-1 Gag, Protease, RT p66 in the First Replication Cycle and How They Appear: Insights from an In Vitro Assay on Mutation Rates and Types

While drug resistant mutations in HIV-1 are largely credited to its error prone HIV-1 RT, the time point in the infection cycle that these mutations can arise and if they appear spontaneously without selection pressures both remained enigmatic. Many HIV-1 RT mutational in vitro studies utilized reporter genes (LacZ) as a template to investigate these questions, thereby not accounting for the possible contribution of viral codon usage. To address this gap, we investigated HIV-1 RT mutation rates and biases on its own Gag, protease, and RT p66 genes in an in vitro selection pressure free system. We found rare clinical mutations with a general avoidance of crucial functional sites in the background mutations rates for Gag, protease, and RT p66 at 4.71 × 10-5, 6.03 × 10-5, and 7.09 × 10-5 mutations/bp, respectively. Gag and p66 genes showed a large number of 'A to G' mutations. Comparisons with silently mutated p66 sequences showed an increase in mutation rates (1.88 × 10-4 mutations/bp) and that 'A to G' mutations occurred in regions reminiscent of ADAR neighbor sequence preferences. Mutational free energies of the 'A to G' mutations revealed an avoidance of destabilizing effects, with the natural p66 gene codon usage providing barriers to disruptive amino acid changes. Our study demonstrates the importance of studying mutation emergence in HIV genes in a RT-PCR in vitro selection pressure free system to understand how fast drug resistance can emerge, providing transferable applications to how new viral diseases and drug resistances can emerge.

Determining the Origins of Human Immunodeficiency Virus Type 1 Drug-resistant Minority Variants in People Who Are Recently Infected Using Phylogenetic Reconstruction

BACKGROUND

Drug-resistant minority variants (DRMinVs) detected in patients who recently acquired human immunodeficiency virus type 1 (HIV-1) can be transmitted, generated de novo through virus replication, or technical errors. The first form is likely to persist and result in treatment failure, while the latter two could be stochastic and transient.

METHODS

Ultradeep sequencing of plasma samples from 835 individuals with recent HIV-1 infection in the United Kingdom was performed to detect DRMinVs at a mutation frequency between 2% and 20%. Sequence alignments including >110 000 HIV-1 partial pol consensus sequences from the UK HIV Drug Resistance Database (UK-HDRD), linked to epidemiological and clinical data from the HIV and AIDS Reporting System, were used for transmission cluster analysis. Transmission clusters were identified using Cluster Picker with a clade support of >90% and maximum genetic distances of 4.5% or 1.5%, the latter to limit detection to likely direct transmission events.

RESULTS

Drug-resistant majority variants (DRMajVs) were detected in 66 (7.9%) and DRMinVs in 84 (10.1%) of the recently infected individuals. High levels of clustering to sequences in UK-HDRD were observed for both DRMajV (n = 48; 72.7%) and DRMinV (n = 63; 75.0%) sequences. Of these, 43 (65.2%) with DRMajVs were in a transmission cluster with sequences that harbored the same DR mutation compared to only 3 (3.6%) sequences with DRMinVs (P < .00001, Fisher exact test). Evidence of likely direct transmission of DRMajVs was observed for 25/66 (37.9%), whereas none were observed for the DRMinVs (P < .00001).

CONCLUSIONS

Using a densely sampled HIV-infected population, we show no evidence of DRMinV transmission among recently infected individuals.

Structural investigation of 2-naphthyl phenyl ether inhibitors bound to WT and Y181C reverse transcriptase highlights key features of the NNRTI binding site

Human immunodeficiency virus (HIV)-1 remains as a global health issue that is primarily treated with highly active antiretroviral therapy, a combination of drugs that target the viral life cycle. One class of these drugs are non-nucleoside reverse transcriptase inhibitors (NNRTIs) that target the viral reverse transcriptase (RT). First generation NNRTIs were troubled with poor pharmacological properties and drug resistance, incentivizing the development of improved compounds. One class of developed compounds are the 2-naphthyl phenyl ethers, showing promising efficacy against the Y181C RT mutation. Further biochemical and structural work demonstrated differences in potency against the Y181C mutation and binding mode of the compounds. This work aims to understand the relationship between the binding mode and ability to overcome drug resistance using macromolecular x-ray crystallography. Comparison of 2-naphthyl phenyl ethers bound to Y181C RT reveal that compounds that interact with the invariant W229 are more capable of retaining efficacy against the resistance mutation. Additional modifications to these compounds at the 4-position, computationally designed to compensate for the Y181C mutation, do not demonstrate improved potency. Ultimately, we highlight important considerations for the development of future HIV-1 drugs that are able to combat drug resistance.

Deficient uracil base excision repair leads to persistent dUMP in HIV proviruses during infection of monocytes and macrophages

Non-dividing cells of the myeloid lineage such as monocytes and macrophages are target cells of HIV that have low dNTP pool concentrations and elevated levels of dUTP, which leads to frequent incorporation of dUMP opposite to A during reverse transcription ("uracilation"). One factor determining the fate of dUMP in proviral DNA is the host cell uracil base excision repair (UBER) system. Here we explore the relative UBER capacity of monocytes (MC) and monocyte-derived macrophages (MDM) and the fate of integrated uracilated viruses in both cell types to understand the implications of viral dUMP on HIV diversification and infectivity. We find that the kinetics for MC infection is compatible with their lifetime in vivo and their near absence of hUNG2 activity is consistent with the retention of viral dUMP at high levels at least until differentiation into macrophages, where UBER becomes possible. Overexpression of human uracil DNA glycosylase in MDM prior to infection resulted in rapid removal of dUMP from HIV cDNA and near complete depletion of dUMP-containing viral copies. This finding establishes that the low hUNG2 expression level in these cells limits UBER but that hUNG2 is restrictive against uracilated viruses. In contrast, overexpression of hUNG2 after viral integration did not accelerate the excision of uracils, suggesting that they may poorly accessible in the context of chromatin. We found that viral DNA molecules with incorporated dUMP contained unique (+) strand transversion mutations that were not observed when dUMP was absent (G→T, T→A, T→G, A→C). These observations and other considerations suggest that dUMP introduces errors predominantly during (-) strand synthesis when the template is RNA. Overall, the likelihood of producing a functional virus from in vitro infection of MC is about 50-fold and 300-fold reduced as compared to MDM and activated T cells. The results implicate viral dUMP incorporation in MC and MDM as a potential viral diversification and restriction pathway during human HIV infection.

Drug Resistance Evolution in HIV in the Late 1990s: Hard Sweeps, Soft Sweeps, Clonal Interference and the Accumulation of Drug Resistance Mutations

The evolution of drug resistance in pathogens such as HIV is an important and widely known example in the field of evolutionary medicine. Here, we focus on a unique data set from the late 1990s with multiple viral sequences from multiple time points in 118 patients. We study patterns of evolutionary dynamics in the viral populations in these patients who were treated with Reverse Transcriptase Inhibitors and Protease Inhibitors in the late 1990s. Specifically, we aim to visualize and analyze examples of population genetic processes such as selective sweeps and clonal interference. The figures and descriptions in this paper can be used in evolution and population genetics classes. We show and analyze a wide variety of patterns, specifically: soft sweeps, hard sweeps, softening sweeps and hardening sweeps, simultaneous sweeps, accumulation of mutations and clonal interference.

The Determination of HIV-1 RT Mutation Rate, Its Possible Allosteric Effects, and Its Implications on Drug Resistance

The high mutation rate of the human immunodeficiency virus type 1 (HIV-1) plays a major role in treatment resistance, from the development of vaccines to therapeutic drugs. In addressing the crux of the issue, various attempts to estimate the mutation rate of HIV-1 resulted in a large range of 10⁻⁵–10⁻³ errors/bp/cycle due to the use of different types of investigation methods. In this review, we discuss the different assay methods, their findings on the mutation rates of HIV-1 and how the locations of mutations can be further analyzed for their allosteric effects to allow for new inhibitor designs. Given that HIV is one of the fastest mutating viruses, it serves as a good model for the comprehensive study of viral mutations that can give rise to a more horizontal understanding towards overall viral drug resistance as well as emerging viral diseases.

Analysis of unusual and signature APOBEC-mutations in HIV-1 pol next-generation sequences

Introduction

At low mutation-detection thresholds, next generation sequencing (NGS) for HIV-1 genotypic resistance testing is susceptible to artifactual detection of mutations arising from PCR error and APOBEC-mediated G-to-A hypermutation.

Methods

We analyzed published HIV-1 pol Illumina NGS data to characterize the distribution of mutations at eight NGS mutation detection thresholds: 20%, 10%, 5%, 2%, 1%, 0.5%, 0.2%, and 0.1%. At each threshold, we determined proportions of amino acid mutations that were unusual (defined as having a prevalence <0.01% in HIV-1 group M sequences) or signature APOBEC mutations.

Results

Eight studies, containing 855 samples, in the NCBI Sequence Read Archive were analyzed. As detection thresholds were lowered, there was a progressive increase in the proportion of positions with usual and unusual mutations and in the proportion of all mutations that were unusual. The median proportion of positions with an unusual mutation increased gradually from 0% at the 20% threshold to 0.3% at the 1% threshold and then exponentially to 1.3% (0.5% threshold), 6.9% (0.2% threshold), and 23.2% (0.1% threshold). In two of three studies with available plasma HIV-1 RNA levels, the proportion of positions with unusual mutations was negatively associated with virus levels. Although the complete set of signature APOBEC mutations was much smaller than that of unusual mutations, the former outnumbered the latter in one-sixth of samples at the 0.5%, 1%, and 2% thresholds.

Conclusions

The marked increase in the proportion of positions with unusual mutations at thresholds below 1% and in samples with lower virus loads suggests that, at low thresholds, many unusual mutations are artifactual, reflecting PCR error or G-to-A hypermutation. Profiling the numbers of unusual and signature APOBEC pol mutations at different NGS mutation detection thresholds may be useful to avoid selecting a threshold that is too low and poses an unacceptable risk of identifying artifactual mutations.

Intact HIV Proviruses Persist in Children Seven to Nine Years after Initiation of Antiretroviral Therapy in the First Year of Life

In adults starting antiretroviral therapy (ART) during acute infection, 2% of proviruses that persist on ART are genetically intact by sequence analysis. In contrast, a recent report in children treated early failed to detect sequence-intact proviruses. In another cohort of children treated early, we sought to detect and characterize proviral sequences after 6 to 9 years on suppressive ART. Peripheral blood mononuclear cells (PBMC) from perinatally infected children from the Children with HIV Early antiRetroviral (CHER) study were analyzed. Nearly full-length proviral amplification and sequencing (NFL-PAS) were performed at one time point after 6 to 9 years on ART. Amplicons with large internal deletions were excluded (<9 kb). All amplicons of ≥9 kb were sequenced and analyzed through a bioinformatic pipeline to detect indels, frameshifts, or hypermutations that would render them defective. In eight children who started ART at a median age of 5.4 months (range, 2.0 to 11.1 months), 733 single NFL-PAS amplicons were generated. Of these, 534 (72.9%) had large internal deletions, 174 (23.7%) had hypermutations, 15 (1.4%) had small internal deletions, 3 (1.0%) had deletions in the packaging signal/major splice donor site, and 7 (1.0%) were sequence intact. These 7 intact sequences were from three children who initiated ART after 2.3 months of age, one of whom had two identical intact sequences, suggestive of a cell clone harboring a replication-competent provirus. No intact proviruses were detected in four children who initiated ART before 2.3 months of age. Rare, intact proviruses can be detected in children who initiate ART after 2.3 months of age and are probably, as in adults, maintained by clonal expansion of cells infected before ART initiation.IMPORTANCE There are limited data about the proviral landscape in children exhibiting long-term suppression after early treatment, particularly in Sub-Saharan Africa where HIV-1 subtype C predominates. Investigating the sequence-intact reservoir could provide insight on the mechanisms by which intact proviruses persist and inform ongoing cure efforts. Through nearly full-length proviral amplification and sequencing (NFL-PAS), we generated 733 NFL-PAS amplicons from eight children. We showed that rare, genetically intact proviruses could be detected in children who initiated ART after 2.3 months of age. The frequency of intact proviruses was lower (P < 0.05) than that reported for HIV subtype B-infected adults treated during early HIV infection. We show that cells harboring genetically intact HIV proviruses are rare in children exhibiting long-term suppression after early treatment and may require the processing of a large number of cells to assess reservoir size. This points to the need for efficient methods to accurately quantify latent reservoirs, particularly in pediatric studies where sample availability is limited.

Picomolar to Micromolar: Elucidating the Role of Distal Mutations in HIV-1 Protease in Conferring Drug Resistance

Drug resistance continues to be a growing global problem. The efficacy of small molecule inhibitors is threatened by pools of genetic diversity in all systems, including antibacterials, antifungals, cancer therapeutics, and antivirals. Resistant variants often include combinations of active site mutations and distal "secondary" mutations, which are thought to compensate for losses in enzymatic activity. HIV-1 protease is the ideal model system to investigate these combinations and underlying molecular mechanisms of resistance. Darunavir (DRV) binds wild-type (WT) HIV-1 protease with a potency of <5 pM, but we have identified a protease variant that loses potency to DRV 150 000-fold, with 11 mutations in and outside the active site. To elucidate the roles of these mutations in DRV resistance, we used a multidisciplinary approach, combining enzymatic assays, crystallography, and molecular dynamics simulations. Analysis of protease variants with 1, 2, 4, 8, 9, 10, and 11 mutations showed that the primary active site mutations caused ∼50-fold loss in potency (2 mutations), while distal mutations outside the active site further decreased DRV potency from 13 nM (8 mutations) to 0.76 μM (11 mutations). Crystal structures and simulations revealed that distal mutations induce subtle changes that are dynamically propagated through the protease. Our results reveal that changes remote from the active site directly and dramatically impact the potency of the inhibitor. Moreover, we find interdependent effects of mutations in conferring high levels of resistance. These mechanisms of resistance are likely applicable to many other quickly evolving drug targets, and the insights may have implications for the design of more robust inhibitors.

Evolutionary Dynamics in Structured Populations Under Strong Population Genetic Forces

In the long-term neutral equilibrium, high rates of migration between subpopulations result in little population differentiation. However, in the short-term, even very abundant migration may not be enough for subpopulations to equilibrate immediately. In this study, we investigate dynamical patterns of short-term population differentiation in adapting populations via stochastic and analytical modeling through time. We characterize a regime in which selection and migration interact to create non-monotonic patterns of population differentiation over time when migration is weaker than selection, but stronger than drift. We demonstrate how these patterns can be leveraged to estimate high migration rates using approximate Bayesian computation. We apply this approach to estimate fast migration in a rapidly adapting intra-host Simian-HIV population sampled from different anatomical locations. We find differences in estimated migration rates between different compartments, even though all are above [Formula: see text] = 1. This work demonstrates how studying demographic processes on the timescale of selective sweeps illuminates processes too fast to leave signatures on neutral timescales.

Circulation of multiple subtypes (A, G and CRFs 02_AG) of human immunodeficiency virus type 1 (HIV-1) in selected districts of Punjab province, Pakistan

Owing to consistent genetic mutation and recombination, various escape mutants and/or drug-resistant mutants of human immunodeficiency virus (HIV-1) are now emerging worldwide. Therefore, an understanding of the genetic characteristics of prevailing strains, particularly with regard to drug-resistance-associated substitutions, is essential for devising and implementing treatments and disease control interventions in endemic settings such as Pakistan. We processed a total of 130 plasma samples originating from HIV-treatment centers in selected districts of Punjab province, Pakistan. The samples were first screened using an HIV-1 Ag/Ab Combo test followed by amplification of the pol gene (1084 bp) from samples that were positive either for the antigen or for both the antigen and antibodies simultaneously. Screening revealed that a total of 45 samples were positive (34.62%; 95% CI: 26.99-43.13) for either antigen or both antigen and antibodies (n = 18, 40%; 95% CI: 27.02-54.55) or for antibodies alone (n = 27, 60%; 95% CI: 45.45-72.98). A largest number of positive samples was from the district of Lahore (n = 19/43, 44.18%; 95% CI: 30.44-58.9) followed by Faisalabad (n= 12/36, 33.33%; 95% CI: 20.21-49.66), Gujranwala (n = 05/23, 21.7%; 95% CI: 9.66-41.9) and Sargodha (n = 09/28, 32.1%; 95% CI: 17.93-50.66). The probability of occurrence of HIV infection was significantly associated with individuals having a history of injecting drug use (68.08%; OR = 11.15; 95% CI: 53.84-79.61, p = 0.0001). Phylogenetic analysis based on the pol gene showed that the sequences from this study clustered into three distinct clades representing recombinant form 02_AG (n = 14, 77.0%; 95% CI: 54.79-91.00), and subtypes A (n = 2, 11.1%; 95% CI: 3.1-32.8) and G (n = 2, 11.1%; 95% CI: 3.1-32.8). Although we screened 18 samples for drug-resistance-associated mutations, except for an accessory mutation (M46K) in the protease (PR) region in one subject, we found a lack of drug-resistance-associated substitutions in the PR region. On the other hand, we found two subjects (2/18) carrying a resistance-associated mutation (V106I) conferring a low level of resistance against non-nucleoside reverse transcriptase inhibitors. The present study shows that multiple subtypes of HIV-1 are present in the affected population. Continuous disease surveillance coupled with evaluation of drug resistance at higher resolution should be done in future studies.

Reviewing HIV-1 Gag Mutations in Protease Inhibitors Resistance: Insights for Possible Novel Gag Inhibitor Designs

HIV protease inhibitors against the viral protease are often hampered by drug resistance mutations in protease and in the viral substrate Gag. To overcome this drug resistance and inhibit viral maturation, targeting Gag alongside protease rather than targeting protease alone may be more efficient. In order to successfully inhibit Gag, understanding of its drug resistance mutations and the elicited structural changes on protease binding needs to be investigated. While mutations on Gag have already been mapped to protease inhibitor resistance, there remain many mutations, particularly the non-cleavage mutations, that are not characterized. Through structural studies to unravel how Gag mutations contributes to protease drug resistance synergistically, it is thus possible to glean insights to design novel Gag inhibitors. In this review, we discuss the structural role of both novel and previously reported Gag mutations in PI resistance, and how new Gag inhibitors can be designed.

Molecular Determinants of Epistasis in HIV-1 Protease: Elucidating the Interdependence of L89V and L90M Mutations in Resistance

Protease inhibitors have the highest potency among antiviral therapies against HIV-1 infections, yet the virus can evolve resistance. Darunavir (DRV), currently the most potent Food and Drug Administration-approved protease inhibitor, retains potency against single-site mutations. However, complex combinations of mutations can confer resistance to DRV. While the interdependence between mutations within HIV-1 protease is key for inhibitor potency, the molecular mechanisms that underlie this control remain largely unknown. In this study, we investigated the interdependence between the L89V and L90M mutations and their effects on DRV binding. These two mutations have been reported to be positively correlated with one another in HIV-1 patient-derived protease isolates, with the presence of one mutation making the probability of the occurrence of the second mutation more likely. The focus of our investigation is a patient-derived isolate, with 24 mutations that we call "KY"; this variant includes the L89V and L90M mutations. Three additional KY variants with back-mutations, KY(V89L), KY(M90L), and the KY(V89L/M90L) double mutation, were used to experimentally assess the individual and combined effects of these mutations on DRV inhibition and substrate processing. The enzymatic assays revealed that the KY(V89L) variant, with methionine at residue 90, is highly resistant, but its catalytic function is compromised. When a leucine to valine mutation at residue 89 is present simultaneously with the L90M mutation, a rescue of catalytic efficiency is observed. Molecular dynamics simulations of these DRV-bound protease variants reveal how the L90M mutation induces structural changes throughout the enzyme that undermine the binding interactions.

A scale-free analysis of the HIV-1 genome demonstrates multiple conserved regions of structural and functional importance

HIV-1 replicates via a low-fidelity polymerase with a high mutation rate; strong conservation of individual nucleotides is highly indicative of the presence of critical structural or functional properties. Identifying such conservation can reveal novel insights into viral behaviour. We analysed 3651 publicly available sequences for the presence of nucleic acid conservation beyond that required by amino acid constraints, using a novel scale-free method that identifies regions of outlying score together with a codon scoring algorithm. Sequences with outlying score were further analysed using an algorithm for producing local RNA folds whilst accounting for alignment properties. 11 different conserved regions were identified, some corresponding to well-known cis-acting functions of the HIV-1 genome but also others whose conservation has not previously been noted. We identify rational causes for many of these, including cis functions, possible additional reading frame usage, a plausible mechanism by which the central polypurine tract primes second-strand DNA synthesis and a conformational stabilising function of a region at the 5' end of env.

Management of Virologic Failure and HIV Drug Resistance

Approximately 20% of people with HIV in the United States prescribed antiretroviral therapy are not virally suppressed. Thus, optimal management of virologic failure has a critical role in the ability to improve viral suppression rates to improve long-term health outcomes for those infected and to achieve epidemic control. This article discusses the causes of virologic failure, the use of resistance testing to guide management after failure, interpretation and relevance of HIV drug resistance patterns, considerations for selection of second-line and salvage therapies, and management of virologic failure in special populations.