The causes and consequences of HIV evolution

Easing genomic surveillance: A comprehensive performance evaluation of long-read assemblers across multi-strain mixture data of HIV-1 and Other pathogenic viruses for constructing a user-friendly bioinformatic pipeline

Background

Determining the appropriate computational requirements and software performance is essential for efficient genomic surveillance. The lack of standardized benchmarking complicates software selection, especially with limited resources.

Methods

We developed a containerized benchmarking pipeline to evaluate seven long-read assemblers-Canu, GoldRush, MetaFlye, Strainline, HaploDMF, iGDA, and RVHaplo-for viral haplotype reconstruction, using both simulated and experimental Oxford Nanopore sequencing data of HIV-1 and other viruses. Benchmarking was conducted on three computational systems to assess each assembler's performance, utilizing QUAST and BLASTN for quality assessment.

Results

Our findings show that assembler choice significantly impacts assembly time, with CPU and memory usage having minimal effect. Assembler selection also influences the size of the contigs, with a minimum read length of 2,000 nucleotides required for quality assembly. A 4,000-nucleotide read length improves quality further. Canu was efficient among de novo assemblers but not suitable for multi-strain mixtures, while GoldRush produced only consensus assemblies. Strainline and MetaFlye were suitable for metagenomic sequencing data, with Strainline requiring high memory and MetaFlye operable on low-specification machines. Among reference-based assemblers, iGDA had high error rates, RVHaplo showed the best runtime and accuracy but became ineffective with similar sequences, and HaploDMF, utilizing machine learning, had fewer errors with a slightly longer runtime.

Conclusions

The HIV-64148 pipeline, containerized using Docker, facilitates easy deployment and offers flexibility to select from a range of assemblers to match computational systems or study requirements. This tool aids in genome assembly and provides valuable information on HIV-1 sequences, enhancing viral evolution monitoring and understanding.

Enforcing Temporal Consistency in Migration History Inference

In addition to undergoing evolution, members of biological populations may also migrate between locations. Examples include the spread of tumor cells from the primary tumor to distant metastases or the spread of pathogens from one host to another. One may represent migration histories by assigning a location label to each vertex of a given phylogenetic tree such that an edge connecting vertices with distinct locations represents a migration. Some biological populations undergo comigration, a phenomenon where multiple taxa from distinct lineages simultaneously comigrate from one location to another. In this work, we show that a previous problem statement for inferring migration histories that are parsimonious in terms of migrations and comigrations may lead to temporally inconsistent solutions. To remedy this deficiency, we introduce precise definitions of temporal consistency of comigrations in a phylogenetic tree, leading to three successive problems. First, we formulate the temporally consistent comigration problem to check if a set of comigrations is temporally consistent and provide a linear time algorithm for solving this problem. Second, we formulate the parsimonious consistent comigrations (PCC) problem, which aims to find comigrations given a location labeling of a phylogenetic tree. We show that PCC is NP-hard. Third, we formulate the parsimonious consistent comigration history (PCCH) problem, which infers the migration history given a phylogenetic tree and locations of its extant vertices only. We show that PCCH is NP-hard as well. On the positive side, we propose integer linear programming models to solve the PCC and PCCH problems. We demonstrate our algorithms on simulated and real data.

Multitasking Pharmacophores Support Cabotegravir-Based Long-Acting HIV Pre-Exposure Prophylaxis (PrEP)

Cabotegravir is an integrase strand transfer inhibitor (INSTI) for HIV treatment and prevention. Cabotegravir-based long-acting pre-exposure prophylaxis (PrEP) presents an emerging paradigm for infectious disease control. In this scheme, a combination of a high efficacy and low solubility of anti-infection drugs permits the establishment of a pharmaceutical firewall in HIV-vulnerable groups over a long period. Although the structure-activity-relationship (SAR) of cabotegravir as an INSTI is known, the structural determinants of its low solubility have not been identified. In this work, we have integrated multiple experimental and computational methods, namely X-ray diffraction, solid-state NMR (SSNMR) spectroscopy, solution NMR spectroscopy, automated fragmentation (AF)-QM/MM and density functional theory (DFT) calculations, to address this question. The molecular organization of cabotegravir in crystal lattice has been determined. The combination of very-fast magic-angle-sample-spinning (VF MAS) SSNMR and solution NMR, as supported by AF-QM/MM and DFT calculations, permits the identification of structural factors that contribute to the low aqueous solubility of cabotegravir. Our study reveals the multitasking nature of pharmacophores in cabotegravir, which controls the drug solubility and, meanwhile, the biological activity. By unraveling these function-defining molecular features, our work could inspire further development of long-acting HIV PrEP drugs.

Molecular dynamics simulations show how antibodies may rescue HIV-1 mutants incapable of infecting host cells

High mutation and replication rates of HIV-1 result in the continuous generation of variants, allowing it to adapt to changing host environments. Mutations often have deleterious effects, but variants carrying them are rapidly purged. Surprisingly, a particular variant incapable of entering host cells was found to be rescued by host antibodies targeting HIV-1. Understanding the molecular mechanism of this rescue is important to develop and improve antibody-based therapies. To unravel the underlying mechanisms, we performed fully atomistic molecular dynamics simulations of the HIV-1 gp41 trimer responsible for viral entry into host cells, its entry-deficient variant, and its complex with the rescuing antibody. We find that the Q563R mutation, which the entry-deficient variant carries, prevents the native conformation of the gp41 6-helix bundle required for entry and stabilizes an alternative conformation instead. This is the consequence of substantial changes in the secondary structure and interactions between the domains of gp41. Binding of the antibody F240 to gp41 reverses these changes and re-establishes the native conformation, resulting in rescue. To test the generality of this mechanism, we performed simulations with the entry-deficient L565A variant and antibody 3D6. We find that 3D6 binding was able to reverse structural and interaction changes introduced by the mutation and restore the native gp41 conformation. Viral variants may not only escape antibodies but be aided by them in their survival, potentially compromising antibody-based therapies, including vaccination and passive immunization. Our simulation framework could serve as a tool to assess the likelihood of such resistance against specific antibodies.Communicated by Ramaswamy H. SarmaCommunicated by Ramaswamy H. Sarma.

Effect of the SARS-CoV-2 Delta-associated G15U mutation on the s2m element dimerization and its interactions with miR-1307-3p

The s2m, a highly conserved 41-nt hairpin structure in the SARS-CoV-2 genome, serves as an attractive therapeutic target that may have important roles in the virus life cycle or interactions with the host. However, the conserved s2m in Delta SARS-CoV-2, a previously dominant variant characterized by high infectivity and disease severity, has received relatively less attention than that of the original SARS-CoV-2 virus. The focus of this work is to identify and define the s2m changes between Delta and SARS-CoV-2 and the subsequent impact of those changes upon the s2m dimerization and interactions with the host microRNA miR-1307-3p. Bioinformatics analysis of the GISAID database targeting the s2m element reveals a >99% correlation of a single nucleotide mutation at the 15th position (G15U) in Delta SARS-CoV-2. Based on 1H NMR spectroscopy assignments comparing the imino proton resonance region of s2m and the s2m G15U at 19°C, we show that the U15-A29 base pair closes, resulting in a stabilization of the upper stem without overall secondary structure deviation. Increased stability of the upper stem did not affect the chaperone activity of the viral N protein, as it was still able to convert the kissing dimers formed by s2m G15U into a stable duplex conformation, consistent with the s2m reference. However, we show that the s2m G15U mutation drastically impacts the binding of host miR-1307-3p. These findings demonstrate that the observed G15U mutation alters the secondary structure of s2m with subsequent impact on viral binding of host miR-1307-3p, with potential consequences on immune responses.

A Bayesian walker coupled with a computational workflow that generates the micro-evolution of SARS-CoV-2 and makes predictions of new mutations that can emerge

In this work, a Bayesian walker was constructed that generates mutations that are more prone as per UNIPROT variant data. The Bayesian walker was used to search the mutational landscape of Severe Acute Respiratory Syndrome Corona Virus 2 (SARS-CoV-2) and a computational workflow was followed to evaluate whether a particular mutation would satisfy natural selection's fitness criteria. For SARS-CoV-2, the empirical known fitness criteria derived from SARS-CoV-2 micro-evolution data is 3-fold criteria. Mutations that have emerged on the Spike protein of SARS-CoV-2 are ones that preserve the structural integrity of the Spike, retain, or increase infectivity and are Immune evasive as per literature reports. Based on the molecular mechanism of infectivity and Immune evasion of SARS-CoV-2, a molecular modelling workflow was adopted to investigate the evolutionary feasibility of the mutations generated by the walker, to check whether the mutations satisfy the 3-fold fitness criteria for SARS-CoV-2 micro-evolution. It was found that the walker (mutation generator) coupled with the computational workflow to evaluate the evolutionary fitness of the generated mutations, re-generated the mutants corresponding to the Alpha, Beta and Gamma variants of SARS-CoV-2 demonstrating the ability of the methodology to generate the micro-evolution of SARS-CoV-2. Having demonstrated the ability of the methodology to generate the micro-evolution of SARS-CoV-2, the computational methodology was used to make predictions for new mutants that could emerge.Communicated by Ramaswamy H. Sarma.

Genetic and antiretroviral drug resistance mutations analysis of reverse transcriptase and protease gene from Pakistani people living with HIV-1

BACKGROUND

Antiretroviral therapy (ART) effectiveness is compromised by the emergence of HIV drug resistance mutations (DRM) and can lead to the failure of ART. Apart from intrinsic viral factors, non-compliance with drugs and/or the use of sub-optimum therapy can lead to the emergence of DRMs. In Pakistan HIV currently exists as a concentrated epidemic, however, ART coverage is very low, and drug adherence is poor. ART is selected assuming without baseline genotyping. Pakistan has recently seen a rise in treatment failures, but the country's actual burden of DRM is still unknown. In this study, we perform the genetic and drug resistance analysis of the pol gene from Pakistani HIV-positive ART-naïve and ART-experienced individuals.

METHODS

In this study, HIV-1 pol was sequenced from 146 HIV-1 positive individuals, divided into ART-naïve (n = 37) and ART-experienced (n = 109). The sequences were also used to determine HIV-1 subtypes, the prevalence of DRM, and pol genetic variability.

RESULTS

DRM analysis identified numerous DRMs against reverse transcriptase inhibitors in both ART-naïve and ART-experienced groups, including a few that are classified as rare. Additionally, the ART-experienced group showed mutations associated with resistance to protease inhibitors. Genetic analysis showed negative selection pressure in both groups, but a higher rate of evolution in the ART-naïve group.

CONCLUSION

High prevalence of DRMs, especially against previous first-line treatment in ART- naïve and the accumulation of DRMs in ART-experienced groups is concerning and warrants that a more extensive DRM survey be carried out to inform first-line and second-line ART regimen recommendations.

HIV infection

The AIDS epidemic has been a global public health issue for more than 40 years and has resulted in ~40 million deaths. AIDS is caused by the retrovirus, HIV-1, which is transmitted via body fluids and secretions. After infection, the virus invades host cells by attaching to CD4 receptors and thereafter one of two major chemokine coreceptors, CCR5 or CXCR4, destroying the host cell, most often a T lymphocyte, as it replicates. If unchecked this can lead to an immune-deficient state and demise over a period of ~2-10 years. The discovery and global roll-out of rapid diagnostics and effective antiretroviral therapy led to a large reduction in mortality and morbidity and to an expanding group of individuals requiring lifelong viral suppressive therapy. Viral suppression eliminates sexual transmission of the virus and greatly improves health outcomes. HIV infection, although still stigmatized, is now a chronic and manageable condition. Ultimate epidemic control will require prevention and treatment to be made available, affordable and accessible for all. Furthermore, the focus should be heavily oriented towards long-term well-being, care for multimorbidity and good quality of life. Intense research efforts continue for therapeutic and/or preventive vaccines, novel immunotherapies and a cure.

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.

Inspecting the interaction between human immunodeficiency virus and the immune system through genetic turnover

Chronic infections of the human immunodeficiency virus (HIV) create a very complex coevolutionary process, where the virus tries to escape the continuously adapting host immune system. Quantitative details of this process are largely unknown and could help in disease treatment and vaccine development. Here we study a longitudinal dataset of ten HIV-infected people, where both the B-cell receptors and the virus are deeply sequenced. We focus on simple measures of turnover, which quantify how much the composition of the viral strains and the immune repertoire change between time points. At the single-patient level, the viral-host turnover rates do not show any statistically significant correlation, however, they correlate if one increases the amount of statistics by aggregating the information across patients. We identify an anti-correlation: large changes in the viral pool composition come with small changes in the B-cell receptor repertoire. This result seems to contradict the naïve expectation that when the virus mutates quickly, the immune repertoire needs to change to keep up. However, a simple model of antagonistically evolving populations can explain this signal. If it is sampled at intervals comparable with the sweep time, one population has had time to sweep while the second cannot start a counter-sweep, leading to the observed anti-correlation. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

HIV-1 Non-Group M Strains and ART

To eliminate HIV infection, there are several elements to take into account to limit transmission and break viral replication, such as epidemiological, preventive or therapeutic management. The UNAIDS goals of screening, treatment and efficacy should allow for this elimination if properly followed. For some infections, the difficulty is linked to the strong genetic divergence of the viruses, which can impact the virological and therapeutic management of patients. To completely eliminate HIV by 2030, we must therefore also be able to act on these atypical variants (HIV-1 non-group M) which are distinct from the group M pandemic viruses. While this diversity has had an impact on the efficacy of antiretroviral treatment in the past, recent data show that there is real hope of eliminating these forms, while maintaining vigilance and constant surveillance, so as not to allow more divergent and resistant forms to emerge. The aim of this work is therefore to share an update on the current knowledge on epidemiology, diagnosis and antiretroviral agent efficacy of HIV-1 non-M variants.

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic feline coronavirus phenotypes

Feline coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed feline enteric coronavirus [FECV]), with around 12 per cent developing into deadly feline infectious peritonitis (FIP; feline infectious peritonitis virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV- and FECV-specific signals of positive selection. We analyzed full-length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site (FCS) and the other within the fusion domain of Spike. We also found fifteen sites subject to positive selection associated with FIPV within Spike, eleven of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were fourteen sites (twelve novel sites) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 FCS and adjacent to C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype and included twenty-four positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV-wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that it is unlikely to be one singular 'switch' mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.

Estimating linkage disequilibrium and selection from allele frequency trajectories

Genetic sequences collected over time provide an exciting opportunity to study natural selection. In such studies, it is important to account for linkage disequilibrium to accurately measure selection and to distinguish between selection and other effects that can cause changes in allele frequencies, such as genetic hitchhiking or clonal interference. However, most high-throughput sequencing methods cannot directly measure linkage due to short-read lengths. Here we develop a simple method to estimate linkage disequilibrium from time-series allele frequencies. This reconstructed linkage information can then be combined with other inference methods to infer the fitness effects of individual mutations. Simulations show that our approach reliably outperforms inference that ignores linkage disequilibrium and, with sufficient sampling, performs similarly to inference using the true linkage information. We also introduce two regularization methods derived from random matrix theory that help to preserve its performance under limited sampling effects. Overall, our method enables the use of linkage-aware inference methods even for data sets where only allele frequency time series are available.

Effect of the SARS-CoV-2 Delta-associated G15U mutation on the s2m element dimerization and its interactions with miR-1307-3p

The stem loop 2 motif (s2m), a highly conserved 41-nucleotide hairpin structure in the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome, serves as an attractive therapeutic target that may have important roles in the virus life cycle or interactions with the host. However, the conserved s2m in Delta SARS-CoV-2, a previously dominant variant characterized by high infectivity and disease severity, has received relatively less attention than that of the original SARS-CoV-2 virus. The focus of this work is to identify and define the s2m changes between Delta and SARS-CoV-2 and subsequent impact of those changes upon the s2m dimerization and interactions with the host microRNA miR-1307-3p. Bioinformatics analysis of the GISAID database targeting the s2m element reveals a greater than 99% correlation of a single nucleotide mutation at the 15 th position (G15U) in Delta SARS-CoV-2. Based on 1 H NMR assignments comparing the imino proton resonance region of s2m and the G15U at 19°C, we find that the U15-A29 base pair closes resulting in a stabilization of the upper stem without overall secondary structure deviation. Increased stability of the upper stem did not affect the chaperone activity of the viral N protein, as it was still able to convert the kissing dimers formed by s2m G15U into a stable duplex conformation, consistent with the s2m reference. However, we find that the s2m G15U mutation drastically reduces the binding affinity of the host miR-1307-3p. These findings demonstrate that the observed G15U mutation alters the secondary structure of s2m with subsequent impact on viral binding of host miR-1307-3p, with potential consequences on the immune response.

A Repeat Pattern of Founder Events for SARS-CoV-2 Variants in Alaska

Alaska is a unique US state because of its large size, geographically disparate population density, and physical distance from the contiguous United States. Here, we describe a pattern of SARS-CoV-2 variant emergence across Alaska reflective of these differences. Using genomic data, we found that in Alaska, the Omicron sublineage BA.2.3 overtook BA.1.1 by the week of 27 February 2022, reaching 48.5% of sequenced cases. On the contrary, in the contiguous United States, BA.1.1 dominated cases for longer, eventually being displaced by BA.2 sublineages other than BA.2.3. BA.2.3 only reached a prevalence of 10.9% in the contiguous United States. Using phylogenetics, we found evidence of potential origins of the two major clades of BA.2.3 in Alaska and with logistic regression estimated how it emerged and spread throughout the state. The combined evidence is suggestive of founder events in Alaska and is reflective of how Alaska's unique dynamics influence the emergence of SARS-CoV-2 variants.

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic Feline Coronavirus phenotypes

Feline Coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed Feline Enteric Coronavirus [FECV]), with around 12% developing into deadly Feline Infectious Peritonitis (FIP; Feline Infectious Peritonitis Virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV and FECV specific signals of positive selection. We analyzed full length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site, and the other within the fusion domain of Spike. We also found 15 sites subject to positive selection associated with FIPV within Spike, 11 of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were 14 sites (12 novel) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 furin cleavage site and adjacent C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype, and included 24 positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that is unlikely to be one singular "switch" mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.

Mechanisms and Consequences of Genetic Variation in Hepatitis C Virus (HCV)

Chronic infection with hepatitis C virus (HCV) is an important contributor to the global incidence of liver diseases, including liver cirrhosis and hepatocellular carcinoma. Although common for single-stranded RNA viruses, HCV displays a remarkable high level of genetic diversity, produced primarily by the error-prone viral polymerase and host immune pressure. The high genetic heterogeneity of HCV has led to the evolution of several distinct genotypes and subtypes, with important consequences for pathogenesis, and clinical outcomes. Genetic variability constitutes an evasion mechanism against immune suppression, allowing the virus to evolve epitope escape mutants that avoid immune recognition. Thus, heterogeneity and variability of the HCV genome represent a great hindrance for the development of vaccines against HCV. In addition, the high genetic plasticity of HCV allows the virus to rapidly develop antiviral resistance mutations, leading to treatment failure and potentially representing a major hindrance for the cure of chronic HCV patients. In this chapter, we will present the central role that genetic diversity has in the viral life cycle and epidemiology of HCV. Incorporation errors and recombination, both the result of HCV polymerase activity, represent the main mechanisms of HCV evolution. The molecular details of both mechanisms have been only partially clarified and will be presented in the following sections. Finally, we will discuss the major consequences of HCV genetic diversity, namely its capacity to rapidly evolve antiviral and immunological escape variants that represent an important limitation for clearance of acute HCV, for treatment of chronic hepatitis C and for broadly protective vaccines.

A Computational Framework for Determining the Breadth of Antibodies Against Highly Mutable Pathogens

Highly mutable pathogens pose daunting challenges for antibody design. The usual criteria of high potency and specificity are often insufficient to design antibodies that provide long-lasting protection. This is due, in part, to the ability of the pathogen to rapidly acquire mutations that permit them to evade the designed antibodies. To overcome these limitations, design of antibodies with a larger neutralizing breadth can be pursued. Such broadly neutralizing antibodies (bnAbs) should remain targeted to a specific epitope, yet show robustness against pathogen mutability, thereby neutralizing a higher number of antigens. This is particularly important for highly mutable pathogens, like the influenza virus and the human immunodeficiency virus (HIV). The protocol describes a method for computing the "breadth" of a given antibody, an essential aspect of antibody design.

A large population sample of African HIV genomes from the 1980s reveals a reduction in subtype D over time associated with propensity for CXCR4 tropism

We present 109 near full-length HIV genomes amplified from blood serum samples obtained during early 1986 from across Uganda, which to our knowledge is the earliest and largest population sample from the initial phase of the HIV epidemic in Africa. Consensus sequences were made from paired-end Illumina reads with a target-capture approach to amplify HIV material following poor success with standard approaches. In comparisons with a smaller 'intermediate' genome dataset from 1998 to 1999 and a 'modern' genome dataset from 2007 to 2016, the proportion of subtype D was significantly higher initially, dropping from 67% (73/109), to 57% (26/46) to 17% (82/465) respectively (p < 0.0001). Subtype D has previously been shown to have a faster rate of disease progression than other subtypes in East African population studies, and to have a higher propensity to use the CXCR4 co-receptor ("X4 tropism"); associated with a decrease in time to AIDS. Here we find significant differences in predicted tropism between A1 and D subtypes in all three sample periods considered, which is particularly striking the 1986 sample: 66% (53/80) of subtype D env sequences were predicted to be X4 tropic compared with none of the 24 subtype A1. We also analysed the frequency of subtype in the envelope region of inter-subtype recombinants, and found that subtype A1 is over-represented in env, suggesting recombination and selection have acted to remove subtype D env from circulation. The reduction of subtype D frequency over three decades therefore appears to be a result of selective pressure against X4 tropism and its higher virulence. Lastly, we find a subtype D specific codon deletion at position 24 of the V3 loop, which may explain the higher propensity for subtype D to utilise X4 tropism.

Myomedin replicas of gp120 V3 loop glycan epitopes recognized by PGT121 and PGT126 antibodies as non-cognate antigens for stimulation of HIV-1 broadly neutralizing antibodies

Introduction

Imprinting broadly neutralizing antibody (bNAb) paratopes by shape complementary protein mimotopes represents a potential alternative for developing vaccine immunogens. This approach, designated as a Non-Cognate Ligand Strategy (NCLS), has recently been used for the identification of protein variants mimicking CD4 binding region epitope or membrane proximal external region (MPER) epitope of HIV-1 envelope (Env) glycoprotein. However, the potential of small binding proteins to mimic viral glycan-containing epitopes has not yet been verified.

Methods

In this work, we employed a highly complex combinatorial Myomedin scaffold library to identify variants recognizing paratopes of super candidate bNAbs, PGT121 and PGT126, specific for HIV-1 V3 loop epitopes.

Results

In the collection of Myomedins called MLD variants targeted to PGT121, three candidates competed with gp120 for binding to this bNAb in ELISA, thus suggesting an overlapping binding site and epitope-mimicking potential. Myomedins targeted to PGT126 designated MLB also provided variants that competed with gp120. Immunization of mice with MLB or MLD binders resulted in the production of anti-gp120 and -Env serum antibodies. Mouse hyper-immune sera elicited with MLB036, MLB041, MLB049, and MLD108 moderately neutralized 8-to-10 of 22 tested HIV-1-pseudotyped viruses of A, B, and C clades in vitro.

Discussion

Our data demonstrate that Myomedin-derived variants can mimic particular V3 glycan epitopes of prominent anti-HIV-1 bNAbs, ascertain the potential of particular glycans controlling neutralizing sensitivity of individual HIV-1 pseudoviruses, and represent promising prophylactic candidates for HIV-1 vaccine development.

Epidemic spreading under mutually independent intra- and inter-host pathogen evolution

The dynamics of epidemic spreading is often reduced to the single control parameter R₀ (reproduction-rate), whose value, above or below unity, determines the state of the contagion. If, however, the pathogen evolves as it spreads, R₀ may change over time, potentially leading to a mutation-driven spread, in which an initially sub-pandemic pathogen undergoes a breakthrough mutation. To predict the boundaries of this pandemic phase, we introduce here a modeling framework to couple the inter-host network spreading patterns with the intra-host evolutionary dynamics. We find that even in the extreme case when these two process are driven by mutually independent selection forces, mutations can still fundamentally alter the pandemic phase-diagram. The pandemic transitions, we show, are now shaped, not just by R₀, but also by the balance between the epidemic and the evolutionary timescales. If mutations are too slow, the pathogen prevalence decays prior to the appearance of a critical mutation. On the other hand, if mutations are too rapid, the pathogen evolution becomes volatile and, once again, it fails to spread. Between these two extremes, however, we identify a broad range of conditions in which an initially sub-pandemic pathogen can breakthrough to gain widespread prevalence.

A Single Early Introduction Governed Viral Diversity in the Second Wave of SARS-CoV-2 Epidemic in Hungary

Retrospective evaluation of past waves of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epidemic is key for designing optimal interventions against future waves and novel pandemics. Here, we report on analysing genome sequences of SARS-CoV-2 from the first two waves of the epidemic in 2020 in Hungary, mirroring a suppression and a mitigation strategy, respectively. Our analysis reveals that the two waves markedly differed in viral diversity and transmission patterns. Specifically, unlike in several European areas or in the USA, we have found no evidence for early introduction and cryptic transmission of the virus in the first wave of the pandemic in Hungary. Despite the introduction of multiple viral lineages, extensive community spread was prevented by a timely national lockdown in March 2020. In sharp contrast, the majority of the cases in the much larger second wave can be linked to a single transmission lineage of the pan-European B.1.160 variant. This lineage was introduced unexpectedly early, followed by a 2-month-long cryptic transmission before a soar of detected cases in September 2020. Epidemic analysis has revealed that the dominance of this lineage in the second wave was not associated with an intrinsic transmission advantage. This finding is further supported by the rapid replacement of B.1.160 by the alpha variant (B.1.1.7) that launched the third wave of the epidemic in February 2021. Overall, these results illustrate how the founder effect in combination with the cryptic transmission, instead of repeated international introductions or higher transmissibility, can govern viral diversity.

Evolution of Multiple Domains of the HIV-1 Envelope Glycoprotein during Coreceptor Switch with CCR5 Antagonist Therapy

HIV-1 uses CD4 as a receptor and chemokine receptors CCR5 and/or CXCR4 as coreceptors. CCR5 antagonists are a class of antiretrovirals used to inhibit viral entry. Phenotypic prediction algorithms such as Geno2Pheno are used to assess CCR5 antagonist eligibility, for which the V3 region is screened. However, there exist scenarios where the algorithm cannot give an accurate prediction of tropism. The current study examined coreceptor shift of HIV-1 from CCR5-tropic strains to CXCR4-tropic or dual-tropic strains among five subjects in a clinical trial of the CCR5 antagonist vicriviroc. Envelope gene amplicon libraries were constructed and subjected to next-generation sequencing, as well as single-clone sequencing and functional analyses. Approximately half of the amplified full-length single envelope-encoding clones had no significant activity for infection of cells expressing high levels of CD4 and CCR5 or CXCR4. Functional analysis of 9 to 21 individual infectious clones at baseline and at the time of VF were used to construct phylogenetic trees and sequence alignments. These studies confirmed that specific residues and the overall charge of the V3 loop were the major determinants of coreceptor use, in addition to specific residues in other domains of the envelope protein in V1/V2, V4, C3, and C4 domains that may be important for coreceptor shift. These results provide greater insight into the viral genetic determinants of coreceptor shift. IMPORTANCE This study is novel in combining single-genome sequence analysis and next-generation sequencing to characterize HIV-1 quasispecies. The work highlights the importance of mutants present at frequencies of 1% or less in development of drug resistance. This study highlights a critical role of specific amino acid substitutions outside V3 that contribute to coreceptor shift as well as important roles of the V1/V2, V4, C3, and C4 domain residues.

The Origin, Epidemiology, and Phylodynamics of Human Immunodeficiency Virus Type 1 CRF47_BF

CRF47_BF is a circulating recombinant form (CRF) of the human immunodeficiency virus type 1 (HIV-1), the etiological agent of AIDS. CRF47_BF represents one of 19 CRFx_BFs and has a geographic focus in Spain, where it was first identified in 2010. Since its discovery, CRF47_BF has expanded considerably in Spain, predominantly through heterosexual contact (∼56% of the infections). Little is known, however, about the origin and diversity of this CRF or its epidemiological correlates, as very few samples have been available so far. This study conducts a phylogenetic analysis with representatives of all CRFx_BF sequence types along with HIV-1 M Group subtypes to validate that the CRF47_BF sequences share a unique evolutionary history. The CRFx_BF sequences cluster into a single, not well supported, clade that includes their dominant parent subtypes (B and F). This clade also includes subtype D and excludes sub-subtype F2. However, the CRF47_BF sequences all share a most recent common ancestor. Further analysis of this clade couples CRF47_BF protease-reverse transcriptase sequences and epidemiological data from an additional 87 samples collected throughout Spain, as well as additional CRF47_BF database sequences from Brazil and Spain to investigate the origin and phylodynamics of CRF47_BF. The Spanish region with the highest proportion of CRF47_BF samples in the data set was the Basque Country (43.7%) with Navarre next highest at 19.5%. We include in our analysis epidemiological data on host sex, mode of transmission, time of collection, and geographic region. The phylodynamic analysis indicates that CRF47_BF originated in Brazil around 1999–2000 and spread to Spain from Brazil in 2002–2003. The virus spread rapidly throughout Spain with an increase in population size from 2011 to 2015 and leveling off more recently. Three strongly supported clusters associated with Spanish regions (Basque Country, Navarre, and Aragon), together comprising 60.8% of the Spanish samples, were identified, one of which was also associated with transmission among men who have sex with men. The expansion in Spain of CRF47_BF, together with that of other CRFs and subtype variants of South American origin, previously reported, reflects the increasing relationship between the South American and European HIV-1 epidemics.

Fc receptors and the diversity of antibody responses to HIV infection and vaccination

The development of an effective vaccine against HIV is desperately needed. The successive failures of HIV vaccine efficacy trials in recent decades have shown the difficulty of inducing an appropriate protective immune response to fight HIV. Different correlates of antibody parameters associated with a decreased risk of HIV-1 acquisition have been identified. However, these parameters are difficult to reproduce and improve, possibly because they have an intricate and combined action. Here, we describe the numerous antibody (Ab) functions associated with HIV-1 protection and report the interrelated parameters regulating their complex functions. Indeed, besides neutralizing and Fc-mediated activity, additional factors such as Ab type, concentration and kinetics of induction, and Fc-receptor expression and binding capacity also influence the protective effect conferred by Abs. As these parameters were described to be associated with ethnicity, age and sex, these additional factors must be considered for the development of an effective immune response. Therefore, future vaccine designs need to consider these multifaceted Ab functions together with the demographic attributes of the patient populations.

The phylodynamics of SARS-CoV-2 during 2020 in Finland

Background

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused millions of infections and fatalities globally since its emergence in late 2019. The virus was first detected in Finland in January 2020, after which it rapidly spread among the populace in spring. However, compared to other European nations, Finland has had a low incidence of SARS-CoV-2. To gain insight into the origins and turnover of SARS-CoV-2 lineages circulating in Finland in 2020, we investigated the phylogeographic and -dynamic history of the virus.

Methods

The origins of SARS-CoV-2 introductions were inferred via Travel-aware Bayesian time-measured phylogeographic analyses. Sequences for the analyses included virus genomes belonging to the B.1 lineage and with the D614G mutation from countries of likely origin, which were determined utilizing Google mobility data. We collected all available sequences from spring and fall peaks to study lineage dynamics.

Results

We observed rapid turnover among Finnish lineages during this period. Clade 20C became the most prevalent among sequenced cases and was replaced by other strains in fall 2020. Bayesian phylogeographic reconstructions suggested 42 independent introductions into Finland during spring 2020, mainly from Italy, Austria, and Spain.

Conclusions

A single introduction from Spain might have seeded one-third of cases in Finland during spring in 2020. The investigations of the original introductions of SARS-CoV-2 to Finland during the early stages of the pandemic and of the subsequent lineage dynamics could be utilized to assess the role of transboundary movements and the effects of early intervention and public health measures.

A repeat pattern of founder events for SARS-CoV-2 variants in Alaska

Alaska is a unique US state because of its large size, geographically disparate population density, and physical distance from the contiguous United States. Here, we describe a pattern of SARS-CoV-2 variant emergence across Alaska reflective of these differences. Using genomic data, we found that in Alaska the Omicron sublineage BA.2.3 overtook BA.1.1 by the week of 2022-02-27, reaching 48.5% of sequenced cases. On the contrary in the contiguous United States, BA.1.1 dominated cases for longer, eventually being displaced by BA.2 sublineages other than BA.2.3. BA.2.3 only reached a prevalence of 10.9% in the contiguous United States. Using phylogenetics, we found evidence of potential origins of the two major clades of BA.2.3 in Alaska and with logistic regression estimated how it emerged and spread throughout the state. The combined evidence is suggestive of founder events in Alaska and is reflective of how Alaska’s unique dynamics influence the emergence of SARS-CoV-2 variants.

Containing the spread of COVID-19 virus facing to its high mutation rate: approach to intervention using a nonspecific way of blocking its entry into the cells

Viruses have multiple mutation rates that are higher than any other member of the kingdom of life. This gives them the ability to evolve, even within the course of a single infection, and to evade multiple host defenses, thereby impacting pathogenesis. Additionally, there are also interplays between mutation and recombination and the high multiplicity of infection (MOI) that enhance viral adaptability and increase levels of recombination leading to complex and conflicting effects on genome selection, and the net results is difficult to predict. Recently, the outbreak of COVID-19 virus represents a pandemic threat that has been declared a public health emergency of international concern. Up to present, however, due to the high mutation rate of COVID-19 virus, there are no effective procedures to contain the spread of this virus across the globe. For such a purpose, there is then an urgent need to explore new approaches. As an opinion, the present approach emphasizes on (a) the use of a nonspecific way of blocking the entry of COVID-19 virus as well as its variants into the cells via a therapeutic biocompatible compound (ideally, "in a pill") targeting its spike (S) glycoprotein; and (b) the construction of expression vectors via the glycosyl-phosphatidylinositol, GPI, anchor for studying intermolecular interactions between the spike S of COVID-19 virus as well as its variants and the angiotensin-converting enzyme 2 (ACE2) of its host receptor for checking the efficacy of any therapeutic biocompatible compound of the nonspecific way of blocking. Such antiviral drug would be safer than the ACE1 and ACE2 inhibitors/angiotensin receptor blockers, and recombinant human ACE2 as well as nucleoside analogs or protease inhibitors used for fighting the spread of the virus inside the cells, and it would also be used as a universal one for any eventual future pandemic related to viruses, especially the RNA viruses with high mutation rates.

Social, Epidemiological, and Virological Characteristics from Peruvian Subjects Living with HIV-1/AIDS with Different Sexual Risk Behavior

HIV-1 genetic diversity and resistance profile might change according to the risky sexual behavior of the host. To show this, we recruited 134 individuals between the years 2015 and 2017 identified as transgender women sex workers (TWSW, n = 73) and Heterosexual Military Officers (HET-MO, n = 61). After obtaining informed consent, we collected a blood sample to perform the HIV genotyping, CD4 cell count, and viral load. We used bioinformatics approaches for detecting resistance mutations and recombination events. Epidemiological data showed that both groups reported sexually transmitted diseases and they were widespread among TWSW, especially syphilis and herpes virus (35.6%). Illegal drugs consumption was higher among TWSW (71.2%), whereas condom use was inconsistent for both HET-MO (57.4%) and TWSW (74.0%). TWSW showed the shortest time exposition to antiretroviral therapy (ART) (3.5 years) and the lowest access to ART (34.2%) that conducted treatment failure (>4 logs). HIV-1 sequences from TWSW and HET-MO were analyzed to determine the genetic diversity and antiretroviral drug resistance. Phylogeny analysis revealed 125 (93%) cases of subtype B, 01 subtype A (0.76%), 07 (5.30%) BF recombinants, and 01 (0.76%) AG recombinant. Also, TWSW showed a higher recombination index (9.5%, 7/73) than HET-MO (1.5%, 1/68). HET-MO only showed acquired resistance (26.23%, 16/61), whereas TWSW showed both acquired as transmitted resistance (9.59% for each). In conclusion, TWSW and HET-MO showed significant differences considering the epidemiological characteristics, genetic diversity, recombination events, and HIV resistance profile.

Immunoescape of HIV-1 in Env-EL9 CD8 + T cell response restricted by HLA-B*14:02 in a Non progressor who lost twenty-seven years of HIV-1 control

Background

Long-Term Non-Progressors (LTNPs) are untreated Human Immunodeficiency virus type 1 (HIV-1) infected individuals able to control disease progression for prolonged periods. However, the LTNPs status is temporary, as viral load increases followed by decreases in CD4 + T-cell counts. Control of HIV-1 infection in LTNPs viremic controllers, have been associated with effective immunodominant HIV-1 Gag-CD8 + T-cell responses restricted by protective HLA-B alleles. Individuals carrying HLA-B*14:02 control HIV-1 infection is related to an immunodominant Env-CD8 + T-cell response. Limited data are available on the contribution of HLA-B*14:02 CD8 + T -cells in LTNPs.

Results

In this study, we performed a virological and immunological detailed analysis of an HLA-B*14:02 LNTP individual that lost viral control (LVC) 27 years after HIV-1 diagnosis. We analysed viral evolution and immune escape in HLA-B*14:02 restricted CD8 + T -cell epitopes and identified viral evolution at the Env-EL9 epitope selecting the L592R mutation. By IFN-γ ELISpot and immune phenotype, we characterized HLA- B*14:02 HIV-1 CD8 + T cell responses targeting, Gag-DA9 and Env-EL9 epitopes before and after LVC. We observed an immunodominant response against the Env-EL9 epitope and a decreased of the CD8 T + cell response over time with LVC. Loss of Env-EL9 responses was concomitant with selecting K588R + L592R mutations at Env-EL9. Finally, we evaluated the impact of Env-EL9 escape mutations on HIV-1 infectivity and Env protein structure. The K588R + L592R escape variant was directly related to HIV-1 increase replicative capacity and stability of Env at the LVC.

Conclusions

These findings support the contribution of immunodominant Env-EL9 CD8 + T-cell responses and the imposition of immune escape variants with higher replicative capacity associated with LVC in this LNTP. These data highlight the importance of Env-EL9 specific-CD8 + T-cell responses restricted by the HLA-B*14:02 and brings new insights into understanding long-term HIV-1 control mediated by Env mediated CD8 + T-cell responses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12977-022-00591-7.

HIV-1-Specific Immunodominant T-Cell Responses Drive the Dynamics of HIV-1 Recombination Following Superinfection

A series of HIV-1 CRF01_AE/CRF07_BC recombinants were previously found to have emerged gradually in a superinfected patient (patient LNA819). However, the extent to which T-cell responses influenced the development of these recombinants after superinfection is unclear. In this study, we undertook a recombination structure analysis of the gag, pol, and nef genes from longitudinal samples of patient LNA819. A total of 9 pol and 5 nef CRF01_AE/CRF07_BC recombinants were detected. The quasispecies makeup and the composition of the pol and nef gene recombinants changed continuously, suggestive of continuous evolution in vivo. T-cell responses targeting peptides of the primary strain and the recombination regions were screened. The results showed that Pol-LY10, Pol-RY9, and Nef-GL9 were the immunodominant epitopes. Pol-LY10 overlapped with the recombination breakpoints in multiple recombinants. For the LY10 epitope, escape from T-cell responses was mediated by both recombination with a CRF07_BC insertion carrying the T467E/T472V variants and T467N/T472V mutations originating in the CRF01_AE strain. In pol recombinants R8 and R9, the recombination breakpoints were located ~23 amino acids upstream of the RY9 epitope. The appearance of new recombination breakpoints harboring a CRF07_BC insertion carrying a R984K variant was associated with escape from RY9-specific T-cell responses. Although the Nef-GL9 epitope was located either within or 10~11 amino acids downstream of the recombination breakpoints, no variant of this epitope was observed in the nef recombinants. Instead, a F85V mutation originating in the CRF01_AE strain was the main immune escape mechanism. Understanding the cellular immune pressure on recombination is critical for monitoring the new circulating recombinant forms of HIV and designing epitope-based vaccines. Vaccines targeting antigens that are less likely to escape immune pressure by recombination and/or mutation are likely to be of benefit to patients with HIV-1.

Threat, challenges, and preparedness for future pandemics: A descriptive review of phylogenetic analysis based predictions

For centuries the world has been confronted with many infectious diseases, with a potential to turn into a pandemic posing a constant threat to human lives. Some of these pandemics occurred due to the emergence of new disease or re-emergence of previously known diseases with a few mutations. In such scenarios their optimal prevention and control options were not adequately developed. Most of these diseases are highly contagious and for their timely control, knowledge about the pathogens and disease progression is the basic necessity. In this review, we have presented a documented chronology of the earlier pandemics, evolutionary analysis of the infectious disease with pandemic potential, the role of RNA, difficulties in controlling pandemics, and the likely pathogens that could trigger future pandemics. In this study, the evolutionary history of the pathogens was identified by carrying out phylogenetic analysis. The percentage similarity between different infectious diseases is critically analysed for the identification of their correlation using online sequence matcher tools. The Baltimore classification system was used for finding the genomic nature of the viruses. It was observed that most of the infectious pathogens rise from their animal hosts with some mutations in their genome composition. The phylogenetic tree shows that the single-stranded RNA diseases have a common origin and many of them are having high similarity percentage. The outcomes of this study will help in the identification of potential pathogens that can cause future pandemics. This information will be helpful in the development of early detection techniques, devising preventive mechanism to limit their spread, prophylactic measures, Infection control and therapeutic options, thereby, strengthening our approach towards global preparedness against future pandemics.

Limits to detecting epistasis in the fitness landscape of HIV

The rapid evolution of HIV is constrained by interactions between mutations which affect viral fitness. In this work, we explore the role of epistasis in determining the mutational fitness landscape of HIV for multiple drug target proteins, including Protease, Reverse Transcriptase, and Integrase. Epistatic interactions between residues modulate the mutation patterns involved in drug resistance, with unambiguous signatures of epistasis best seen in the comparison of the Potts model predicted and experimental HIV sequence “prevalences” expressed as higher-order marginals (beyond triplets) of the sequence probability distribution. In contrast, experimental measures of fitness such as viral replicative capacities generally probe fitness effects of point mutations in a single background, providing weak evidence for epistasis in viral systems. The detectable effects of epistasis are obscured by higher evolutionary conservation at sites. While double mutant cycles in principle, provide one of the best ways to probe epistatic interactions experimentally without reference to a particular background, we show that the analysis is complicated by the small dynamic range of measurements. Overall, we show that global pairwise interaction Potts models are necessary for predicting the mutational landscape of viral proteins.

Analysis of SARS-CoV-2 Variants From 24,181 Patients Exemplifies the Role of Globalization and Zoonosis in Pandemics

After the end of the first epidemic episode of SARS-CoV-2 infections, as cases began to rise again during the summer of 2020, we at IHU Méditerranée Infection in Marseille, France, intensified the genomic surveillance of SARS-CoV-2, and described the first viral variants. In this study, we compared the incidence curves of SARS-CoV-2-associated deaths in different countries and reported the classification of SARS-CoV-2 variants detected in our institute, as well as the kinetics and sources of the infections. We used mortality collected from a COVID-19 data repository for 221 countries. Viral variants were defined based on ≥5 hallmark mutations along the whole genome shared by ≥30 genomes. SARS-CoV-2 genotype was determined for 24,181 patients using next-generation genome and gene sequencing (in 47 and 11% of cases, respectively) or variant-specific qPCR (in 42% of cases). Sixteen variants were identified by analyzing viral genomes from 9,788 SARS-CoV-2-diagnosed patients. Our data show that since the first SARS-CoV-2 epidemic episode in Marseille, importation through travel from abroad was documented for seven of the new variants. In addition, for the B.1.160 variant of Pangolin classification (a.k.a. Marseille-4), we suspect transmission from farm minks. In conclusion, we observed that the successive epidemic peaks of SARS-CoV-2 infections are not linked to rebounds of viral genotypes that are already present but to newly introduced variants. We thus suggest that border control is the best mean of combating this type of introduction, and that intensive control of mink farms is also necessary to prevent the emergence of new variants generated in this animal reservoir.

High Recombination Rate of Hepatitis C Virus Revealed by a Green Fluorescent Protein Reconstitution Cell System

Genetic recombination is an important evolutionary mechanism for RNA viruses and can facilitate escape from immune and drug pressure. Recombinant hepatitis C virus (HCV) variants have rarely been detected in patients, suggesting that HCV has intrinsic low recombination rate. Recombination of HCV has been demonstrated in vitro between non-functional genomes, but its frequency and relevance for viral evolution and life cycle has not been clarified. We developed a cell-based assay to detect and quantify recombination between fully viable HCV genomes, using the reconstitution of green fluorescent protein (GFP) as a surrogate marker for recombination. Here, two GFP-expressing HCV genomes carrying different inactivating GFP mutations can produce a virus carrying a functional GFP by recombining within the GFP region. Generated constructs allowed quantification of recombination rates between markers spaced 603 and 553 nucleotides apart by flow cytometry and next-generation sequencing (NGS). Viral constructs showed comparable spread kinetics and reached similar infectivity titers in Huh7.5 cells, allowing their use in co-transfections and co-infections. Single-cycle co-transfection experiments, performed in CD81-deficient S29 cells, showed GFP expression in double-infected cells, demonstrating genome mixing and occurrence of recombination. Quantification of recombinant genomes by NGS revealed an average rate of 6.1 per cent, corresponding to 49 per cent of maximum detectable recombination (MDR). Experiments examining recombination during the full replication cycle of HCV, performed in Huh7.5 cells, demonstrated average recombination rates of 5.0 per cent (40.0 per cent MDR) and 3.6 per cent (28.8 per cent MDR) for markers spaced by 603 and 553 nucleotides, respectively, supporting a linear relationship between marker distance and recombination rates. First passage infections using recombinant virus supernatant resulted in comparable recombination rates of 5.9 per cent (47.2 per cent MDR) and 3.5 per cent (28.0 per cent MDR), respectively, for markers spaced by 603 and 553 nucleotides. We developed a functional cell-based assay that, to the best of our knowledge, allows for the first time detailed quantification of recombination rates using fully viable HCV constructs. Our data indicate that HCV recombines at high frequency between highly similar genomes and that the frequency of recombination increases with the distance between marker sites. These results have implication for our understanding of HCV evolution and emphasize the importance of recombination in the reassortment of mutations in the HCV genome.

Integration features of intact latent HIV-1 in CD4+ T cell clones contribute to viral persistence

Latent intact HIV-1 proviruses persist in a small subset of long-lived CD4+ T cells that can undergo clonal expansion in vivo. Expanded clones of CD4+ T cells dominate latent reservoirs in individuals on long-term antiretroviral therapy (ART) and represent a major barrier to HIV-1 cure. To determine how integration landscape might contribute to latency, we analyzed integration sites of near full length HIV-1 genomes from individuals on long-term ART, focusing on individuals whose reservoirs are highly clonal. We find that intact proviruses in expanded CD4+ T cell clones are preferentially integrated within Krüppel-associated box (KRAB) domain-containing zinc finger (ZNF) genes. ZNF genes are associated with heterochromatin in memory CD4+ T cells; nevertheless, they are expressed in these cells under steady-state conditions. In contrast to genes carrying unique integrations, ZNF genes carrying clonal intact integrations are down-regulated upon cellular activation. Together, the data suggest selected genomic sites, including ZNF genes, can be especially permissive for maintaining HIV-1 latency during memory CD4+ T cell expansion.

A Comprehensive Picture of the Structures, Energies, and Bonding in the Alanine Dimers

High level quantum mechanical computations and extensive stochastic searches of the potential energy surfaces of the Alanine dimers uncover rich and complex structural and interaction landscapes. A total of 416 strongly bound (up 13.4 kcal mol^-1 binding energies at the DLPNO-CCSD(T)/6-311++G(d,p) level corrected by the basis set superposition error and by the zero point vibrational energies over B3LYP-D3 geometries), close energy equilibrium structures were located, bonded via 32 specific types of intermolecular contacts including Y⋅⋅⋅H-X primary and Y⋅⋅⋅H-C secondary hydrogen bonds, H⋅⋅⋅H dihydrogen contacts, and non conventional anti-electrostatic Y δ - ⋯ X δ - interactions. The putative global minimum is triply degenerate, corresponding to the structure of the common dimer of a carboxylic acid. All quantum descriptors of chemical bonding point to a multitude of weak individual interactions within each dimer, whose cumulative effect results in large binding energies and in an attractive fluxional wall of non-covalent interactions in the interstitial region between the monomers.

Evolutionary analysis and lineage designation of SARS-CoV-2 genomes

The pandemic due to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the etiological agent of coronavirus disease 2019 (COVID-19), has caused immense global disruption. With the rapid accumulation of SARS-CoV-2 genome sequences, however, thousands of genomic variants of SARS-CoV-2 are now publicly available. To improve the tracing of the viral genomes' evolution during the development of the pandemic, we analyzed single nucleotide variants (SNVs) in 121,618 high-quality SARS-CoV-2 genomes. We divided these viral genomes into two major lineages (L and S) based on variants at sites 8782 and 28144, and further divided the L lineage into two major sublineages (L1 and L2) using SNVs at sites 3037, 14408, and 23403. Subsequently, we categorized them into 130 sublineages (37 in S, 35 in L1, and 58 in L2) based on marker SNVs at 201 additional genomic sites. This lineage/sublineage designation system has a hierarchical structure and reflects the relatedness among the subclades of the major lineages. We also provide a companion website (www.covid19evolution.net) that allows users to visualize sublineage information and upload their own SARS-CoV-2 genomes for sublineage classification. Finally, we discussed the possible roles of compensatory mutations and natural selection during SARS-CoV-2's evolution. These efforts will improve our understanding of the temporal and spatial dynamics of SARS-CoV-2's genome evolution.

Human B cell lineages associated with germinal centers following influenza vaccination are measurably evolving

The poor efficacy of seasonal influenza virus vaccines is often attributed to pre-existing immunity interfering with the persistence and maturation of vaccine-induced B cell responses. We previously showed that a subset of vaccine-induced B cell lineages are recruited into germinal centers (GCs) following vaccination, suggesting that affinity maturation of these lineages against vaccine antigens can occur. However, it remains to be determined whether seasonal influenza vaccination stimulates additional evolution of vaccine-specific lineages, and previous work has found no significant increase in somatic hypermutation among influenza-binding lineages sampled from the blood following seasonal vaccination in humans. Here, we investigate this issue using a phylogenetic test of measurable immunoglobulin sequence evolution. We first validate this test through simulations and survey measurable evolution across multiple conditions. We find significant heterogeneity in measurable B cell evolution across conditions, with enrichment in primary response conditions such as HIV infection and early childhood development. We then show that measurable evolution following influenza vaccination is highly compartmentalized: while lineages in the blood are rarely measurably evolving following influenza vaccination, lineages containing GC B cells are frequently measurably evolving. Many of these lineages appear to derive from memory B cells. We conclude from these findings that seasonal influenza virus vaccination can stimulate additional evolution of responding B cell lineages, and imply that the poor efficacy of seasonal influenza vaccination is not due to a complete inhibition of vaccine-specific B cell evolution.

HIV Proviral Burden, Genetic Diversity, and Dynamics in Viremic Controllers Who Subsequently Initiated Suppressive Antiretroviral Therapy

Curing HIV will require eliminating the reservoir of integrated, replication-competent proviruses that persist despite antiretroviral therapy (ART). Understanding the burden, genetic diversity, and longevity of persisting proviruses in diverse individuals with HIV is critical to this goal, but these characteristics remain understudied in some groups. Among them are viremic controllers—individuals who naturally suppress HIV to low levels but for whom therapy is nevertheless recommended. We reconstructed within-host HIV evolutionary histories from longitudinal single-genome amplified viral sequences in four viremic controllers who eventually initiated ART and used this information to characterize the age and diversity of proviruses persisting on therapy. We further leveraged these within-host proviral age distributions to estimate rates of proviral turnover prior to ART. This is an important yet understudied metric, since pre-ART proviral turnover dictates reservoir composition at ART initiation (and thereafter), which is when curative interventions, once developed, would be administered. Despite natural viremic control, all participants displayed significant within-host HIV evolution pretherapy, where overall on-ART proviral burden and diversity broadly reflected the extent of viral replication and diversity pre-ART. Consistent with recent studies of noncontrollers, the proviral pools of two participants were skewed toward sequences that integrated near ART initiation, suggesting dynamic proviral turnover during untreated infection. In contrast, proviruses recovered from the other two participants dated to time points that were more evenly spread throughout infection, suggesting slow or negligible proviral decay following deposition. HIV cure strategies will need to overcome within-host proviral diversity, even in individuals who naturally controlled HIV replication before therapy.

A Possible Role of Remdesivir and Plasma Therapy in the Selective Sweep and Emergence of New SARS-CoV-2 Variants

Since summer 2020, SARS-CoV-2 strains at the origin of the COVID-19 pandemic have suddenly been replaced by new SARS-CoV-2 variants, some of which are highly transmissible and spread at a high rate. These variants include the Marseille-4 lineage (Nextclade 20A.EU2) in Europe, the 20I/501Y.V1 variant first detected in the UK, the 20H/501Y.V2 variant first detected in South Africa, and the 20J/501Y.V3 variant first detected in Brazil. These variants are characterized by multiple mutations in the viral spike protein that is targeted by neutralizing antibodies elicited in response to infection or vaccine immunization. The usual coronavirus mutation rate through genetic drift alone cannot account for such rapid changes. Recent reports of the occurrence of such mutations in immunocompromised patients who received remdesivir and/or convalescent plasma or monoclonal antibodies to treat prolonged SARS-CoV-2 infections led us to hypothesize that experimental therapies that fail to cure the patients from COVID-19 could favor the emergence of immune escape SARS-CoV-2 variants. We review here the data that support this hypothesis and urge physicians and clinical trial promoters to systematically monitor viral mutations by whole-genome sequencing for patients who are administered these treatments.

HAPHPIPE: Haplotype Reconstruction and Phylodynamics for Deep Sequencing of Intrahost Viral Populations

Deep sequencing of viral populations using next-generation sequencing (NGS) offers opportunities to understand and investigate evolution, transmission dynamics, and population genetics. Currently, the standard practice for processing NGS data to study viral populations is to summarize all the observed sequences from a sample as a single consensus sequence, thus discarding valuable information about the intrahost viral molecular epidemiology. Furthermore, existing analytical pipelines may only analyze genomic regions involved in drug resistance, thus are not suited for full viral genome analysis. Here, we present HAPHPIPE, a HAplotype and PHylodynamics PIPEline for genome-wide assembly of viral consensus sequences and haplotypes. The HAPHPIPE protocol includes modules for quality trimming, error correction, de novo assembly, alignment, and haplotype reconstruction. The resulting consensus sequences, haplotypes, and alignments can be further analyzed using a variety of phylogenetic and population genetic software. HAPHPIPE is designed to provide users with a single pipeline to rapidly analyze sequences from viral populations generated from NGS platforms and provide quality output properly formatted for downstream evolutionary analyses.

Predicting HIV drug resistance using weighted machine learning method at target protein sequence-level

Acquired immune deficiency syndrome (AIDS) is a fatal disease caused by human immunodeficiency virus (HIV). Although 23 different drugs have been available, the treatment of AIDS remains challenging because the virus mutates very quickly which can lead to drug resistance. Therefore, predicting drug resistance before treatment is crucial for individual treatments. Here, based on HIV target protein sequence information, we analyzed 21-drug resistance caused by mutated residues using machine learning (ML) methods. To transform target sequences into numeric vectors, seven physicochemical properties were used, which can well represent the interacting characteristics of target proteins. Then, principal component analysis (PCA) method was adopted to reduce the feature dimensionality. Random forest (RF) and support vector machine (SVM) based on three different kernel functions, including linear, polynomial and radial basis function (RBF), were all employed. By comparisons, we found that RBF-based SVM method gives a comparative performance with RF model. Further, we added the weight information to RBF-based SVM method by four different weight evaluation methods of RF, eXtreme Gradient Boosting (XGB), CfsSubsetEval and ReliefFAttributeEval, respectively. Results show that the RF-weighted RBF-based SVM yield the superior performance and 13 out of 21 drug models provide the correlation coefficients (R²) over 0.8 and 3 of them are higher than 0.9. Finally, position-specific importance analysis indicates that most of the mutation residues with high RF weight scores are proved to be closely related with drug resistance, which has been revealed in previous reports. Overall, we can expect that this method can be a supplementary tool for predicting HIV drug resistance for newly discovered mutations. Here, based on HIV target protein sequence information, we analyzed 21-drug resistance caused by mutated residues using machine learning (ML) methods by fusing the weight information of different mutation positions.

SARS-CoV-2 Portrayed against HIV: Contrary Viral Strategies in Similar Disguise

SARS-CoV-2 and HIV are zoonotic viruses that rapidly reached pandemic scale, causing global losses and fear. The COVID-19 and AIDS pandemics ignited massive efforts worldwide to develop antiviral strategies and characterize viral architectures, biological and immunological properties, and clinical outcomes. Although both viruses have a comparable appearance as enveloped viruses with positive-stranded RNA and envelope spikes mediating cellular entry, the entry process, downstream biological and immunological pathways, clinical outcomes, and disease courses are strikingly different. This review provides a systemic comparison of both viruses' structural and functional characteristics, delineating their distinct strategies for efficient spread.

Evolutionary dynamics of group A and B respiratory syncytial virus in China, 2009-2018

Respiratory syncytial virus (RSV) is a major cause of acute respiratory tract infections in children and is a public health threat globally. To investigate the spatiotemporal dynamics of RSV evolution, we performed systematic phylogenetic analysis using all available sequences from the GenBank database, together with sequences from Shanghai, China. Both RSV-A and RSV-B appear to have originated in North America, with an inferred origin time of 1954.0 (1938.7-1967.6) and 1969.7 (1962.6-1975.5), respectively. BA-like strains of RSV-B, with a 60-nt insertion, and the ON1 strain of RSV-A, with a 72-nt insertion, emerged in 1997.6 (1996.2-1998.6) and 2010.1 (2009.1-2010.3), respectively. Since their origin, both genotypes have gradually replaced the former circulating genotypes to become the dominant strain. The population dynamic of RSV-A showed a seasonal epidemic pattern with obvious expansion in the periods of 2006-2007, 2010-2011, 2011-2012, and 2013-2014. Thirty fixed amino acid substitutions were identified during the divergence of NA4 from GA1 genotypes of RSV-A, and 13 were found during the divergence of SAB4 from GB1 of RSV-B. Importantly, ongoing evolution has occurred since the emergence of ON1, including four amino acid substitutions (I208L, E232G, T253K, and P314L). RSV-A genotypes GA5, NA4, NA1, and ON1 and RSV-B genotypes CB1, SAB4, BA-C, BA10, BA7, and BA9 were co-circulating in China from 2005 to 2015. In particular, RSV-A genotype ON1 was first detected in China in 2011, and it completely replaced GA2 to become the predominant strain after 2016. These data provide important insights into the evolution and epidemiology of RSV.

Genetic source completeness of HIV-1 circulating recombinant forms (CRFs) predicted by multi-label learning

MOTIVATION

Infection with strains of different subtypes and the subsequent crossover reading between the two strands of genomic RNAs by host cells' reverse transcriptase are the main causes of the vast HIV-1 sequence diversity. Such inter-subtype genomic recombinants can become circulating recombinant forms (CRFs) after widespread transmissions in a population. Complete prediction of all the subtype sources of a CRF strain is a complicated machine learning problem. It is also difficult to understand whether a strain is an emerging new subtype and if so, how to accurately identify the new components of the genetic source.

RESULTS

We introduce a multi-label learning algorithm for the complete prediction of multiple sources of a CRF sequence as well as the prediction of its chronological number. The prediction is strengthened by a voting of various multi-label learning methods to avoid biased decisions. In our steps, frequency and position features of the sequences are both extracted to capture signature patterns of pure subtypes and CRFs. The method was applied to 7185 HIV-1 sequences, comprising 5530 pure subtype sequences and 1655 CRF sequences. Results have demonstrated that the method can achieve very high accuracy (reaching 99%) in the prediction of the complete set of labels of HIV-1 recombinant forms. A few wrong predictions are actually incomplete predictions, very close to the complete set of genuine labels.

AVAILABILITY AND IMPLEMENTATION

https://github.com/Runbin-tang/The-source-of-HIV-CRFs-prediction.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Trends and strategies to combat viral infections: A review on FDA approved antiviral drugs

The infectious microscopic viruses invade living cells to reproduce themselves, and causes chronic infections such as HIV/AIDS, hepatitis B and C, flu, etc. in humans which may lead to death if not treated. Different strategies have been utilized to develop new and superior antiviral drugs to counter the viral infections. The FDA approval of HIV nucleoside reverse transcriptase inhibitor, zidovudine in 1987 boosted the development of antiviral agents against different viruses. Currently, there are a number of combination drugs developed against various viral infections to arrest the activity of same or different viral macromolecules at multiple stages of its life cycle; among which majority are targeted to interfere with the replication of viral genome. Besides these, other type of antiviral molecules includes entry inhibitors, integrase inhibitors, protease inhibitors, interferons, immunomodulators, etc. The antiviral drugs can be toxic to human cells, particularly in case of administration of combination drugs, and on the other hand viruses can grow resistant to the antiviral drugs. Furthermore, emergence of new viruses like Ebola, coronaviruses (SARS-CoV, SARS-CoV-2) emphasizes the need for more innovative strategies to develop better antiviral drugs to fight the existing and the emerging viral infections. Hence, we reviewed the strategic enhancements in developing antiviral drugs for the treatment of different viral infections over the years.

On the founder effect in COVID-19 outbreaks: how many infected travelers may have started them all?

How many incoming travelers (I0 at time 0, equivalent to the 'founders' in evolutionary genetics) infected with SARS-CoV-2 who visit or return to a region could have started the epidemic of that region? I0 would be informative about the initiation and progression of epidemics. To obtain I0 , we analyze the genetic divergence among viral populations of different regions. By applying the 'individual-output' model of genetic drift to the SARS-CoV-2 diversities, we obtain I0 < 10, which could have been achieved by one infected traveler in a long-distance flight. The conclusion is robust regardless of the source population, the continuation of inputs (It for t > 0) or the fitness of the variants. With such a tiny trickle of human movement igniting many outbreaks, the crucial stage of repressing an epidemic in any region should, therefore, be the very first sign of local contagion when positive cases first become identifiable. The implications of the highly 'portable' epidemics, including their early evolution prior to any outbreak, are explored in the companion study (Ruan et al., personal communication).

A Single Amino Acid Substitution at the HIV-1 Protease Termini Dimer Interface Significantly Reduces Viral Particles Processing Efficiency

The dimeric form of HIV-1 protease (PR) is required for its full proteolytic activity. The stability of the dimer primarily depends on the termini interface, with N-terminal residues 1-4 of one monomer encountering C-terminal residues 96-99 of another. We made an alanine substitution for valine 3 (V3) or leucine 97 (L97) at the termini dimer interface and tested their proteolytic activity. We found that an alanine substitution for L97 (PRL97A) completely inhibited the proteolytic activity of the PR. However, an alanine substitution for V3 (PRV3A) partially impaired the proteolytic activity. We then introduced two forced-dimerization systems involving nucleocapsid (NC) replacement or the addition of 1-2 leucine zippers to determine whether the proteolytic activity of dimer-defective PRs could be restored. We found that two forced-dimerization systems compensated for the defect in PRV3A, but not in PRL97A. This implies that PRV3A and PRL97A potentially impair the PR via different mechanisms or cause defects in PR activity to different extents. These novel findings will likely serve as a foundation for developing new PR inhibitors for treating drug-resistant HIV-1 infections in the future.

HIV-1 p24Gag adaptation to modern and archaic HLA-allele frequency differences in ethnic groups contributes to viral subtype diversification

Pathogen-driven selection and past interbreeding with archaic human lineages have resulted in differences in human leukocyte antigen (HLA)-allele frequencies between modern human populations. Whether or not this variation affects pathogen subtype diversification is unknown. Here we show a strong positive correlation between ethnic diversity in African countries and both human immunodeficiency virus (HIV)-1 p24gag and subtype diversity. We demonstrate that ethnic HLA-allele differences between populations have influenced HIV-1 subtype diversification as the virus adapted to escape common antiviral immune responses. The evolution of HIV Subtype B (HIV-B), which does not appear to be indigenous to Africa, is strongly affected by immune responses associated with Eurasian HLA variants acquired through adaptive introgression from Neanderthals and Denisovans. Furthermore, we show that the increasing and disproportionate number of HIV-infections among African Americans in the USA drive HIV-B evolution towards an Africa-centric HIV-1 state. Similar adaptation of other pathogens to HLA variants common in affected populations is likely.