Revisiting the recombinant history of HIV-1 group M with dynamic network community detection.

The prevailing abundance of full-length HIV type 1 (HIV-1) genome sequences provides an opportunity to revisit the standard model of HIV-1 group M (HIV-1/M) diversity that clusters genomes into largely nonrecombinant subtypes, which is not consistent with recent evidence of deep recombinant histories for simian immunodeficiency virus (SIV) and other HIV-1 groups. Here we develop an unsupervised nonparametric clustering approach, which does not rely on predefined nonrecombinant genomes, by adapting a community detection method developed for dynamic social network analysis. We show that this method (dynamic stochastic block model [DSBM]) attains a significantly lower mean error rate in detecting recombinant breakpoints in simulated data (quasibinomial generalized linear model (GLM), P<8×10−8), compared to other reference-free recombination detection programs (genetic algorithm for recombination detection [GARD], recombination detection program 4 [RDP4], and RDP5). When this method was applied to a representative sample of n = 525 actual HIV-1 genomes, we determined k = 29 as the optimal number of DSBM clusters and used change-point detection to estimate that at least 95% of these genomes are recombinant. Further, we identified both known and undocumented recombination hotspots in the HIV-1 genome and evidence of intersubtype recombination in HIV-1 subtype reference genomes. We propose that clusters generated by DSBM can provide an informative framework for HIV-1 classification.

Phylogenetic Reconstruction and Functional Characterization of the Ancestral Nef Protein of Primate Lentiviruses.

Nef is an accessory protein unique to the primate HIV-1, HIV-2, and SIV lentiviruses. During infection, Nef functions by interacting with multiple host proteins within infected cells to evade the immune response and enhance virion infectivity. Notably, Nef can counter immune regulators such as CD4 and MHC-I, as well as the SERINC5 restriction factor in infected cells. In this study, we generated a posterior sample of time-scaled phylogenies relating SIV and HIV Nef sequences, followed by reconstruction of ancestral sequences at the root and internal nodes of the sampled trees up to the HIV-1 Group M ancestor. Upon expression of the ancestral primate lentivirus Nef protein within CD4+ HeLa cells, flow cytometry analysis revealed that the primate lentivirus Nef ancestor robustly downregulated cell-surface SERINC5, yet only partially downregulated CD4 from the cell surface. Further analysis revealed that the Nef-mediated CD4 downregulation ability evolved gradually, while Nef-mediated SERINC5 downregulation was recovered abruptly in the HIV-1/M ancestor. Overall, this study provides a framework to reconstruct ancestral viral proteins and enable the functional characterization of these proteins to delineate how functions could have changed throughout evolutionary history.

Enhancement of CD4 Binding, Host Cell Entry, and Sensitivity to CD4bs Antibody Inhibition Conferred by a Natural but Rare Polymorphism in the HIV-1 Envelope.

A rare but natural polymorphism in the HIV-1 envelope (Env) glycoprotein, lysine at position 425 was selected as a mutation conferring resistance to maraviroc (MVC) in vitro. N425K has not been identified in HIV-infected individuals failing an MVC-based treatment. This study reports that the rare K425 polymorphism in an HIV-1 subtype A Env has increased affinity for CD4, resulting in faster host cell entry kinetics and the ability to scavenge for low cell surface expression of CD4 to mediate entry. Whereas the subtype A wild-type isolate-74 Env (N425) is inhibited by soluble (s) CD4, HIV-1 with K425 A74 Env shows enhanced infection and the ability to infect CCR5+ cells when pretreated with sCD4. Upon adding K425 or N425 HIV-1 to CD4+/CCR5+ cells along with RANTES/CCL3, only K425 HIV-1 was able to infect cells when CCR5 recycled/returned to the cell surface at 12 h post-treatment. These findings suggest that upon binding to CD4, K425 Env may maintain a stable State 2 "open" conformation capable of engaging CCR5 for entry. Only K425 was significantly more sensitivity than wild-type N425 A74 to inhibition by the CD4 binding site (bs) compound, BMS-806, the CD4bs antibody, VRC01 and N6, and the single-chain CD4i antibody, SCm9. K425 A74 was also capable of activating B cells expressing the VRC01 surface immunoglobulin. In summary, despite increased replicative fitness, we propose that K425 HIV-1 may be counterselected within infected individuals if K425 HIV-1 is rapidly eliminated by CD4bs-neutralizing antibodies. IMPORTANCE Typically, a natural amino acid polymorphism is found as the wild-type sequence in the HIV-1 population if it provides a selective advantage to the virus. The natural K425 polymorphism in HIV-1 Env results in higher host cell entry efficiency and greater replicative fitness by virtue of its high binding affinity to CD4. The studies presented herein suggest that the rare K425 HIV-1, compared to the common N425 HIV-1, may be more sensitive to inhibition by CD4bs-neutralizing antibodies (i.e., antibodies that bind to the CD4 binding pocket on the HIV-1 envelope glycoprotein). If CD4bs antibodies did emerge in an infected individual, the K425 HIV-1 may be hypersensitive to inhibition, and thus this K425 virus variant may be removed from the HIV-1 swarm despite its higher replication fitness. Studies are now underway to determine whether addition of the K425 polymorphism into the Envelope-based HIV-1 vaccines could enhance protective immunity.

High-level resistance to bictegravir and cabotegravir in subtype A- and D-infected HIV-1 patients failing raltegravir with multiple resistance mutations.

OBJECTIVES: The second-generation integrase strand transfer inhibitor (INSTI) bictegravir is becoming accessible in low- and middle-income countries (LMICs), and another INSTI, cabotegravir, has recently been approved as a long-acting injectable. Data on bictegravir and cabotegravir susceptibility in raltegravir-experienced HIV-1 subtype A- and D-infected patients carrying drug resistance mutations (DRMs) remain very scarce in LMICs. PATIENTS AND METHODS: HIV-1 integrase (IN)-recombinant viruses from eight patients failing raltegravir-based third-line therapy in Uganda were genotypically and phenotypically tested for susceptibility to bictegravir and cabotegravir. Ability of these viruses to integrate into human genomes was assessed in MT-4 cells. RESULTS: HIV-1 IN-recombinant viruses harbouring single primary mutations (N155H or Y143R/S) or in combination with secondary INSTI mutations (T97A, M50I, L74IM, E157Q, G163R or V151I) were susceptible to both bictegravir and cabotegravir. However, combinations of primary INSTI-resistance mutations such as E138A/G140A/G163R/Q148R or E138K/G140A/S147G/Q148K led to decreased susceptibility to both cabotegravir (fold change in EC50 values from 429 to 1000×) and bictegravir (60 to 100×), exhibiting a high degree of cross-resistance. However, these same IN-recombinant viruses showed impaired integration capacity (14% to 48%) relative to the WT HIV-1 NL4-3 strain in the absence of drug. CONCLUSIONS: Though not currently widely accessible in most LMICs, bictegravir and cabotegravir offer a valid alternative to HIV-infected individuals harbouring subtype A and D HIV-1 variants with reduced susceptibility to first-generation INSTIs but previous exposure to raltegravir may reduce efficacy, more so with cabotegravir.

Accumulation of integrase strand transfer inhibitor resistance mutations confers high-level resistance to dolutegravir in non-B subtype HIV-1 strains from patients failing raltegravir in Uganda.

BACKGROUND: Increasing first-line treatment failures in low- and middle-income countries (LMICs) have led to increased use of integrase strand transfer inhibitors (INSTIs) such as dolutegravir. However, HIV-1 susceptibility to INSTIs in LMICs, especially with previous raltegravir exposure, is poorly understood due to infrequent reporting of INSTI failures and testing for INSTI drug resistance mutations (DRMs). METHODS: A total of 51 non-subtype B HIV-1 infected patients failing third-line (raltegravir-based) therapy in Uganda were initially selected for the study. DRMs were detected using Sanger and deep sequencing. HIV integrase genes of 13 patients were cloned and replication capacities (RCs) and phenotypic susceptibilities to dolutegravir, raltegravir and elvitegravir were determined with TZM-bl cells. Spearman's correlation coefficient was used to determine cross-resistance between INSTIs. RESULTS: INSTI DRMs were detected in 47% of patients. HIV integrase-recombinant virus carrying one primary INSTI DRM (N155H or Y143R/S) was susceptible to dolutegravir but highly resistant to raltegravir and elvitegravir (>50-fold change). Two patients, one with E138A/G140A/Q148R/G163R and one with E138K/G140A/S147G/Q148K, displayed the highest reported resistance to raltegravir, elvitegravir and even dolutegravir. The former multi-DRM virus had WT RC whereas the latter had lower RCs than WT. CONCLUSIONS: In HIV-1 subtype A- and D-infected patients failing raltegravir and harbouring INSTI DRMs, there is high-level resistance to elvitegravir and raltegravir. More routine monitoring of INSTI treatment may be advised in LMICs, considering that multiple INSTI DRMs may have accumulated during prolonged exposure to raltegravir during virological failure, leading to high-level INSTI resistance, including dolutegravir resistance.

Effective and targeted latency reversal in CD4+ T cells from individuals on long term combined antiretroviral therapy initiated during chronic HIV-1 infection.

To date, an affordable, effective treatment for an HIV-1 cure remains only a concept with most "latency reversal" agents (LRAs) lacking specificity for the latent HIV-1 reservoir and failing in early clinical trials. We assessed HIV-1 latency reversal using a multivalent HIV-1-derived virus-like particle (HLP) to treat samples from 32 people living with HIV-1 (PLWH) in Uganda, US and Canada who initiated combined antiretroviral therapy (cART) during chronic infection. Even after 5-20 years on stable cART, HLP could target CD4+ T cells harbouring latent HIV-1 reservoir resulting in 100-fold more HIV-1 release into culture supernatant than by common recall antigens, and 1000-fold more than by chemotherapeutic LRAs. HLP induced release of a divergent and replication-competent HIV-1 population from PLWH on cART. These findings suggest HLP provides a targeted approach to reactivate the majority of latent HIV-1 proviruses among individuals infected with HIV-1.

Tracking SARS-CoV-2 variants of concern in wastewater: an assessment of nine computational tools using simulated genomic data.

Wastewater-based surveillance (WBS) is an important epidemiological and public health tool for tracking pathogens across the scale of a building, neighbourhood, city, or region. WBS gained widespread adoption globally during the SARS-CoV-2 pandemic for estimating community infection levels by qPCR. Sequencing pathogen genes or genomes from wastewater adds information about pathogen genetic diversity, which can be used to identify viral lineages (including variants of concern) that are circulating in a local population. Capturing the genetic diversity by WBS sequencing is not trivial, as wastewater samples often contain a diverse mixture of viral lineages with real mutations and sequencing errors, which must be deconvoluted computationally from short sequencing reads. In this study we assess nine different computational tools that have recently been developed to address this challenge. We simulated 100 wastewater sequence samples consisting of SARS-CoV-2 BA.1, BA.2, and Delta lineages, in various mixtures, as well as a Delta-Omicron recombinant and a synthetic 'novel' lineage. Most tools performed well in identifying the true lineages present and estimating their relative abundances and were generally robust to variation in sequencing depth and read length. While many tools identified lineages present down to 1 % frequency, results were more reliable above a 5 % threshold. The presence of an unknown synthetic lineage, which represents an unclassified SARS-CoV-2 lineage, increases the error in relative abundance estimates of other lineages, but the magnitude of this effect was small for most tools. The tools also varied in how they labelled novel synthetic lineages and recombinants. While our simulated dataset represents just one of many possible use cases for these methods, we hope it helps users understand potential sources of error or bias in wastewater sequencing analysis and to appreciate the commonalities and differences across methods.

The Canadian VirusSeq Data Portal & Duotang: open resources for SARS-CoV-2 viral sequences and genomic epidemiology.

The COVID-19 pandemic led to a large global effort to sequence SARS-CoV-2 genomes from patient samples to track viral evolution and inform public health response. Millions of SARS-CoV-2 genome sequences have been deposited in global public repositories. The Canadian COVID-19 Genomics Network (CanCOGeN - VirusSeq), a consortium tasked with coordinating expanded sequencing of SARS-CoV-2 genomes across Canada early in the pandemic, created the Canadian VirusSeq Data Portal, with associated data pipelines and procedures, to support these efforts. The goal of VirusSeq was to allow open access to Canadian SARS-CoV-2 genomic sequences and enhanced, standardized contextual data that were unavailable in other repositories and that meet FAIR standards (Findable, Accessible, Interoperable and Reusable). In addition, the Portal data submission pipeline contains data quality checking procedures and appropriate acknowledgement of data generators that encourages collaboration. From inception to execution, the portal was developed with a conscientious focus on strong data governance principles and practices. Extensive efforts ensured a commitment to Canadian privacy laws, data security standards, and organizational processes. This Portal has been coupled with other resources like Viral AI and was further leveraged by the Coronavirus Variants Rapid Response Network (CoVaRR-Net) to produce a suite of continually updated analytical tools and notebooks. Here we highlight this Portal, including its contextual data not available elsewhere, and the 'Duotang', a web platform that presents key genomic epidemiology and modeling analyses on circulating and emerging SARS-CoV-2 variants in Canada. Duotang presents dynamic changes in variant composition of SARS-CoV-2 in Canada and by province, estimates variant growth, and displays complementary interactive visualizations, with a text overview of the current situation. The VirusSeq Data Portal and Duotang resources, alongside additional analyses and resources computed from the Portal (COVID-MVP, CoVizu), are all open-source and freely available. Together, they provide an updated picture of SARS-CoV-2 evolution to spur scientific discussions, inform public discourse, and support communication with and within public health authorities. They also serve as a framework for other jurisdictions interested in open, collaborative sequence data sharing and analyses.

Early Warning Measurement of SARS-CoV-2 Variants of Concern in Wastewaters by Mass Spectrometry.

Wastewater surveillance has rapidly emerged as an early warning tool to track COVID-19. However, the early warning measurement of new SARS-CoV-2 variants of concern (VOCs) in wastewaters remains a major challenge. We herein report a rapid analytical strategy for quantitative measurement of VOCs, which couples nested polymerase chain reaction and liquid chromatography-mass spectrometry (nPCR-LC-MS). This method showed a greater selectivity than the current allele-specific quantitative PCR (AS-qPCR) for tracking new VOC and allowed the detection of multiple signature mutations in a single measurement. By measuring the Omicron variant in wastewaters across nine Ontario wastewater treatment plants serving over a three million population, the nPCR-LC-MS method demonstrated a better quantification accuracy than next-generation sequencing (NGS), particularly at the early stage of community spreading of Omicron. This work addresses a major challenge for current SARS-CoV-2 wastewater surveillance by rapidly and accurately measuring VOCs in wastewaters for early warning.

Tuning intrinsic disorder predictors for virus proteins.

Many virus-encoded proteins have intrinsically disordered regions that lack a stable, folded three-dimensional structure. These disordered proteins often play important functional roles in virus replication, such as down-regulating host defense mechanisms. With the widespread availability of next-generation sequencing, the number of new virus genomes with predicted open reading frames is rapidly outpacing our capacity for directly characterizing protein structures through crystallography. Hence, computational methods for structural prediction play an important role. A large number of predictors focus on the problem of classifying residues into ordered and disordered regions, and these methods tend to be validated on a diverse training set of proteins from eukaryotes, prokaryotes, and viruses. In this study, we investigate whether some predictors outperform others in the context of virus proteins and compared our findings with data from non-viral proteins. We evaluate the prediction accuracy of 21 methods, many of which are only available as web applications, on a curated set of 126 proteins encoded by viruses. Furthermore, we apply a random forest classifier to these predictor outputs. Based on cross-validation experiments, this ensemble approach confers a substantial improvement in accuracy, e.g., a mean 36 per cent gain in Matthews correlation coefficient. Lastly, we apply the random forest predictor to severe acute respiratory syndrome coronavirus 2 ORF6, an accessory gene that encodes a short (61 AA) and moderately disordered protein that inhibits the host innate immune response. We show that disorder prediction methods perform differently for viral and non-viral proteins, and that an ensemble approach can yield more robust and accurate predictions.

CoVizu: Rapid analysis and visualization of the global diversity of SARS-CoV-2 genomes.

Phylogenetics has played a pivotal role in the genomic epidemiology of severe acute respiratory syndrome coronavirus 2, such as tracking the emergence and global spread of variants and scientific communication. However, the rapid accumulation of genomic data from around the world-with over two million genomes currently available in the Global Initiative on Sharing All Influenza Data database-is testing the limits of standard phylogenetic methods. Here, we describe a new approach to rapidly analyze and visualize large numbers of SARS-CoV-2 genomes. Using Python, genomes are filtered for problematic sites, incomplete coverage, and excessive divergence from a strict molecular clock. All differences from the reference genome, including indels, are extracted using minimap2 and compactly stored as a set of features for each genome. For each Pango lineage (https://cov-lineages.org), we collapse genomes with identical features into 'variants', generate 100 bootstrap samples of the feature set union to generate weights, and compute the symmetric differences between the weighted feature sets for every pair of variants. The resulting distance matrices are used to generate neighbor-joining trees in RapidNJ that are converted into a majority-rule consensus tree for each lineage. Branches with support values below 50 per cent or mean lengths below 0.5 differences are collapsed, and tip labels on affected branches are mapped to internal nodes as directly sampled ancestral variants. Currently, we process about 2 million genomes in approximately 9 h on 52 cores. The resulting trees are visualized using the JavaScript framework D3.js as 'beadplots', in which variants are represented by horizontal line segments, annotated with beads representing samples by collection date. Variants are linked by vertical edges to represent branches in the consensus tree. These visualizations are published at https://filogeneti.ca/CoVizu. All source code was released under an MIT license at https://github.com/PoonLab/covizu.

Evidence for a recombinant origin of HIV-1 Group M from genomic variation.

Reconstructing the early dynamics of the HIV-1 pandemic can provide crucial insights into the socioeconomic drivers of emerging infectious diseases in human populations, including the roles of urbanization and transportation networks. Current evidence indicates that the global pandemic comprising almost entirely of HIV-1/M originated around the 1920s in central Africa. However, these estimates are based on molecular clock estimates that are assumed to apply uniformly across the virus genome. There is growing evidence that recombination has played a significant role in the early history of the HIV-1 pandemic, such that different regions of the HIV-1 genome have different evolutionary histories. In this study, we have conducted a dated-tip analysis of all near full-length HIV-1/M genome sequences that were published in the GenBank database. We used a sliding window approach similar to the 'bootscanning' method for detecting breakpoints in inter-subtype recombinant sequences. We found evidence of substantial variation in estimated root dates among windows, with an estimated mean time to the most recent common ancestor of 1922. Estimates were significantly autocorrelated, which was more consistent with an early recombination event than with stochastic error variation in phylogenetic reconstruction and dating analyses. A piecewise regression analysis supported the existence of at least one recombination breakpoint in the HIV-1/M genome with interval-specific means around 1929 and 1913, respectively. This analysis demonstrates that a sliding window approach can accommodate early recombination events outside the established nomenclature of HIV-1/M subtypes, although it is difficult to incorporate the earliest available samples due to their limited genome coverage.

Adaptive HIV-1 evolutionary trajectories are constrained by protein stability.

Despite the use of combination antiretroviral drugs for the treatment of HIV-1 infection, the emergence of drug resistance remains a problem. Resistance may be conferred either by a single mutation or a concerted set of mutations. The involvement of multiple mutations can arise due to interactions between sites in the amino acid sequence as a consequence of the need to maintain protein structure. To better understand the nature of such epistatic interactions, we reconstructed the ancestral sequences of HIV-1's Pol protein, and traced the evolutionary trajectories leading to mutations associated with drug resistance. Using contemporary and ancestral sequences we modelled the effects of mutations (i.e. amino acid replacements) on protein structure to understand the functional effects of residue changes. Although the majority of resistance-associated sequences tend to destabilise the protein structure, we find there is a general tendency for protein stability to decrease across HIV-1's evolutionary history. That a similar pattern is observed in the non-drug resistance lineages indicates that non-resistant mutations, for example, associated with escape from the immune response, also impacts on protein stability. Maintenance of optimal protein structure therefore represents a major constraining factor to the evolution of HIV-1.

Ebolavirus is evolving but not changing: No evidence for functional change in EBOV from 1976 to the 2014 outbreak.

The 2014 epidemic of Ebola virus disease (EVD) has had a devastating impact in West Africa. Sequencing of ebolavirus (EBOV) from infected individuals has revealed extensive genetic variation, leading to speculation that the virus may be adapting to humans, accounting for the scale of the 2014 outbreak. We computationally analyze the variation associated with all EVD outbreaks, and find none of the amino acid replacements lead to identifiable functional changes. These changes have minimal effect on protein structure, being neither stabilizing nor destabilizing, are not found in regions of the proteins associated with known functions and tend to cluster in poorly constrained regions of proteins, specifically intrinsically disordered regions. We find no evidence that the difference between the current and previous outbreaks is due to evolutionary changes associated with transmission to humans. Instead, epidemiological factors are likely to be responsible for the unprecedented spread of EVD.