Mutation rate of SARS-CoV-2 and emergence of mutators during experimental evolution

Combining Short- and Long-Read Sequencing Technologies to Identify SARS-CoV-2 Variants in Wastewater

During the COVID-19 pandemic, the monitoring of SARS-CoV-2 RNA in wastewater was used to track the evolution and emergence of variant lineages and gauge infection levels in the community, informing appropriate public health responses without relying solely on clinical testing. As more sublineages were discovered, it increased the difficulty in identifying distinct variants in a mixed population sample, particularly those without a known lineage. Here, we compare the sequencing technology from Illumina and from Oxford Nanopore Technologies, in order to determine their efficacy at detecting variants of differing abundance, using 248 wastewater samples from various Quebec and Ontario cities. Our study used two analytical approaches to identify the main variants in the samples: the presence of signature and marker mutations and the co-occurrence of signature mutations within the same amplicon. We observed that each sequencing method detected certain variants at different frequencies as each method preferentially detects mutations of distinct variants. Illumina sequencing detected more mutations with a predominant lineage that is in low abundance across the population or unknown for that time period, while Nanopore sequencing had a higher detection rate of mutations that are predominantly found in the high abundance B.1.1.7 (Alpha) lineage as well as a higher sequencing rate of co-occurring mutations in the same amplicon. We present a workflow that integrates short-read and long-read sequencing to improve the detection of SARS-CoV-2 variant lineages in mixed population samples, such as wastewater.

A Brighton Collaboration standardized template with key considerations for a benefit-risk assessment for the Comirnaty COVID-19 mRNA vaccine

The Brighton Collaboration Benefit-Risk Assessment of VAccines by TechnolOgy (BRAVATO) Working Group evaluates the safety and other key features of new platform technology vaccines, including nucleic acid (RNA and DNA) vaccines. This manuscript uses the BRAVATO template to report the key considerations for a benefit-risk assessment of the coronavirus disease 2019 (COVID-19) mRNA-based vaccine BNT162b2 (Comirnaty®, or Pfizer-BioNTech COVID-19 vaccine) including the subsequent Original/Omicron BA.1, Original/Omicron BA.4-5 and Omicron XBB.1.5 variant-adapted vaccines developed by BioNTech and Pfizer to protect against COVID-19 caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Initial Emergency Use Authorizations or conditional Marketing Authorizations for the original BNT162b2 vaccine were granted based upon a favorable benefit-risk assessment taking into account clinical safety, immunogenicity, and efficacy data, which was subsequently reconfirmed for younger age groups, and by real world evidence data. In addition, the favorable benefit-risk assessment was maintained for the bivalent vaccines, developed against newly arising SARS-CoV-2 variants, with accumulating clinical trial data.

Worldwide SARS-CoV-2 Omicron variant infection: Emerging sub-variants and future vaccination perspectives

The coronavirus disease 2019 (COVID-19) pandemic has sparked widespread concern globally, particularly with the Omicron variant and its sub-lineages emerging as the predominant cause of infection for nearly two years. Taiwan's successful containment of COVID-19, underscored by broad vaccine coverage, the utilization of anti-viral therapeutics, and timely response strategies, has resulted in reduced excess mortality. Moreover, there is a crucial need for a phased exit strategy, balancing efforts to curtail disease transmission with the mitigation of socioeconomic impacts from rigorous measures. In this review, we examined the evolution and the epidemiological landscape of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron sub-variants in Taiwan as well as other countries of the world. We also critically evaluated the effectiveness of COVID-19 vaccines against various SARS-CoV-2 variants. Additionally, we addressed the advantages of heterologous immunization strategies, fluctuations in neutralizing antibody titers, and complexities in establishing protective correlates among swiftly mutating viral variants.

Rapid intra-host diversification and evolution of SARS-CoV-2 in advanced HIV infection

Previous studies have linked the evolution of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) genetic variants to persistent infections in people with immunocompromising conditions, but the processes responsible for these observations are incompletely understood. Here we use high-throughput, single-genome amplification and sequencing (HT-SGS) to sequence SARS-CoV-2 spike genes from people with HIV (PWH, n = 22) and people without HIV (PWOH, n = 25). In PWOH and PWH with CD4 T cell counts (i.e., CD4 counts) ≥ 200 cells/μL, we find that most SARS-CoV-2 genomes sampled in each person share one spike sequence. By contrast, in people with advanced HIV infection (i.e., CD4 counts < 200 cells/μL), HT-SGS reveals a median of 46 distinct linked groupings of spike mutations per person. Elevated intra-host spike diversity in people with advanced HIV infection is detected immediately after COVID-19 symptom onset, and early intra-host spike diversity predicts SARS-CoV-2 shedding duration among PWH. Analysis of longitudinal timepoints reveals rapid fluctuations in spike sequence populations, replacement of founder sequences by groups of new haplotypes, and positive selection at functionally important residues. These findings demonstrate remarkable intra-host genetic diversity of SARS-CoV-2 in advanced HIV infection and suggest that adaptive intra-host SARS-CoV-2 evolution in this setting may contribute to the emergence of new variants of concern.

A novel immunoassay technique using principal component analysis for enhanced detection of emerging viral variants

Rapid diagnostics are critical infectious disease tools that are designed to detect a known biomarker using antibodies specific to that biomarker. However, a way to detect unknown disease variants has not yet been achieved in a paper test format. We describe here a route to make an adaptable paper immunoassay that can detect an unknown biomarker, demonstrating it on SARS-CoV-2 variants. The immunoassay repurposes cross reactive antibodies raised against the alpha variant. Gold nanoparticles of two different colors conjugated to two different antibodies create a colorimetric signal, and machine learning of the resulting colorimetric pattern is used to train the assay to discriminate between variants of alpha and Omicron BA.5. By using principal component analysis, the colorimetric test patterns can pick up and discriminate an unknown variant that it has not encountered before, Omicron BA.1. The test has an accuracy of 100% and a potential calculated discriminatory power of 900. We show that it can be used adaptively and that it can be used to pick up emerging variants without the need to raise new antibodies.

Preamplification-free viral RNA diagnostics with single-nucleotide resolution using MARVE, an origami paper-based colorimetric nucleic acid test

The evolution and mutation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are urgent concerns as they pose the risk of vaccine failure and increased viral transmission. However, affordable and scalable tools allowing rapid identification of SARS-CoV-2 variants are not readily available, which impedes diagnosis and epidemiological surveillance. Here we present a colorimetric nucleic acid assay named MARVE (multiplexed, preamplification-free, single-nucleotide-resolved viral evolution) that is convenient to perform and yields single-nucleotide resolution. The assay integrates nucleic acid strand displacement reactions with enzymatic amplification to colorimetrically sense viral RNA using a metal ion-incorporated DNA probe (TEprobe). We provide detailed guidelines to design TEprobes for discriminating single-nucleotide variations in viral RNAs, and to fabricate a test paper for the detection of SARS-CoV-2 variants of concern. Compared with other nucleic acid assays, our assay is preamplification-free, single-nucleotide-resolvable and results are visible via a color change. Besides, it is smartphone readable, multiplexed, quick and cheap ($0.30 per test). The protocol takes ~2 h to complete, from the design and preparation of the DNA probes and test papers (~1 h) to the detection of SARS-CoV-2 or its variants (30-45 min). The design of the TEprobes requires basic knowledge of molecular biology and familiarity with NUPACK or the Python programming language. The fabrication of the origami papers requires access to a wax printer using the CAD and PDF files provided or requires users to be familiar with AutoCAD to design new origami papers. The protocol is also applicable for designing assays to detect other pathogens and their variants.

Inhibitors of SARS-CoV-2 Main Protease (Mpro) as Anti-Coronavirus Agents

The main protease (Mpro) of SARS-CoV-2 is an essential enzyme that plays a critical part in the virus's life cycle, making it a significant target for developing antiviral drugs. The inhibition of SARS-CoV-2 Mpro has emerged as a promising approach for developing therapeutic agents to treat COVID-19. This review explores the structure of the Mpro protein and analyzes the progress made in understanding protein-ligand interactions of Mpro inhibitors. It focuses on binding kinetics, origin, and the chemical structure of these inhibitors. The review provides an in-depth analysis of recent clinical trials involving covalent and non-covalent inhibitors and emerging dual inhibitors targeting SARS-CoV-2 Mpro. By integrating findings from the literature and ongoing clinical trials, this review captures the current state of research into Mpro inhibitors, offering a comprehensive understanding of challenges and directions in their future development as anti-coronavirus agents. This information provides new insights and inspiration for medicinal chemists, paving the way for developing more effective Mpro inhibitors as novel COVID-19 therapies.

A Dynamic and Effective Peptide-Based Strategy for Promptly Addressing Emerging SARS-CoV-2 Variants of Concern

Genomic surveillance based on sequencing the entire genetic code of SARS-CoV-2 involves monitoring and studying genetic changes and variations in disease-causing organisms such as viruses and bacteria. By tracing the virus, it is possible to prevent epidemic spread in the community, ensuring a 'precision public health' strategy. A peptide-based design was applied to provide an efficacious strategy that is able to counteract any emerging viral variant of concern dynamically and promptly to affect the outcomes of a pandemic at an early stage while waiting for the production of the anti-variant-specific vaccine, which require longer times. The inhibition of the interaction between the receptor-binding domain (RBD) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and one of the cellular receptors (DPP4) that its receptors routinely bind to infect human cells is an intriguing therapeutic approach to prevent the virus from entering human cells. Among the other modalities developed for this purpose, peptides surely offer unique advantages, including ease of synthesis, serum stability, low immunogenicity and toxicity, and small production and distribution chain costs. Here, we obtained a potent new inhibitor based on the rearrangement of a previously identified peptide that has been rationally designed on a cell dipeptidyl peptidase 4 (DPP4) sequence, a ubiquitous membrane protein known to bind the RBD-SPIKE domain of the virus. This novel peptide (named DPP4-derived), conceived as an endogenous "drug", is capable of targeting the latest tested variants with a high affinity, reducing the VSV* DG-Fluc pseudovirus Omicron's infection capacity by up to 14%, as revealed by in vitro testing in human Calu-3 cells. Surface plasmon resonance (SPR) confirmed the binding affinity of the new DPP4-derived peptide with Omicron variant RBD.

Barcoded SARS-CoV-2 viruses define the impact of time and route of transmission on the transmission bottleneck in a Syrian hamster model

The transmission bottleneck, defined as the number of viruses that transmit from one host to infect another, is an important determinant of the rate of virus evolution and the level of immunity required to protect against virus transmission. Despite its importance, SARS-CoV-2's transmission bottleneck remains poorly characterized, in part due to a lack of quantitative measurement tools. To address this, we adapted a SARS-CoV-2 reverse genetics system to generate a pool of >200 isogenic SARS-CoV-2 viruses harboring specific 6-nucleotide barcodes inserted in ORF10, a non-translated ORF. We directly inoculated donor Syrian hamsters intranasally with this barcoded virus pool and exposed a paired naïve contact hamster to each donor. Following exposure, the nasal turbinates, trachea, and lungs were collected, viral titers were measured, and the number of barcodes in each tissue were enumerated to quantify the transmission bottleneck. The duration and route (airborne, direct contact, and fomite) of exposure were varied to assess their impact on the transmission bottleneck. In airborne-exposed hamsters, the transmission bottleneck increased with longer exposure durations. We found that direct contact exposure produced the largest transmission bottleneck (average 27 BCs), followed by airborne exposure (average 16 BCs) then fomite exposure (average 8 BCs). Interestingly, we detected unique BCs in both the upper and lower respiratory tract of contact animals from all routes of exposure, suggesting that SARS-CoV-2 can directly infect hamster lungs. Altogether, these findings highlight the utility of barcoded viruses as tools to rigorously study virus transmission. In the future, barcoded SARS-CoV-2 will strengthen studies of immune factors that influence virus transmission.

Comparative Bioinformatic Analysis Reveals Conserved Regions in SARS-CoV-2 Genome for RAPID Pandemic Response

In the face of the SARS-CoV-2 pandemic, characterized by the virus's rapid mutation rates, developing timely and targeted therapeutic and diagnostic interventions presents a significant challenge. This study utilizes bioinformatic analyses to pinpoint conserved genomic regions within SARS-CoV-2, offering a strategic advantage in the fight against this and future pathogens. Our approach has enabled the creation of a diagnostic assay that is not only rapid, reliable, and cost-effective but also possesses a remarkable capacity to detect a wide array of current and prospective variants with unmatched precision. The significance of our findings lies in the demonstration that focusing on these conserved genomic sequences can significantly enhance our preparedness for and response to emerging infectious diseases. By providing a blueprint for the development of versatile diagnostic tools and therapeutics, this research paves the way for a more effective global pandemic response strategy.

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 spike (S) protein. With this approach, we first identified candidate adaptive polymorphisms (CAPs) of the SARS-CoV-2 S protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

Role of SARS-CoV-2 mutations in the evolution of the COVID-19 pandemic

OBJECTIVES

The SARS-CoV-2 pandemic and large-scale genomic surveillance provided an exceptional opportunity to analyze mutations that appeared over three years in viral genomes. Here we studied mutations and their epidemic consequences for SARS-CoV-2 genomes from our center.

METHODS

We analyzed 61,397 SARS-CoV-2 genomes we sequenced from respiratory samples for genomic surveillance. Mutations frequencies were calculated using Nextclade, Microsoft Excel, and an in-house Python script.

RESULTS

A total of 22,225 nucleotide mutations were identified, 220 (1.0%) being each at the root of ≥836 genomes, classifying mutations as 'hyperfertile'. Two seeded the European pandemic: P323L in RNA polymerase, associated with an increased mutation rate, and D614G in spike that improved fitness. Most 'hyperfertile' mutations occurred in areas not predicted with increased virulence. Their mean number was 8±6 (0-22) per 1000 nucleotides per gene. They were 3.7-times more frequent in accessory than informational genes (13.8 versus 3.7/1000 nucleotides). Particularly, they were 4.1-times more frequent in ORF8 than in the RNA polymerase gene. Interestingly, stop codons were present in 97 positions, almost only in accessory genes, including ORF8 (21/100 codons).

CONCLUSIONS

most 'hyperfertile' mutations did not predict emergence of a new epidemic, and some were stop codons indicating the existence of so-named 'non-virulence' genes.

Performance of amplicon and capture based next-generation sequencing approaches for the epidemiological surveillance of Omicron SARS-CoV-2 and other variants of concern

To control the SARS-CoV-2 pandemic, healthcare systems have focused on ramping up their capacity for epidemiological surveillance through viral whole genome sequencing. In this paper, we tested the performance of two protocols of SARS-CoV-2 nucleic acid enrichment, an amplicon enrichment using different versions of the ARTIC primer panel and a hybrid-capture method using KAPA RNA Hypercap. We focused on the challenge of the Omicron variant sequencing, the advantages of automated library preparation and the influence of the bioinformatic analysis in the final consensus sequence. All 94 samples were sequenced using Illumina iSeq 100 and analysed with two bioinformatic pipelines: a custom-made pipeline and an Illumina-owned pipeline. We were unsuccessful in sequencing six samples using the capture enrichment due to low reads. On the other hand, amplicon dropout and mispriming caused the loss of mutation G21987A and the erroneous addition of mutation T15521A respectively using amplicon enrichment. Overall, we found high sequence agreement regardless of method of enrichment, bioinformatic pipeline or the use of automation for library preparation in eight different SARS-CoV-2 variants. Automation and the use of a simple app for bioinformatic analysis can simplify the genotyping process, making it available for more diagnostic facilities and increasing global vigilance.

INSaFLU-TELEVIR: an open web-based bioinformatics suite for viral metagenomic detection and routine genomic surveillance

BACKGROUND

Implementation of clinical metagenomics and pathogen genomic surveillance can be particularly challenging due to the lack of bioinformatics tools and/or expertise. In order to face this challenge, we have previously developed INSaFLU, a free web-based bioinformatics platform for virus next-generation sequencing data analysis. Here, we considerably expanded its genomic surveillance component and developed a new module (TELEVIR) for metagenomic virus identification.

RESULTS

The routine genomic surveillance component was strengthened with new workflows and functionalities, including (i) a reference-based genome assembly pipeline for Oxford Nanopore technologies (ONT) data; (ii) automated SARS-CoV-2 lineage classification; (iii) Nextclade analysis; (iv) Nextstrain phylogeographic and temporal analysis (SARS-CoV-2, human and avian influenza, monkeypox, respiratory syncytial virus (RSV A/B), as well as a "generic" build for other viruses); and (v) algn2pheno for screening mutations of interest. Both INSaFLU pipelines for reference-based consensus generation (Illumina and ONT) were benchmarked against commonly used command line bioinformatics workflows for SARS-CoV-2, and an INSaFLU snakemake version was released. In parallel, a new module (TELEVIR) for virus detection was developed, after extensive benchmarking of state-of-the-art metagenomics software and following up-to-date recommendations and practices in the field. TELEVIR allows running complex workflows, covering several combinations of steps (e.g., with/without viral enrichment or host depletion), classification software (e.g., Kaiju, Kraken2, Centrifuge, FastViromeExplorer), and databases (RefSeq viral genome, Virosaurus, etc.), while culminating in user- and diagnosis-oriented reports. Finally, to potentiate real-time virus detection during ONT runs, we developed findONTime, a tool aimed at reducing costs and the time between sample reception and diagnosis.

CONCLUSIONS

The accessibility, versatility, and functionality of INSaFLU-TELEVIR are expected to supply public and animal health laboratories and researchers with a user-oriented and pan-viral bioinformatics framework that promotes a strengthened and timely viral metagenomic detection and routine genomics surveillance. INSaFLU-TELEVIR is compatible with Illumina, Ion Torrent, and ONT data and is freely available at https://insaflu.insa.pt/ (online tool) and https://github.com/INSaFLU (code).

Co-Mutations and Possible Variation Tendency of the Spike RBD and Membrane Protein in SARS-CoV-2 by Machine Learning

Since the onset of the coronavirus disease 2019 (COVID-19) pandemic, SARS-CoV-2 variants capable of breakthrough infections have attracted global attention. These variants have significant mutations in the receptor-binding domain (RBD) of the spike protein and the membrane (M) protein, which may imply an enhanced ability to evade immune responses. In this study, an examination of co-mutations within the spike RBD and their potential correlation with mutations in the M protein was conducted. The EVmutation method was utilized to analyze the distribution of the mutations to elucidate the relationship between the mutations in the spike RBD and the alterations in the M protein. Additionally, the Sequence-to-Sequence Transformer Model (S2STM) was employed to establish mapping between the amino acid sequences of the spike RBD and M proteins, offering a novel and efficient approach for streamlined sequence analysis and the exploration of their interrelationship. Certain mutations in the spike RBD, G339D-S373P-S375F and Q493R-Q498R-Y505, are associated with a heightened propensity for inducing mutations at specific sites within the M protein, especially sites 3 and 19/63. These results shed light on the concept of mutational synergy between the spike RBD and M proteins, illuminating a potential mechanism that could be driving the evolution of SARS-CoV-2.

Vero cell-adapted SARS-CoV-2 strain shows increased viral growth through furin-mediated efficient spike cleavage

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) utilizes several host proteases to cleave the spike (S) protein to enter host cells. SARS-CoV-2 S protein is cleaved into S1 and S2 subunits by furin, which is closely involved in the pathogenicity of SARS-CoV-2. However, the effects of the modulated protease cleavage activity due to S protein mutations on viral replication and pathogenesis remain unclear. Herein, we serially passaged two SARS-CoV-2 strains in Vero cells and characterized the cell-adapted SARS-CoV-2 strains in vitro and in vivo. The adapted strains showed high viral growth, effective S1/S2 cleavage of the S protein, and low pathogenicity compared with the wild-type strain. Furthermore, the viral growth and S1/S2 cleavage were enhanced by the combination of the Δ68-76 and H655Y mutations using recombinant SARS-CoV-2 strains generated by the circular polymerase extension reaction. The recombinant SARS-CoV-2 strain, which contained the mutation of the adapted strain, showed increased susceptibility to the furin inhibitor, suggesting that the adapted SARS-CoV-2 strain utilized furin more effectively than the wild-type strain. Pathogenicity was attenuated by infection with effectively cleaved recombinant SARS-CoV-2 strains, suggesting that the excessive cleavage of the S proteins decreases virulence. Finally, the high-growth-adapted SARS-CoV-2 strain could be used as the seed for a low-cost inactivated vaccine; immunization with this vaccine can effectively protect the host from SARS-CoV-2 variants. Our findings provide novel insights into the growth and pathogenicity of SARS-CoV-2 in the evolution of cell-cell transmission.

IMPORTANCE

The efficacy of the S protein cleavage generally differs among the SARS-CoV-2 variants, resulting in distinct viral characteristics. The relationship between a mutation and the entry of SARS-CoV-2 into host cells remains unclear. In this study, we analyzed the sequence of high-growth Vero cell-adapted SARS-CoV-2 and factors determining the enhancement of the growth of the adapted virus and confirmed the characteristics of the adapted strain by analyzing the recombinant SARS-CoV-2 strain. We successfully identified mutations Δ68-76 and H655Y, which enhance viral growth and the S protein cleavage by furin. Using recombinant viruses enabled us to conduct a virus challenge experiment in vivo. The pathogenicity of SARS-CoV-2 introduced with the mutations Δ68-76, H655Y, P812L, and Q853L was attenuated in hamsters, indicating the possibility of the attenuation of excessive cleaved SARS-CoV-2. These findings provide novel insights into the infectivity and pathogenesis of SARS-CoV-2 strains, thereby significantly contributing to the field of virology.

Modeling the emergence of viral resistance for SARS-CoV-2 during treatment with an anti-spike monoclonal antibody

To mitigate the loss of lives during the COVID-19 pandemic, emergency use authorization was given to several anti-SARS-CoV-2 monoclonal antibody (mAb) therapies for the treatment of mild-to-moderate COVID-19 in patients with a high risk of progressing to severe disease. Monoclonal antibodies used to treat SARS-CoV-2 target the spike protein of the virus and block its ability to enter and infect target cells. Monoclonal antibody therapy can thus accelerate the decline in viral load and lower hospitalization rates among high-risk patients with variants susceptible to mAb therapy. However, viral resistance has been observed, in some cases leading to a transient viral rebound that can be as large as 3-4 orders of magnitude. As mAbs represent a proven treatment choice for SARS-CoV-2 and other viral infections, evaluation of treatment-emergent mAb resistance can help uncover underlying pathobiology of SARS-CoV-2 infection and may also help in the development of the next generation of mAb therapies. Although resistance can be expected, the large rebounds observed are much more difficult to explain. We hypothesize replenishment of target cells is necessary to generate the high transient viral rebound. Thus, we formulated two models with different mechanisms for target cell replenishment (homeostatic proliferation and return from an innate immune response antiviral state) and fit them to data from persons with SARS-CoV-2 treated with a mAb. We showed that both models can explain the emergence of resistant virus associated with high transient viral rebounds. We found that variations in the target cell supply rate and adaptive immunity parameters have a strong impact on the magnitude or observability of the viral rebound associated with the emergence of resistant virus. Both variations in target cell supply rate and adaptive immunity parameters may explain why only some individuals develop observable transient resistant viral rebound. Our study highlights the conditions that can lead to resistance and subsequent viral rebound in mAb treatments during acute infection.

Targeting G9a translational mechanism of SARS-CoV-2 pathogenesis for multifaceted therapeutics of COVID-19 and its sequalae

By largely unknown mechanism(s), SARS-CoV-2 hijacks the host translation apparatus to promote COVID-19 pathogenesis. We report that the histone methyltransferase G9a noncanonically regulates viral hijacking of the translation machinery to bring about COVID-19 symptoms of hyperinflammation, lymphopenia, and blood coagulation. Chemoproteomic analysis of COVID-19 patient peripheral mononuclear blood cells (PBMC) identified enhanced interactions between SARS-CoV-2-upregulated G9a and distinct translation regulators, particularly the N 6 -methyladenosine (m 6 A) RNA methylase METTL3. These interactions with translation regulators implicated G9a in translational regulation of COVID-19. Inhibition of G9a activity suppressed SARS-CoV-2 replication in human alveolar epithelial cells. Accordingly, multi-omics analysis of the same alveolar cells identified SARS-CoV-2-induced changes at the transcriptional, m 6 A-epitranscriptional, translational, and post-translational (phosphorylation or secretion) levels that were reversed by inhibitor treatment. As suggested by the aforesaid chemoproteomic analysis, these multi-omics-correlated changes revealed a G9a-regulated translational mechanism of COVID-19 pathogenesis in which G9a directs translation of viral and host proteins associated with SARS-CoV-2 replication and with dysregulation of host response. Comparison of proteomic analyses of G9a inhibitor-treated, SARS-CoV-2 infected cells, or ex vivo culture of patient PBMCs, with COVID-19 patient data revealed that G9a inhibition reversed the patient proteomic landscape that correlated with COVID-19 pathology/symptoms. These data also indicated that the G9a-regulated, inhibitor-reversed, translational mechanism outperformed G9a-transcriptional suppression to ultimately determine COVID-19 pathogenesis and to define the inhibitor action, from which biomarkers of serve symptom vulnerability were mechanistically derived. This cell line-to-patient conservation of G9a-translated, COVID-19 proteome suggests that G9a inhibitors can be used to treat patients with COVID-19, particularly patients with long-lasting COVID-19 sequelae.

Incipient functional SARS-CoV-2 diversification identified through neural network haplotype maps

Since its introduction in the human population, SARS-CoV-2 has evolved into multiple clades, but the events in its intrahost diversification are not well understood. Here, we compare three-dimensional (3D) self-organized neural haplotype maps (SOMs) of SARS-CoV-2 from thirty individual nasopharyngeal diagnostic samples obtained within a 19-day interval in Madrid (Spain), at the time of transition between clades 19 and 20. SOMs have been trained with the haplotype repertoire present in the mutant spectra of the nsp12- and spike (S)-coding regions. Each SOM consisted of a dominant neuron (displaying the maximum frequency), surrounded by a low-frequency neuron cloud. The sequence of the master (dominant) neuron was either identical to that of the reference Wuhan-Hu-1 genome or differed from it at one nucleotide position. Six different deviant haplotype sequences were identified among the master neurons. Some of the substitutions in the neural clouds affected critical sites of the nsp12-nsp8-nsp7 polymerase complex and resulted in altered kinetics of RNA synthesis in an in vitro primer extension assay. Thus, the analysis has identified mutations that are relevant to modification of viral RNA synthesis, present in the mutant clouds of SARS-CoV-2 quasispecies. These mutations most likely occurred during intrahost diversification in several COVID-19 patients, during an initial stage of the pandemic, and within a brief time period.

Exploration of low-frequency allelic variants of SARS-CoV-2 genomes reveals coinfections in Mexico occurred during periods of VOCs turnover

A total of 14 973 alleles in 29 661 sequenced samples collected between March 2021 and January 2023 by the Mexican Consortium for Genomic Surveillance (CoViGen-Mex) and collaborators were used to construct a thorough map of mutations of the Mexican SARS-CoV-2 genomic landscape containing Intra-Patient Minor Allelic Variants (IPMAVs), which are low-frequency alleles not ordinarily present in a genomic consensus sequence. This additional information proved critical in identifying putative coinfecting variants included alongside the most common variants, B.1.1.222, B.1.1.519, and variants of concern (VOCs) Alpha, Gamma, Delta, and Omicron. A total of 379 coinfection events were recorded in the dataset (a rate of 1.28 %), resulting in the first such catalogue in Mexico. The most common putative coinfections occurred during the spread of Delta or after the introduction of Omicron BA.2 and its descendants. Coinfections occurred constantly during periods of variant turnover when more than one variant shared the same niche and high infection rate was observed, which was dependent on the local variants and time. Coinfections might occur at a higher frequency than customarily reported, but they are often ignored as only the consensus sequence is reported for lineage identification.

Evolution of a Distinct SARS-CoV-2 Lineage Identified during an Investigation of a Hospital Outbreak

The SARS-CoV-2 virus steadily evolves, and numerous antigenically distinct variants have emerged over the past three years. Tracking the evolution of the virus would help us understand the process that generates the diverse variants and predict the future evolutionary trajectory of SARS-CoV-2. Here, we report the evolutionary trajectory of a unique Omicron lineage identified during an outbreak investigation that occurred in a residence unit in the healthcare system. The new lineage had four distinct non-synonymous and two distinct synonymous mutations apart from its parental lineage. Since this lineage of virus was exclusively found during the outbreak, we were able to track the detailed evolutionary history of the entire lineage along the transmission path. Furthermore, we estimated the evolutionary rate of the SARS-CoV-2 Omicron variant from the analysis of the evolution of the lineage. This new Omicron sub-lineage acquired 3 mutations in a 12-day period, and the evolutionary rate was estimated as 3.05 × 10-3 subs/site/year. This study provides more insight into an ever-evolving virus.

AlphaFold2 Reveals Structural Patterns of Seasonal Haplotype Diversification in SARS-CoV-2 Spike Protein Variants

The slow experimental acquisition of high-quality atomic structures of the rapidly changing proteins of the COVID-19 virus challenges vaccine and therapeutic drug development efforts. Fortunately, deep learning tools such as AlphaFold2 can quickly generate reliable models of atomic structure at experimental resolution. Current modeling studies have focused solely on definitions of mutant constellations of Variants of Concern (VOCs), leaving out the impact of haplotypes on protein structure. Here, we conduct a thorough comparative structural analysis of S-proteins belonging to major VOCs and corresponding latitude-delimited haplotypes that affect viral seasonal behavior. Our approach identified molecular regions of importance as well as patterns of structural recruitment. The S1 subunit hosted the majority of structural changes, especially those involving the N-terminal domain (NTD) and the receptor-binding domain (RBD). In particular, structural changes in the NTD were much greater than just translations in three-dimensional space, altering the sub-structures to greater extents. We also revealed a notable pattern of structural recruitment with the early VOCs Alpha and Delta behaving antagonistically by suppressing regions of structural change introduced by their corresponding haplotypes, and the current VOC Omicron behaving synergistically by amplifying or collecting structural change. Remarkably, haplotypes altering the galectin-like structure of the NTD were major contributors to seasonal behavior, supporting its putative environmental-sensing role. Our results provide an extensive view of the evolutionary landscape of the S-protein across the COVID-19 pandemic. This view will help predict important regions of structural change in future variants and haplotypes for more efficient vaccine and drug development.

Within-host evolution of SARS-CoV-2: how often are de novo mutations transmitted from symptomatic infections?

Despite a relatively low mutation rate, the large number of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections has allowed for substantial genetic change, leading to a multitude of emerging variants. Using a recently determined mutation rate (per site replication), as well as within-host parameter estimates for symptomatic SARS-CoV-2 infection, we apply a stochastic transmission-bottleneck model to describe the survival probability of de novo SARS-CoV-2 mutations as a function of bottleneck size and selection coefficient. For narrow bottlenecks, we find that mutations affecting per-target-cell attachment rate (with phenotypes associated with fusogenicity and ACE2 binding) have similar transmission probabilities to mutations affecting viral load clearance (with phenotypes associated with humoral evasion). We further find that mutations affecting the eclipse rate (with phenotypes associated with reorganization of cellular metabolic processes and synthesis of viral budding precursor material) are highly favoured relative to all other traits examined. We find that mutations leading to reduced removal rates of infected cells (with phenotypes associated with innate immune evasion) have limited transmission advantage relative to mutations leading to humoral evasion. Predicted transmission probabilities, however, for mutations affecting innate immune evasion are more consistent with the range of clinically estimated household transmission probabilities for de novo mutations. This result suggests that although mutations affecting humoral evasion are more easily transmitted when they occur, mutations affecting innate immune evasion may occur more readily. We examine our predictions in the context of a number of previously characterized mutations in circulating strains of SARS-CoV-2. Our work offers both a null model for SARS-CoV-2 mutation rates and predicts which aspects of viral life history are most likely to successfully evolve, despite low mutation rates and repeated transmission bottlenecks.

Temporal Dynamics and Genomic Landscape of SARS-CoV-2 After Four Years of Evolution

Introduction Since its emergence, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has undergone extensive genomic evolution, impacting public health policies, diagnosis, medication, and vaccine development. This study leverages advanced bioinformatics to assess the virus's temporal and regional genomic evolution from December 2019 to October 2023. Methods Our analysis incorporates 16,575 complete SARS-CoV-2 sequences collected from 214 countries. These samples were comparatively analyzed, with a detailed characterization of nucleic mutations, lineages, distribution, and evolutionary patterns during each year, using the Wuhan-Hu-1 strain as the reference. Results Our analysis has identified a total of 21,580 mutations that we classified into transient mutations, which diminished over time, and persistent mutations with steadily increasing frequencies. This mutation landscape led to a notable surge in the evolutionary rate, rising from 13 mutations per sample in 2020 to 96 by 2023, with minor geographic variations. The phylogenetic analysis unveiled three distinct evolutionary branches, each representing unique viral evolution pathways. These lineages exhibited a tendency for a reduced duration of dominance with a shortening prevalence period over time, as dominant strains were consistently replaced by more fit variants. Notably, the emergence of the Alpha and Delta variants in 2021 was followed by the subsequent dominance of Omicron clade variants that have branched into several recombinant variants in 2022, marking a significant shift in the viral landscape. Conclusion This study sheds light on the dynamic nature of SARS-CoV-2 evolution, emphasizing the importance of continuous and vigilant genomic surveillance. The dominance of recombinant lineages, coupled with the disappearance of local variants, underscores the virus's adaptability.

SARS-CoV-2 and innate immunity: the good, the bad, and the "goldilocks"

An ancient conflict between hosts and pathogens has driven the innate and adaptive arms of immunity. Knowledge about this interplay can not only help us identify biological mechanisms but also reveal pathogen vulnerabilities that can be leveraged therapeutically. The humoral response to SARS-CoV-2 infection has been the focus of intense research, and the role of the innate immune system has received significantly less attention. Here, we review current knowledge of the innate immune response to SARS-CoV-2 infection and the various means SARS-CoV-2 employs to evade innate defense systems. We also consider the role of innate immunity in SARS-CoV-2 vaccines and in the phenomenon of long COVID.

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 Spike (S) protein. With this approach, we first identified Candidate Adaptive Polymorphisms (CAPs) of the SARS-CoV-2 Spike protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

Differential laboratory passaging of SARS-CoV-2 viral stocks impacts the in vitro assessment of neutralizing antibodies

Viral populations in natural infections can have a high degree of sequence diversity, which can directly impact immune escape. However, antibody potency is often tested in vitro with a relatively clonal viral populations, such as laboratory virus or pseudotyped virus stocks, which may not accurately represent the genetic diversity of circulating viral genotypes. This can affect the validity of viral phenotype assays, such as antibody neutralization assays. To address this issue, we tested whether recombinant virus carrying SARS-CoV-2 spike (VSV-SARS-CoV-2-S) stocks could be made more genetically diverse by passage, and if a stock passaged under selective pressure was more capable of escaping monoclonal antibody (mAb) neutralization than unpassaged stock or than viral stock passaged without selective pressures. We passaged VSV-SARS-CoV-2-S four times concurrently in three cell lines and then six times with or without polyclonal antiserum selection pressure. All three of the monoclonal antibodies tested neutralized the viral population present in the unpassaged stock. The viral inoculum derived from serial passage without antiserum selection pressure was neutralized by two of the three mAbs. However, the viral inoculum derived from serial passage under antiserum selection pressure escaped neutralization by all three mAbs. Deep sequencing revealed the rapid acquisition of multiple mutations associated with antibody escape in the VSV-SARS-CoV-2-S that had been passaged in the presence of antiserum, including key mutations present in currently circulating Omicron subvariants. These data indicate that viral stock that was generated under polyclonal antiserum selection pressure better reflects the natural environment of the circulating virus and may yield more biologically relevant outcomes in phenotypic assays. Thus, mAb assessment assays that utilize a more genetically diverse, biologically relevant, virus stock may yield data that are relevant for prediction of mAb efficacy and for enhancing biosurveillance.

Ultrasensitive and fast detection of SARS-CoV-2 using RT-LAMP without pH-dependent dye

This study investigates the performance of reverse transcription loop-mediated isothermal amplification (RT-LAMP) assay for the colorimetric detection of SARS-CoV-2 using fluorometric dye, namely, calcein. The detection limit (LoD) with the N-ID1 primer set resulted in superior performance, corresponding to ~ 2 copies/reaction or ~ 0.1 copies/μL of the RNA sample. The color development can be observed by the naked eye, using an ultraviolet (UV) transilluminator or a hand-UV light without the requirement of expensive devices. The average time-to-reaction (TTR) value was 26.2 min in high-copy number samples, while it was about 50 min in rRT-PCR. A mobile application was proposed to quantify the positive and negative results based on the three-color spaces (RGB, Lab, and HSB). Compared to rRT-PCR (n = 67), this assay allows fast and sensitive visual detection of SARS-CoV-2, with high sensitivity (90.9%), selectivity (100%), and accuracy (94.03%). Besides, the assay was sensitive regardless of variants. Since this assay uses a fluorescent dye for visual observation, it can be easily adapted in RT-LAMP assays with high sensitivity. Thus, it can be utilized in low-source centers and field testing such as conferences, sports meetings, refugee camps, companies, and schools.

SARS-CoV-2 variant biology and immune evasion

This chapter discusses the SARS-CoV-2 variants and their immune evasion strategies, shedding light on the dynamic nature of the COVID-19 pandemic. The ecological dynamics and viral evolution of SARS-CoV-2 are explored, considering carriers of infection, individual immunity profiles, and human movement as key factors in the emergence and dissemination of variants. The chapter discusses SARS-CoV-2 mutation, including mutation rate, substitution rate, and recombination, influencing genetic diversity and evolution. Transmission bottlenecks are highlighted as determinants of dominant variants during viral spread. The evolution phases of the pandemic are outlined, from limited early evolution to the emergence of notable changes like the D614G substitution and variants with heavy mutations. Variants of Concern (VOCs), including Alpha, Beta, Gamma, and the recent Omicron variant, are examined, with insights into inter-lineage and intra-lineage dynamics. The origin of VOCs and the Omicron variant is explored, alongside the role of the furin cleavage site (FCS) in variant emergence. The impact of structural and non-structural proteins on viral infectivity is assessed, as well as innate immunity evasion strategies employed by SARS-CoV-2 variants. The chapter concludes by considering future possibilities, including ongoing virus evolution, the need for surveillance, vaccine development, and public health measures.

Transmission Bottleneck Size Estimation from De Novo Viral Genetic Variation

Sequencing of viral infections has become increasingly common over the last decade. Deep sequencing data in particular have proven useful in characterizing the roles that genetic drift and natural selection play in shaping within-host viral populations. They have also been used to estimate transmission bottleneck sizes from identified donor-recipient pairs. These bottleneck sizes quantify the number of viral particles that establish genetic lineages in the recipient host and are important to estimate due to their impact on viral evolution. Current approaches for estimating bottleneck sizes exclusively consider the subset of viral sites that are observed as polymorphic in the donor individual. However, these approaches have the potential to substantially underestimate true transmission bottleneck sizes. Here, we present a new statistical approach for instead estimating bottleneck sizes using patterns of viral genetic variation that arise de novo within a recipient individual. Specifically, our approach makes use of the number of clonal viral variants observed in a transmission pair, defined as the number of viral sites that are monomorphic in both the donor and the recipient but carry different alleles. We first test our approach on a simulated dataset and then apply it to both influenza A virus sequence data and SARS-CoV-2 sequence data from identified transmission pairs. Our results confirm the existence of extremely tight transmission bottlenecks for these 2 respiratory viruses.

Domestic and wild animal samples and diagnostic testing for SARS-CoV-2

From the first cases in 2019, COVID-19 infections caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have resulted in over 6 million human deaths in a worldwide pandemic. SARS-CoV-2 is commonly spread from human to human through close contact and is capable of infecting both humans and animals. Worldwide, there have been over 675 animal outbreaks reported that resulted in over 2000 animal infections including domestic and wild animals. As the role of animal infections in the transmission, pathogenesis, and evolution of SARS-CoV-2 is still unfolding, accurate and reliable animal diagnostic tests are critical to aid in managing both human and animal health. This review highlights key animal samples and the three main diagnostic approaches used for animal testing: PCR, serology, and Next Generation Sequencing. Diagnostic results help inform (often difficult) clinical decision-making, but also possible ways to mitigate spread among pets, food supplies, or wildlife. A One Health approach has been key to monitoring the SARS-CoV-2 pandemic, as consistent human-animal interactions can lead to novel variants. Having multiple animal diagnostic tests for SARS-CoV-2 available is critical to ensure human, animal, and environmental health.

On the space of SARS-CoV-2 genetic sequence variants

The coronavirus pandemic caused by the SARS-CoV-2 virus, which humanity resisted using the latest advances in science, left behind, among other things, extensive genetic data. Every day since the end of 2019, samples of the virus genomes have been collected around the world, which makes it possible to trace its evolution in detail from its emergence to the present. The accumulated statistics of testing results showed that the number of confirmed cases of SARS-CoV-2 infection was at least 767.5 million (9.5 % of the current world population, excluding asymptomatic people), and the number of sequenced virus genomes is more than 15.7 million (which is over 2 % of the total number of infected people). These new data potentially contain information about the mechanisms of the variability and spread of the virus, its interaction with the human immune system, the main parameters characterizing the mechanisms of the development of a pandemic, and much more. In this article, we analyze the space of possible variants of SARS-CoV-2 genetic sequences both from a mathematical point of view and taking into account the biological limitations inherent in this system, known both from general biological knowledge and from the consideration of the characteristics of this particular virus. We have developed software capable of loading and analyzing SARS-CoV-2 nucleotide sequences in FASTA format, determining the 5' and 3' UTR positions, the number and location of unidentified nucleotides ("N"), performing alignment with the reference sequence by calling the program designed for this, determining mutations, deletions and insertions, as well as calculating various characteristics of virus genomes with a given time step (days, weeks, months, etc.). The data obtained indicate that, despite the apparent mathematical diversity of possible options for changing the virus over time, the corridor of the evolutionary trajectory that the coronavirus has passed through seems to be quite narrow. Thus it can be assumed that it is determined to some extent, which allows us to hope for a possibility of modeling the evolution of the coronavirus.

Boosted production of antibodies that neutralized different SARS-CoV-2 variants in a COVID-19 convalescent following messenger RNA vaccination - a case study

Vaccinated convalescents do not develop severe COVID-19 after infection with new SARS-CoV-2 variants. We questioned how messenger RNA (mRNA) vaccination of convalescents provides protection from emerging virus variants. From the cohort of 71 convalescent plasma donors, we identified a patient who developed immune response to infection with SARS-CoV-2 variant of 20A clade and who subsequently received mRNA vaccine encoding spike (S) protein of strain of 19A clade. We showed that vaccination increased the production of immune cells and anti-S antibodies in the serum. Serum antibodies neutralized not only 19A and 20A, but also 20B, 20H, 21J, and 21K virus variants. One of the serum antibodies (100F8) completely neutralized 20A, 21J, and partially 21K strains. 100F8 was structurally similar to published Ab188 antibody, which recognized non-conserved epitope on the S protein. We proposed that 100F8 and other serum antibodies of the patient which recognized non- and conserved epitopes of the S protein, could have additive or synergistic effects to neutralize various virus variants. Thus, mRNA vaccination could be beneficial for convalescents because it boosts production of neutralizing antibodies with broad-spectrum activity.

SARS-CoV-2 Orphan Gene ORF10 Contributes to More Severe COVID-19 Disease

The orphan gene of SARS-CoV-2, ORF10, is the least studied gene in the virus responsible for the COVID-19 pandemic. Recent experimentation indicated ORF10 expression moderates innate immunity in vitro. However, whether ORF10 affects COVID-19 in humans remained unknown. We determine that the ORF10 sequence is identical to the Wuhan-Hu-1 ancestral haplotype in 95% of genomes across five variants of concern (VOC). Four ORF10 variants are associated with less virulent clinical outcomes in the human host: three of these affect ORF10 protein structure, one affects ORF10 RNA structural dynamics. RNA-Seq data from 2070 samples from diverse human cells and tissues reveals ORF10 accumulation is conditionally discordant from that of other SARS-CoV-2 transcripts. Expression of ORF10 in A549 and HEK293 cells perturbs immune-related gene expression networks, alters expression of the majority of mitochondrially-encoded genes of oxidative respiration, and leads to large shifts in levels of 14 newly-identified transcripts. We conclude ORF10 contributes to more severe COVID-19 clinical outcomes in the human host.

Low Risk of SARS-CoV-2 Reinfection for Fully or Boosted mRNA Vaccinated Subjects in Sicily: A Population-Based Study Using Real-World Data

Background: Reinfections occur as a response to natural infections wanes and novel strains of SARS-CoV-2 emerge. The present research explored the correlation between sex, age, COVID-19 vaccination, prior infection hospitalization, and SARS-CoV-2 reinfection in Sicily, Italy. Materials and Methods: A population-based retrospective cohort study was articulated using the vaccination flux from a regional registry and the Sicilian COVID-19 monitoring system of the Italian Institute of Health. Only adult Sicilians were included in the study, and hazard ratios were calculated using Cox regression. Results: Partial vaccination provided some protection (adj-HR: 0.92), when compared to unvaccinated individuals; furthermore, reinfection risk was reduced by full vaccination (adj-HR: 0.43), and the booster dose (adj-HR: 0.41). Males had a lower risk than females of reinfection with SARS-CoV-2 (adj-HR: 0.75). Reinfection with SARS-CoV-2 was diminished by hospitalization during the first infection (adj-HR: 0.78). Reinfection risk was higher among those aged 30-39 and 40-49 compared to those aged 18-29, whereas those aged 60-69, 70-79, and 80+ were statistically protected. Reinfection was significantly more frequent during the wild-type-Alpha, Delta, Delta-Omicron, and Omicron dominance/codominance waves compared to the wild type. Conclusions: This study establishes a solid base for comprehending the reinfection phenomenon in Sicily by pinpointing the most urgent policy hurdles and identifying some of the major factors. COVID-19 vaccination, one of the most effective public health tools, protects against reinfection, mostly caused by the Omicron strain. Elderly and hospitalized people's lower risk suggests stricter PPE use.

The Functional Implications of Broad Spectrum Bioactive Compounds Targeting RNA-Dependent RNA Polymerase (RdRp) in the Context of the COVID-19 Pandemic

BACKGROUND

As long as COVID-19 endures, viral surface proteins will keep changing and new viral strains will emerge, rendering prior vaccines and treatments decreasingly effective. To provide durable targets for preventive and therapeutic agents, there is increasing interest in slowly mutating viral proteins, including non-surface proteins like RdRp.

METHODS

A scoping review of studies was conducted describing RdRp in the context of COVID-19 through MEDLINE/PubMed and EMBASE. An iterative approach was used with input from content experts and three independent reviewers, focused on studies related to either RdRp activity inhibition or RdRp mechanisms against SARS-CoV-2.

RESULTS

Of the 205 records screened, 43 studies were included in the review. Twenty-five evaluated RdRp activity inhibition, and eighteen described RdRp mechanisms of existing drugs or compounds against SARS-CoV-2. In silico experiments suggested that RdRp inhibitors developed for other RNA viruses may be effective in disrupting SARS-CoV-2 replication, indicating a possible reduction of disease progression from current and future variants. In vitro, in vivo, and human clinical trial studies were largely consistent with these findings.

CONCLUSIONS

Future risk mitigation and treatment strategies against forthcoming SARS-CoV-2 variants should consider targeting RdRp proteins instead of surface proteins.

Structural understanding of SARS-CoV-2 virus entry to host cells

Coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a major global health concern associated with millions of fatalities worldwide. Mutant variants of the virus have further exacerbated COVID-19 mortality and infection rates, emphasizing the urgent need for effective preventive strategies. Understanding the viral infection mechanism is crucial for developing therapeutics and vaccines. The entry of SARS-CoV-2 into host cells is a key step in the infection pathway and has been targeted for drug development. Despite numerous reviews of COVID-19 and the virus, there is a lack of comprehensive reviews focusing on the structural aspects of viral entry. In this review, we analyze structural changes in Spike proteins during the entry process, dividing the entry process into prebinding, receptor binding, proteolytic cleavage, and membrane fusion steps. By understanding the atomic-scale details of viral entry, we can better target the entry step for intervention strategies. We also examine the impacts of mutations in Spike proteins, including the Omicron variant, on viral entry. Structural information provides insights into the effects of mutations and can guide the development of therapeutics and vaccines. Finally, we discuss available structure-based approaches for the development of therapeutics and vaccines. Overall, this review provides a detailed analysis of the structural aspects of SARS-CoV-2 viral entry, highlighting its significance in the development of therapeutics and vaccines against COVID-19. Therefore, our review emphasizes the importance of structural information in combating SARS-CoV-2 infection.

A Bayesian inference method to estimate transmission trees with multiple introductions; applied to SARS-CoV-2 in Dutch mink farms

Knowledge of who infected whom during an outbreak of an infectious disease is important to determine risk factors for transmission and to design effective control measures. Both whole-genome sequencing of pathogens and epidemiological data provide useful information about the transmission events and underlying processes. Existing models to infer transmission trees usually assume that the pathogen is introduced only once from outside into the population of interest. However, this is not always true. For instance, SARS-CoV-2 is suggested to be introduced multiple times in mink farms in the Netherlands from the SARS-CoV-2 pandemic among humans. Here, we developed a Bayesian inference method combining whole-genome sequencing data and epidemiological data, allowing for multiple introductions of the pathogen in the population. Our method does not a priori split the outbreak into multiple phylogenetic clusters, nor does it break the dependency between the processes of mutation, within-host dynamics, transmission, and observation. We implemented our method as an additional feature in the R-package phybreak. On simulated data, our method correctly identifies the number of introductions, with an accuracy depending on the proportion of all observed cases that are introductions. Moreover, when a single introduction was simulated, our method produced similar estimates of parameters and transmission trees as the existing package. When applied to data from a SARS-CoV-2 outbreak in Dutch mink farms, the method provides strong evidence for independent introductions of the pathogen at 13 farms, infecting a total of 63 farms. Using the new feature of the phybreak package, transmission routes of a more complex class of infectious disease outbreaks can be inferred which will aid infection control in future outbreaks.

Inferring Viral Transmission Pathways from Within-Host Variation

Genome sequencing can offer critical insight into pathogen spread in viral outbreaks, but existing transmission inference methods use simplistic evolutionary models and only incorporate a portion of available genetic data. Here, we develop a robust evolutionary model for transmission reconstruction that tracks the genetic composition of within-host viral populations over time and the lineages transmitted between hosts. We confirm that our model reliably describes within-host variant frequencies in a dataset of 134,682 SARS-CoV-2 deep-sequenced genomes from Massachusetts, USA. We then demonstrate that our reconstruction approach infers transmissions more accurately than two leading methods on synthetic data, as well as in a controlled outbreak of bovine respiratory syncytial virus and an epidemiologically-investigated SARS-CoV-2 outbreak in South Africa. Finally, we apply our transmission reconstruction tool to 5,692 outbreaks among the 134,682 Massachusetts genomes. Our methods and results demonstrate the utility of within-host variation for transmission inference of SARS-CoV-2 and other pathogens, and provide an adaptable mathematical framework for tracking within-host evolution.

Adapting Neutralizing Antibodies to Viral Variants by Structure-Guided Affinity Maturation Using Phage Display Technology

Neutralizing monoclonal antibodies have achieved great efficacy and safety for the treatment of numerous infectious diseases. However, their neutralization potency is often rapidly lost when the target antigen mutates. Instead of isolating new antibodies each time a pathogen variant arises, it can be attractive to adapt existing antibodies, making them active against the new variant. Potential benefits of this approach include reduced development time, cost, and regulatory burden. Here a methodology is described to rapidly evolve neutralizing antibodies of proven activity, improving their function against new pathogen variants without losing efficacy against previous ones. The reported procedure is based on structure-guided affinity maturation using combinatorial mutagenesis and phage display technology. Its use against the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is demonstrated, but it is suitable for any other pathogen. As proof of concept, the method is applied to CoV-X2, a human bispecific antibody that binds with high affinity to the early SARS-CoV-2 variants but lost neutralization potency against Delta. Antibodies emerging from the affinity maturation selection exhibit significantly improved neutralization potency against Delta and no loss of efficacy against the other viral sequences tested. These results illustrate the potential application of structure-guided affinity maturation in facilitating the rapid adaptation of neutralizing antibodies to pathogen variants.

Humoral Responses Elicited by SARS-CoV-2 mRNA Vaccine in People Living with HIV

While mRNA SARS-CoV-2 vaccination elicits strong humoral responses in the general population, humoral responses in people living with HIV (PLWH) remain to be clarified. Here, we conducted a longitudinal study of vaccine immunogenicity elicited after two and three doses of mRNA SARS-CoV-2 vaccine in PLWH stratified by their CD4 count. We measured the capacity of the antibodies elicited by vaccination to bind the Spike glycoprotein of different variants of concern (VOCs). We also evaluated the Fc-mediated effector functions of these antibodies by measuring their ability to eliminate CEM.NKr cells stably expressing SARS-CoV-2 Spikes. Finally, we measured the relative capacity of the antibodies to neutralize authentic SARS-CoV-2 virus after the third dose of mRNA vaccine. We found that after two doses of SARS-CoV-2 mRNA vaccine, PLWH with a CD4 count < 250/mm3 had lower levels of anti-RBD IgG antibodies compared to PLWH with a CD4 count > 250/mm3 (p < 0.05). A third dose increased these levels and importantly, no major differences were observed in their capacity to mediate Fc-effector functions and neutralize authentic SARS-CoV-2. Overall, our work demonstrates the importance of mRNA vaccine boosting in immuno-compromised individuals presenting low levels of CD4.

Within-host diversity improves phylogenetic and transmission reconstruction of SARS-CoV-2 outbreaks

Accurate inference of who infected whom in an infectious disease outbreak is critical for the delivery of effective infection prevention and control. The increased resolution of pathogen whole-genome sequencing has significantly improved our ability to infer transmission events. Despite this, transmission inference often remains limited by the lack of genomic variation between the source case and infected contacts. Although within-host genetic diversity is common among a wide variety of pathogens, conventional whole-genome sequencing phylogenetic approaches exclusively use consensus sequences, which consider only the most prevalent nucleotide at each position and therefore fail to capture low-frequency variation within samples. We hypothesized that including within-sample variation in a phylogenetic model would help to identify who infected whom in instances in which this was previously impossible. Using whole-genome sequences from SARS-CoV-2 multi-institutional outbreaks as an example, we show how within-sample diversity is partially maintained among repeated serial samples from the same host, it can transmitted between those cases with known epidemiological links, and how this improves phylogenetic inference and our understanding of who infected whom. Our technique is applicable to other infectious diseases and has immediate clinical utility in infection prevention and control.

Key Considerations during the Transition from the Acute Phase of the COVID-19 Pandemic: A Narrative Review

The COVID-19 pandemic has been met with an unprecedented response from the scientific community, leading to the development, investigation, and authorization of vaccines and antivirals, ultimately reducing the impact of SARS-CoV-2 on global public health. However, SARS-CoV-2 is far from being eradicated, continues to evolve, and causes substantial health and economic burdens. In this narrative review, we posit essential points on SARS-CoV-2 and its responsible management during the transition from the acute phase of the COVID-19 pandemic. As discussed, despite Omicron (sub)variant(s) causing clinically milder infections, SARS-CoV-2 is far from being a negligible pathogen. It requires continued genomic surveillance, particularly if one considers that its future (sub)lineages do not necessarily have to be milder. Antivirals and vaccines remain the essential elements in COVID-19 management. However, the former could benefit from further development and improvements in dosing, while the seasonal administration of the latter requires simplification to increase interest and tackle vaccine hesitancy. It is also essential to ensure the accessibility of COVID-19 pharmaceuticals and vaccines in low-income countries and improve the understanding of their use in the context of the long-term goals of SARS-CoV-2 management. Regardless of location, the primary role of COVID-19 awareness and education must be played by healthcare workers, who directly communicate with patients and serve as role models for healthy behaviors.

Modeling the emergence of viral resistance for SARS-CoV-2 during treatment with an anti-spike monoclonal antibody

The COVID-19 pandemic has led to over 760 million cases and 6.9 million deaths worldwide. To mitigate the loss of lives, emergency use authorization was given to several anti-SARS-CoV-2 monoclonal antibody (mAb) therapies for the treatment of mild-to-moderate COVID-19 in patients with a high risk of progressing to severe disease. Monoclonal antibodies used to treat SARS-CoV-2 target the spike protein of the virus and block its ability to enter and infect target cells. Monoclonal antibody therapy can thus accelerate the decline in viral load and lower hospitalization rates among high-risk patients with susceptible variants. However, viral resistance has been observed, in some cases leading to a transient viral rebound that can be as large as 3-4 orders of magnitude. As mAbs represent a proven treatment choice for SARS-CoV-2 and other viral infections, evaluation of treatment-emergent mAb resistance can help uncover underlying pathobiology of SARS-CoV-2 infection and may also help in the development of the next generation of mAb therapies. Although resistance can be expected, the large rebounds observed are much more difficult to explain. We hypothesize replenishment of target cells is necessary to generate the high transient viral rebound. Thus, we formulated two models with different mechanisms for target cell replenishment (homeostatic proliferation and return from an innate immune response anti-viral state) and fit them to data from persons with SARS-CoV-2 treated with a mAb. We showed that both models can explain the emergence of resistant virus associated with high transient viral rebounds. We found that variations in the target cell supply rate and adaptive immunity parameters have a strong impact on the magnitude or observability of the viral rebound associated with the emergence of resistant virus. Both variations in target cell supply rate and adaptive immunity parameters may explain why only some individuals develop observable transient resistant viral rebound. Our study highlights the conditions that can lead to resistance and subsequent viral rebound in mAb treatments during acute infection.

Overview of immunological & virological factors driving the evolution & global spread of SARS-CoV-2 variants

The SARS-CoV-2, a highly infectious positive strand RNA virus first identified in December 2019, has produced multiple genetic variants that have rapidly and sequentially spread worldwide during the coronavirus disease 2019 (COVID-19) pandemic. Genetic changes in SARS-CoV-2 for greater infectivity, replication and transmission were selected during the early stages of the pandemic. More recently, after widespread infection and vaccination, SARS-CoV-2 variants that evade antigen-specific adaptive immunity, have begun to be selected. This article provides an overview of the molecular immunological and virological factors underlying the origin and global spread of important SARS-CoV-2 variant lineages.

Transmission bottleneck size estimation from de novo viral genetic variation

Sequencing of viral infections has become increasingly common over the last decade. Deep sequencing data in particular have proven useful in characterizing the roles that genetic drift and natural selection play in shaping within-host viral populations. They have also been used to estimate transmission bottleneck sizes from identified donor-recipient pairs. These bottleneck sizes quantify the number of viral particles that establish genetic lineages in the recipient host and are important to estimate due to their impact on viral evolution. Current approaches for estimating bottleneck sizes exclusively consider the subset of viral sites that are observed as polymorphic in the donor individual. However, allele frequencies can change dramatically over the course of an individual's infection, such that sites that are polymorphic in the donor at the time of transmission may not be polymorphic in the donor at the time of sampling and allele frequencies at donor-polymorphic sites may change dramatically over the course of a recipient's infection. Because of this, transmission bottleneck sizes estimated using allele frequencies observed at a donor's polymorphic sites may be considerable underestimates of true bottleneck sizes. Here, we present a new statistical approach for instead estimating bottleneck sizes using patterns of viral genetic variation that arose de novo within a recipient individual. Specifically, our approach makes use of the number of clonal viral variants observed in a transmission pair, defined as the number of viral sites that are monomorphic in both the donor and the recipient but carry different alleles. We first test our approach on a simulated dataset and then apply it to both influenza A virus sequence data and SARS-CoV-2 sequence data from identified transmission pairs. Our results confirm the existence of extremely tight transmission bottlenecks for these two respiratory viruses, using an approach that does not tend to underestimate transmission bottleneck sizes.

Metabolic Associated Fatty Liver Disease (MAFLD) and COVID-19 Infection: An Independent Predictor of Poor Disease Outcome?

Background and Objectives: Early reports on COVID-19 infection suggested that the SARS-CoV-2 virus solely attacks respiratory tract cells. As the pandemic spread, it became clear that the infection is multiorganic. Metabolic associated fatty liver disease (MAFLD) is a chronic liver disease strongly associated with insulin resistance and diabetes. The aim of this study was to assess a possible interplay between MAFLD and COVID-19 infection and its implication in COVID-19 outcome. Materials and Methods: A retrospective observational study, including 130 COVID-19 positive patients was conducted. MAFLD diagnosis was made based on the International Consensus criteria. Patients were divided into two groups, group A (MAFLD) and group B (nonMAFLD). Anthropometric and laboratory analysis were obtained. COVID-19 severity was assessed using the NEWS2 score. Disease outcome was threefold and regarded as discharged, patients who required mechanical ventilation (MV), and deceased patients. Results: MAFLD prevalence was 42%, 67% of patients were discharged, and 19% needed MV. Mortality rate was 14%. MAFLD patients were significantly younger (p < 0.001), and had higher body mass index (p < 0.05), respiratory rate (p < 0.05) and systolic blood pressure (p < 0.05) than nonMAFLD patients. Regarding metabolic syndrome and inflammatory markers: group A had significantly higher glycemia at admission (p = 0.008), lower HDL-c (p < 0.01), higher triglycerides (p < 0.01), CRP (p < 0.001), IL-6 (p < 0.05) and ferritin (p < 0.05) than group B. MAFLD was associated with more prevalent type 2 diabetes (p = 0.035) and hypertension (p < 0.05). MAFLD patients had a more severe disease course (NEWS2 score, 6.5 ± 0.5 vs. 3 ± 1.0, p < 0.05). MAFLD presence was associated with lower patient discharge (p < 0.01) and increased need for MV (p = 0.024). Multiple regression analysis showed that BMI (p = 0.045), IL-6 (p = 0.03), and MAFLD (p < 0.05) are significant independent risk factors for a poor COVID-19 outcome. Conclusions: The prevalence of MAFLD is relatively high. MAFLD patients had a more severe COVID-19 clinical course and worse disease outcome. Our results imply that early patient stratification and risk assessment are mandatory in order to avoid poor outcomes.

A critical analysis of the current state of virus taxonomy

Taxonomical classification has preceded evolutionary understanding. For that reason, taxonomy has become a battleground fueled by knowledge gaps, technical limitations, and a priorism. Here we assess the current state of the challenging field, focusing on fallacies that are common in viral classification. We emphasize that viruses are crucial contributors to the genomic and functional makeup of holobionts, organismal communities that behave as units of biological organization. Consequently, viruses cannot be considered taxonomic units because they challenge crucial concepts of organismality and individuality. Instead, they should be considered processes that integrate virions and their hosts into life cycles. Viruses harbor phylogenetic signatures of genetic transfer that compromise monophyly and the validity of deep taxonomic ranks. A focus on building phylogenetic networks using alignment-free methodologies and molecular structure can help mitigate the impasse, at least in part. Finally, structural phylogenomic analysis challenges the polyphyletic scenario of multiple viral origins adopted by virus taxonomy, defeating a polyphyletic origin and supporting instead an ancient cellular origin of viruses. We therefore, prompt abandoning deep ranks and urgently reevaluating the validity of taxonomic units and principles of virus classification.

Evolutionary safety of lethal mutagenesis driven by antiviral treatment

Nucleoside analogs are a major class of antiviral drugs. Some act by increasing the viral mutation rate causing lethal mutagenesis of the virus. Their mutagenic capacity, however, may lead to an evolutionary safety concern. We define evolutionary safety as a probabilistic assurance that the treatment will not generate an increased number of mutants. We develop a mathematical framework to estimate the total mutant load produced with and without mutagenic treatment. We predict rates of appearance of such virus mutants as a function of the timing of treatment and the immune competence of patients, employing realistic assumptions about the vulnerability of the viral genome and its potential to generate viable mutants. We focus on the case study of Molnupiravir, which is an FDA-approved treatment against Coronavirus Disease-2019 (COVID-19). We estimate that Molnupiravir is narrowly evolutionarily safe, subject to the current estimate of parameters. Evolutionary safety can be improved by restricting treatment with this drug to individuals with a low immunological clearance rate and, in future, by designing treatments that lead to a greater increase in mutation rate. We report a simple mathematical rule to determine the fold increase in mutation rate required to obtain evolutionary safety that is also applicable to other pathogen-treatment combinations.