Host-directed editing of the SARS-CoV-2 genome

Evolution of drug resistance against antiviral agents that target cellular factors

Antiviral drugs have classically been developed by directly disrupting the functions of viral proteins. However, this strategy has been largely unsuccessful due to the rapid generation of viral escape mutants. It has been well established that as compared to the virus-centric approach, the strategy of developing antiviral drugs by targeting host-dependency factors (HDFs) minimizes drug resistance. However, recent reports have indicated that drug resistance against some of the host-targeting antiviral agents can in fact occur under some circumstances. Long-term selection pressure of a host-targeting antiviral agent may induce the virus to use an alternate cellular factor or alters its affinity towards the target that confers resistance. Alternatively, virus may synchronize its life cycle with the patterns of drug therapy. In addition, virus may subvert host's immune system to perpetuate under the limiting conditions of the targeted cellular factor. This review describes novel potential mechanisms that may account for the acquiring resistance against agents that target HDFs.

Impact of vaccination on SARS-CoV-2 evolution and immune escape variants

Vaccines and host genetic factors can influence the SARS-CoV-2 evolution and emergence of new variants. Even vaccinated cases get affected as virus continues to evolve, raising concerns about vaccine efficacy and the emergence of immune escape variants. Here, we have analyzed 2295 whole-genome sequences of SARS-CoV-2 collected from vaccinated and unvaccinated cases to evaluate the impact of vaccines on virus diversity within hosts. Our comparative analysis revealed a significant higher incidence of intra-host single nucleotides variants (iSNVs) in vaccinated cases compared to unvaccinated ones (p value<0.0001). Furthermore, we have found that specific mutational processes, including APOBEC (C > T) mediated and ADAR1 (A > G) mediated mutations, were found more prevalent in vaccinated cases. Vaccinated cases exhibited higher accumulation of nonsynonymous mutation than unvaccinated cases. Fixed iSNVs were predominantly located in the ORF1ab and spike genes, several key omicron defining immune escape variants S477N, Q493R, Q498R, Y505H, L452R, and N501Y were identified in the RBD domain of spike gene in vaccinated cases. Our findings suggest that vaccine plays an important role in the evolution of the virus genome. The virus genome acquires random mutations due to error-prone replication of the virus, host modification through APOBEC and ADAR1 mediated editing mechanism, and oxidative stress. These mutations become fixed in the viral population due to the selective pressure imposed by vaccination.

Synergism between remdesivir and ribavirin leads to SARS-CoV-2 extinction in cell culture

BACKGROUND AND PURPOSE

There is a need for effective anti-COVID-19 treatments, mainly for individuals at risk of severe disease such as the elderly and the immunosuppressed. Drug repositioning has proved effective in identifying drugs that can find a new application for the control of coronavirus disease, in particular COVID-19. The purpose of the present study was to find synergistic antiviral combinations for COVID-19 based on lethal mutagenesis.

EXPERIMENTAL APPROACH

The effect of combinations of remdesivir and ribavirin on the infectivity of SARS-CoV-2 in cell culture has been tested. Viral populations were monitored by ultra-deep sequencing, and the decrease of infectivity as a result of the treatment was measured.

KEY RESULTS

Remdesivir and ribavirin exerted a synergistic inhibitory activity against SARS-CoV-2, quantified both by CompuSyn (Chou-Talalay method) and Synergy Finder (ZIP-score model). In serial passage experiments, virus extinction was readily achieved with remdesivir-ribavirin combinations at concentrations well below their cytotoxic 50 value, but not with the drugs used individually. Deep sequencing of treated viral populations showed that remdesivir, ribavirin, and their combinations evoked significant increases of the number of viral mutations and haplotypes, as well as modification of diversity indices that characterize viral quasi-species.

CONCLUSION AND IMPLICATIONS

SARS-CoV-2 extinction can be achieved by synergistic combination treatments based on lethal mutagenesis. In addition, the results offer prospects of triple drug treatments for effective SARS-CoV-2 suppression.

SARS-CoV-2 population dynamics in immunocompetent individuals in a closed transmission chain shows genomic diversity over the course of infection

BACKGROUND

SARS-CoV-2 remains rapidly evolving, and many biologically important genomic substitutions/indels have characterised novel SARS-CoV-2 lineages, which have emerged during successive global waves of the pandemic. Worldwide genomic sequencing has been able to monitor these waves, track transmission clusters, and examine viral evolution in real time to help inform healthcare policy. One school of thought is that an apparent greater than average divergence in an emerging lineage from contemporary variants may require persistent infection, for example in an immunocompromised host. Due to the nature of the COVID-19 pandemic and sampling, there were few studies that examined the evolutionary trajectory of SARS-CoV-2 in healthy individuals.

METHODS

We investigated viral evolutionary trends and participant symptomatology within a cluster of 16 SARS-CoV-2 infected, immunocompetent individuals with no co-morbidities in a closed transmission chain. Longitudinal nasopharyngeal swab sampling allowed characterisation of SARS-CoV-2 intra-host variation over time at both the dominant and minor genomic variant levels through Nimagen-Illumina sequencing.

RESULTS

A change in viral lineage assignment was observed in individual infections; however, there was only one indel and no evidence of recombination over the period of an acute infection. Minor and dominant genomic modifications varied between participants, with some minor genomic modifications increasing in abundance to become the dominant viral sequence during infection.

CONCLUSIONS

Data from this cohort of SARS-CoV-2-infected participants demonstrated that long-term persistent infection in an immunocompromised host was not necessarily a prerequisite for generating a greater than average frequency of amino acid substitutions. Amino acid substitutions at both the dominant and minor genomic sequence level were observed in immunocompetent individuals during infection showing that viral lineage changes can occur generating viral diversity.

Transcriptome Analysis in Air-Liquid Interface Porcine Respiratory Epithelial Cell Cultures Reveals That the Betacoronavirus Porcine Encephalomyelitis Hemagglutinating Virus Induces a Robust Interferon Response to Infection

Porcine hemagglutinating encephalomyelitis virus (PHEV) replicates in the upper respiratory tract and tonsils of pigs. Using an air-liquid interface porcine respiratory epithelial cells (ALI-PRECs) culture system, we demonstrated that PHEV disrupts respiratory epithelia homeostasis by impairing ciliary function and inducing antiviral, pro-inflammatory cytokine, and chemokine responses. This study explores the mechanisms driving early innate immune responses during PHEV infection through host transcriptome analysis. Total RNA was collected from ALI-PRECs at 24, 36, and 48 h post inoculation (hpi). RNA-seq analysis was performed using an Illumina Hiseq 600 to generate 100 bp paired-end reads. Differential gene expression was analyzed using DeSeq2. PHEV replicated actively in ALI-PRECs, causing cytopathic changes and progressive mucociliary disruption. Transcriptome analysis revealed downregulation of cilia-associated genes such as CILK1, DNAH11, LRRC-23, -49, and -51, and acidic sialomucin CD164L2. PHEV also activated antiviral signaling pathways, significantly increasing the expression of interferon-stimulated genes (RSAD2, MX1, IFIT, and ISG15) and chemokine genes (CCL5 and CXCL10), highlighting inflammatory regulation. This study contributes to elucidating the molecular mechanisms of the innate immune response to PHEV infection of the airway epithelium, emphasizing the critical roles of the mucociliary, interferon, and chemokine responses.

Long COVID Definition, Symptoms, Risk Factors, Epidemiology and Autoimmunity: A Narrative Review

The virus called SARS-CoV-2 emerged in 2019 and quickly spread worldwide, causing COVID-19. It has greatly impacted on everyday life, healthcare systems, and the global economy. In order to save as many lives as possible, precautions such as social distancing, quarantine, and testing policies were implemented, and effective vaccines were developed. A growing amount of data collected worldwide allowed the characterization of this new disease, which turned out to be more complex than other common respiratory tract infections. An increasing number of convalescents presented with a variety of nonspecific symptoms emerging after the acute infection. This possible new global health problem was identified and labelled as long COVID. Since then, a great effort has been made by clinicians and the scientific community to understand the underlying mechanisms and to develop preventive measures and effective treatment. The role of autoimmunity induced by SARS-CoV-2 infection in the development of long COVID is discussed in this review. We aim to deliver a description of several conditions with an autoimmune background observed in COVID-19 convalescents, including Guillain-Barré syndrome, antiphospholipid syndrome and related thrombosis, and Kawasaki disease highlighting a relationship between SARS-CoV-2 infection and the development of autoimmunity. However, further studies are required to determine its true clinical significance.

Mapping of long stretches of highly conserved sequences in over 6 million SARS-CoV-2 genomes

We identified 11 conserved stretches in over 6.3 million SARS-CoV-2 genomes including all the major variants of concerns. Each conserved stretch is ≥100 nucleotides in length with ≥99.9% conservation at each nucleotide position. Interestingly, six of the eight conserved stretches in ORF1ab overlapped significantly with well-folded experimentally verified RNA secondary structures. Furthermore, two of the conserved stretches were mapped to regions within the S2-subunit that undergo dynamic structural rearrangements during viral fusion. In addition, the conserved stretches were significantly depleted for zinc-finger antiviral protein (ZAP) binding sites, which facilitated the recognition and degradation of viral RNA. These highly conserved stretches in the SARS-CoV-2 genome were poorly conserved at the nucleotide level among closely related β-coronaviruses, thus representing ideal targets for highly specific and discriminatory diagnostic assays. Our findings highlight the role of structural constraints at both RNA and protein levels that contribute to the sequence conservation of specific genomic regions in SARS-CoV-2.

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.

Sequence signatures within the genome of SARS-CoV-2 can be used to predict host source

We conducted an in silico analysis to better understand the potential factors impacting host adaptation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in white-tailed deer, humans, and mink due to the strong evidence of sustained transmission within these hosts. Classification models trained on single nucleotide and amino acid differences between samples effectively identified white-tailed deer-, human-, and mink-derived SARS-CoV-2. For example, the balanced accuracy score of Extremely Randomized Trees classifiers was 0.984 ± 0.006. Eighty-eight commonly identified predictive mutations are found at sites under strong positive and negative selective pressure. A large fraction of sites under selection (86.9%) or identified by machine learning (87.1%) are found in genes other than the spike. Some locations encoded by these gene regions are predicted to be B- and T-cell epitopes or are implicated in modulating the immune response suggesting that host adaptation may involve the evasion of the host immune system, modulation of the class-I major-histocompatibility complex, and the diminished recognition of immune epitopes by CD8+ T cells. Our selection and machine learning analysis also identified that silent mutations, such as C7303T and C9430T, play an important role in discriminating deer-derived samples across multiple clades. Finally, our investigation into the origin of the B.1.641 lineage from white-tailed deer in Canada discovered an additional human sequence from Michigan related to the B.1.641 lineage sampled near the emergence of this lineage. These findings demonstrate that machine-learning approaches can be used in combination with evolutionary genomics to identify factors possibly involved in the cross-species transmission of viruses and the emergence of novel viral lineages.IMPORTANCESevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a highly transmissible virus capable of infecting and establishing itself in human and wildlife populations, such as white-tailed deer. This fact highlights the importance of developing novel ways to identify genetic factors that contribute to its spread and adaptation to new host species. This is especially important since these populations can serve as reservoirs that potentially facilitate the re-introduction of new variants into human populations. In this study, we apply machine learning and phylogenetic methods to uncover biomarkers of SARS-CoV-2 adaptation in mink and white-tailed deer. We find evidence demonstrating that both non-synonymous and silent mutations can be used to differentiate animal-derived sequences from human-derived ones and each other. This evidence also suggests that host adaptation involves the evasion of the immune system and the suppression of antigen presentation. Finally, the methods developed here are general and can be used to investigate host adaptation in viruses other than SARS-CoV-2.

Multiple introductions of monkeypox virus to Ireland during the international mpox outbreak, May 2022 to October 2023

BackgroundMpox, caused by monkeypox virus (MPXV), was considered a rare zoonotic disease before May 2022, when a global epidemic of cases in non-endemic countries led to the declaration of a Public Health Emergency of International Concern. Cases of mpox in Ireland, a country without previous mpox reports, could reflect extended local transmission or multiple epidemiological introductions.AimTo elucidate the origins and molecular characteristics of MPXV circulating in Ireland between May 2022 and October 2023.MethodsWhole genome sequencing of MPXV from 75% of all Irish mpox cases (182/242) was performed and compared to sequences retrieved from public databases (n = 3,362). Bayesian approaches were used to infer divergence time between sequences from different subclades and evaluate putative importation events from other countries.ResultsOf 242 detected mpox cases, 99% were males (median age: 35 years; range: 15-60). All 182 analysed genomes were assigned to Clade IIb and, presence of 12 distinguishable subclades suggests multiple introductions into Ireland. Estimation of time to divergence of subclades further supports the hypothesis for multiple importation events from numerous countries, indicative of extended and sustained international spread of mpox. Further analysis of sequences revealed that 92% of nucleotide mutations were from cytosine to thymine (or from guanine to adenine), leading to a high number of non-synonymous mutations across subclades; mutations associated with tecovirimat resistance were not observed.ConclusionWe provide insights into the international transmission dynamics supporting multiple introductions of MPXV into Ireland. Such information supported the implementation of evidence-informed public health control measures.

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.

Mutational signature dynamics indicate SARS-CoV-2's evolutionary capacity is driven by host antiviral molecules

The COVID-19 pandemic has been characterised by sequential variant-specific waves shaped by viral, individual human and population factors. SARS-CoV-2 variants are defined by their unique combinations of mutations and there has been a clear adaptation to more efficient human infection since the emergence of this new human coronavirus in late 2019. Here, we use machine learning models to identify shared signatures, i.e., common underlying mutational processes and link these to the subset of mutations that define the variants of concern (VOCs). First, we examined the global SARS-CoV-2 genomes and associated metadata to determine how viral properties and public health measures have influenced the magnitude of waves, as measured by the number of infection cases, in different geographic locations using regression models. This analysis showed that, as expected, both public health measures and virus properties were associated with the waves of regional SARS-CoV-2 reported infection numbers and this impact varies geographically. We attribute this to intrinsic differences such as vaccine coverage, testing and sequencing capacity and the effectiveness of government stringency. To assess underlying evolutionary change, we used non-negative matrix factorisation and observed three distinct mutational signatures, unique in their substitution patterns and exposures from the SARS-CoV-2 genomes. Signatures 1, 2 and 3 were biased to C→T, T→C/A→G and G→T point mutations. We hypothesise assignments of these mutational signatures to the host antiviral molecules APOBEC, ADAR and ROS respectively. We observe a shift amidst the pandemic in relative mutational signature activity from predominantly Signature 1 changes to an increasingly high proportion of changes consistent with Signature 2. This could represent changes in how the virus and the host immune response interact and indicates how SARS-CoV-2 may continue to generate variation in the future. Linkage of the detected mutational signatures to the VOC-defining amino acids substitutions indicates the majority of SARS-CoV-2's evolutionary capacity is likely to be associated with the action of host antiviral molecules rather than virus replication errors.

Coronavirus takeover of host cell translation and intracellular antiviral response: a molecular perspective

Coronaviruses are a group of related RNA viruses that cause respiratory diseases in humans and animals. Understanding the mechanisms of translation regulation during coronaviral infections is critical for developing antiviral therapies and preventing viral spread. Translation of the viral single-stranded RNA genome in the host cell cytoplasm is an essential step in the life cycle of coronaviruses, which affects the cellular mRNA translation landscape in many ways. Here we discuss various viral strategies of translation control, including how members of the Betacoronavirus genus shut down host cell translation and suppress host innate immune functions, as well as the role of the viral non-structural protein 1 (Nsp1) in the process. We also outline the fate of viral RNA, considering stress response mechanisms triggered in infected cells, and describe how unique viral RNA features contribute to programmed ribosomal -1 frameshifting, RNA editing, and translation shutdown evasion.

Puzzles, challenges, and information reservoir of SARS-CoV-2 quasispecies

Upon the emergence of SARS-CoV-2 in the human population, it was conjectured that for this coronavirus the dynamic intra-host heterogeneity typical of RNA viruses would be toned down. Nothing of this sort is observed. Here we review the main observations on the complexity and diverse composition of SARS-CoV-2 mutant spectra sampled from infected patients, within the framework of quasispecies dynamics. The analyses suggest that the information provided by myriads of genomic sequences within infected individuals may have a predictive value of the genomic sequences that acquire epidemiological relevance. Possibilities to reconcile the presence of broad mutant spectra in the large RNA coronavirus genome with its encoding a 3' to 5' exonuclease proofreading-repair activity are considered. Indeterminations in the behavior of individual viral genomes provide a benefit for the survival of the ensemble. We propose that this concept falls in the domain of "stochastic thinking," a notion that applies also to cellular processes, as a means for biological systems to face unexpected needs.

Elevated A-to-I RNA editing in COVID-19 infected individuals

Given the current status of coronavirus disease 2019 (COVID-19) as a global pandemic, it is of high priority to gain a deeper understanding of the disease's development and how the virus impacts its host. Adenosine (A)-to-Inosine (I) RNA editing is a post-transcriptional modification, catalyzed by the ADAR family of enzymes, that can be considered part of the inherent cellular defense mechanism as it affects the innate immune response in a complex manner. It was previously reported that various viruses could interact with the host's ADAR enzymes, resulting in epigenetic changes both to the virus and the host. Here, we analyze RNA-seq of nasopharyngeal swab specimens as well as whole-blood samples of COVID-19 infected individuals and show a significant elevation in the global RNA editing activity in COVID-19 compared to healthy controls. We also detect specific coding sites that exhibit higher editing activity. We further show that the increment in editing activity during the disease is temporary and returns to baseline shortly after the symptomatic period. These significant epigenetic changes may contribute to the immune system response and affect adverse outcomes seen in post-viral cases.

Differential RNA editing landscapes in host cell versus the SARS-CoV-2 genome

The SARS-CoV-2 pandemic was defined by the emergence of new variants formed through virus mutation originating from random errors not corrected by viral proofreading and/or the host antiviral response introducing mutations into the viral genome. While sequencing information hints at cellular RNA editing pathways playing a role in viral evolution, here, we use an in vitro human cell infection model to assess RNA mutation types in two SARS-CoV-2 strains representing the original and the alpha variants. The variants showed both different cellular responses and mutation patterns with alpha showing higher mutation frequency with most substitutions observed being C-U, indicating an important role for apolipoprotein B mRNA editing catalytic polypeptide-like editing. Knockdown of select APOBEC3s through RNAi increased virus production in the original virus, but not in alpha. Overall, these data suggest a deaminase-independent anti-viral function of APOBECs in SARS-CoV-2 while the C-U editing itself might function to enhance genetic diversity enabling evolutionary adaptation.

Long-term dynamic shifts in genomic base content and evolutionary trajectories of SARS-CoV-2 variants

The rapid spread and remarkable mutations of SARS-CoV-2 variants, particularly Omicron, necessitate an understanding of their evolutionary characteristics. In this study, we analyzed representative high-quality whole-genome sequences of 2008 SARS-CoV-2 variants to explore long-term dynamic changes in genomic base (especially GC) content and variations during viral evolution. Our results demonstrated a highly negative correlation between GC content and variant emergence time (r = -0.765, p < 2.22e-16). Major gene partitions (S, N, ORF1ab) displayed similar trends. Omicron exhibited a significantly lower GC content than non-Omicron variants (p < 2.22e-16). Notably, we observed a robust negative correlation between C and T content (r = -0.778, p < 2.22e-16) and between G and A content (r = -0.773, p < 2.22e-16). Among all strains, Omicron showed the greatest base variation, with C->T mutations being the most frequent (median [interquartile range [IQR]]: 29 (27, 31), 37.67%), succeeded by G->A mutations (11 (9, 13), 14.63%). Over a 3-year span, an annual decline rate of 0.12% in SARS-CoV-2 GC content was observed and could become more pronounced in future emerging variants. These findings provided insights into the evolutionary trajectory of SARS-CoV-2, underscoring the significance of continuous genomic surveillance for effective prediction of and response to future variants.

SARS-CoV-2 within-host diversity of human hosts and its implications for viral immune evasion

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is continuously evolving, bringing great challenges to the control of the virus. In the present study, we investigated the characteristics of SARS-CoV-2 within-host diversity of human hosts and its implications for immune evasion using about 2,00,000 high-depth next-generation genome sequencing data of SARS-CoV-2. A total of 44% of the samples showed within-host variations (iSNVs), and the average number of iSNVs in the samples with iSNV was 1.90. C-to-U is the dominant substitution pattern for iSNVs. C-to-U/G-to-A and A-to-G/U-to-C preferentially occur in 5'-CG-3' and 5'-AU-3' motifs, respectively. In addition, we found that SARS-CoV-2 within-host variations are under negative selection. About 15.6% iSNVs had an impact on the content of the CpG dinucleotide (CpG) in SARS-CoV-2 genomes. We detected signatures of faster loss of CpG-gaining iSNVs, possibly resulting from zinc-finger antiviral protein-mediated antiviral activities targeting CpG, which could be the major reason for CpG depletion in SARS-CoV-2 consensus genomes. The non-synonymous iSNVs in the S gene can largely alter the S protein's antigenic features, and many of these iSNVs are distributed in the amino-terminal domain (NTD) and receptor-binding domain (RBD). These results suggest that SARS-CoV-2 interacts actively with human hosts and attempts to take different evolutionary strategies to escape human innate and adaptive immunity. These new findings further deepen and widen our understanding of the within-host evolutionary features of SARS-CoV-2. IMPORTANCE Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative pathogen of the coronavirus disease 2019, has evolved rapidly since it was discovered. Recent studies have pointed out that some mutations in the SARS-CoV-2 S protein could confer SARS-CoV-2 the ability to evade the human adaptive immune system. In addition, it is observed that the content of the CpG dinucleotide in SARS-CoV-2 genome sequences has decreased over time, reflecting the adaptation to the human host. The significance of our research is revealing the characteristics of SARS-CoV-2 within-host diversity of human hosts, identifying the causes of CpG depletion in SARS-CoV-2 consensus genomes, and exploring the potential impacts of non-synonymous within-host variations in the S gene on immune escape, which could further deepen and widen our understanding of the evolutionary features of SARS-CoV-2.

Distinguishing Genetic Drift from Selection in Papillomavirus Evolution

Pervasive purifying selection on non-synonymous substitutions is a hallmark of papillomavirus genome history, but the role of selection on and the drift of non-coding DNA motifs on HPV diversification is poorly understood. In this study, more than a thousand complete genomes representing Alphapapillomavirus types, lineages, and SNP variants were examined phylogenetically and interrogated for the number and position of non-coding DNA sequence motifs using Principal Components Analyses, Ancestral State Reconstructions, and Phylogenetic Independent Contrasts. For anciently diverged Alphapapillomavirus types, composition of the four nucleotides (A, C, G, T), codon usage, trimer usage, and 13 established non-coding DNA sequence motifs revealed phylogenetic clusters consistent with genetic drift. Ancestral state reconstruction and Phylogenetic Independent Contrasts revealed ancient genome alterations, particularly for the CpG and APOBEC3 motifs. Each evolutionary analytical method we performed supports the unanticipated conclusion that genetic drift and different evolutionary drivers have structured Alphapapillomavirus genomes in distinct ways during successive epochs, even extending to differences in more recently formed variant lineages.

SARS-CoV-2 entry into and evolution within a skilled nursing facility

SARS-CoV-2 belongs to the family Coronaviridae which includes multiple human pathogens that have an outsized impact on aging populations. As a novel human pathogen, SARS-CoV-2 is undergoing continuous adaptation to this new host species and there is evidence of this throughout the scientific and public literature. However, most investigations of SARS-CoV-2 evolution have focused on large-scale collections of data across diverse populations and/or living environments. Here we investigate SARS-CoV-2 evolution in epidemiologically linked individuals within a single outbreak at a skilled nursing facility beginning with initial introduction of the pathogen. The data demonstrate that SARS-CoV-2 was introduced to the facility multiple times without establishing an interfacility transmission chain, followed by a single introduction that infected many individuals within a week. This large-scale introduction by a single genotype then persisted in the facility. SARS-CoV-2 sequences were investigated at both the consensus and intra-host variation levels. Understanding the variability in SARS-CoV-2 during transmission chains will assist in understanding the spread of this disease and can ultimately inform best practices for mitigation strategies.

The importance of equally accessible genomic surveillance in the age of pandemics

Genomic epidemiology is now a core component in investigating the spread of a disease during an outbreak and for future preparedness to tackle emerging zoonoses. During the last decades, several viral diseases arose and emphasized the importance of molecular epidemiology in tracking the dispersal route, supporting proper mitigation measures, and appropriate vaccine development. In this perspective article, we summarized what has been done so far in the genomic epidemiology field and what should be considered in the future. We traced back the methods and protocols employed over time for zoonotic disease response. Either to small outbreaks such as the severe acute respiratory syndrome (SARS) outbreak identified first in 2002 in Guangdong, China, or to a global pandemic like the one that we are experiencing now since 2019 when the severe acute respiratory syndrome 2 (SARS-CoV-2) virus emerged in Wuhan, China, following several pneumonia cases, and subsequently spread worldwide. We explored both the benefits and shortages encountered when relying on genomic epidemiology, and we clearly present the disadvantages of inequity in accessing these tools around the world, especially in countries with less developed economies. For effectively addressing future pandemics, it is crucial to work for better sequencing equity around the globe.

Molecular transition of SARS-CoV-2 from critical patients during the first year of the COVID-19 pandemic in Mexico City

Background

The SARS-CoV-2 virus has caused unprecedented mortality since its emergence in late 2019. The continuous evolution of the viral genome through the concerted action of mutational forces has produced distinct variants that became dominant, challenging human immunity and vaccine development.

Aim and methods

In this work, through an integrative genomic approach, we describe the molecular transition of SARS-CoV-2 by analyzing the viral whole genome sequences from 50 critical COVID-19 patients recruited during the first year of the pandemic in Mexico City.

Results

Our results revealed differential levels of the evolutionary forces across the genome and specific mutational processes that have shaped the first two epidemiological waves of the pandemic in Mexico. Through phylogenetic analyses, we observed a genomic transition in the circulating SARS-CoV-2 genomes from several lineages prevalent in the first wave to a dominance of the B.1.1.519 variant (defined by T478K, P681H, and T732A mutations in the spike protein) in the second wave.

Conclusion

This work contributes to a better understanding of the evolutionary dynamics and selective pressures that act at the genomic level, the prediction of more accurate variants of clinical significance, and a better comprehension of the molecular mechanisms driving the evolution of SARS-CoV-2 to improve vaccine and drug development.

Host-mediated RNA editing in viruses

Viruses rely on hosts for life and reproduction, cause a variety of symptoms from common cold to AIDS to COVID-19 and provoke public health threats claiming millions of lives around the globe. RNA editing, as a crucial co-/post-transcriptional modification inducing nucleotide alterations on both endogenous and exogenous RNA sequences, exerts significant influences on virus replication, protein synthesis, infectivity and toxicity. Hitherto, a number of host-mediated RNA editing sites have been identified in diverse viruses, yet lacking a full picture of RNA editing-associated mechanisms and effects in different classes of viruses. Here we synthesize the current knowledge of host-mediated RNA editing in a variety of viruses by considering two enzyme families, viz., ADARs and APOBECs, thereby presenting a landscape of diverse editing mechanisms and effects between viruses and hosts. In the ongoing pandemic, our study promises to provide potentially valuable insights for better understanding host-mediated RNA editing on ever-reported and newly-emerging viruses.

L-shaped distribution of the relative substitution rate (c/μ) observed for SARS-COV-2's genome, inconsistent with the selectionist theory, the neutral theory and the nearly neutral theory but a near-neutral balanced selection theory: Implication on "neutralist-selectionist" debate

The genomic substitution rate (GSR) of SARS-CoV-2 exhibits a molecular clock feature and does not change under fluctuating environmental factors such as the infected human population (10°-107), vaccination etc. The molecular clock feature is believed to be inconsistent with the selectionist theory (ST). The GSR shows lack of dependence on the effective population size, suggesting Ohta's nearly neutral theory (ONNT) is not applicable to this virus. Big variation of the substitution rate within its genome is also inconsistent with Kimura's neutral theory (KNT). Thus, all three existing evolution theories fail to explain the evolutionary nature of this virus. In this paper, we proposed a Segment Substitution Rate Model (SSRM) under non-neutral selections and pointed out that a balanced mechanism between negative and positive selection of some segments that could also lead to the molecular clock feature. We named this hybrid mechanism as near-neutral balanced selection theory (NNBST) and examined if it was followed by SARS-CoV-2 using the three independent sets of SARS-CoV-2 genomes selected by the Nextstrain team. Intriguingly, the relative substitution rate of this virus exhibited an L-shaped probability distribution consisting with NNBST rather than Poisson distribution predicted by KNT or an asymmetric distribution predicted by ONNT in which nearly neutral sites are believed to be slightly deleterious only, or the distribution that is lack of nearly neutral sites predicted by ST. The time-dependence of the substitution rates for some segments and their correlation with the vaccination were observed, supporting NNBST. Our relative substitution rate method provides a tool to resolve the long standing "neutralist-selectionist" controversy. Implications of NNBST in resolving Lewontin's Paradox is also discussed.

Correlated substitutions reveal SARS-like coronaviruses recombine frequently with a diverse set of structured gene pools

Quantifying SARS-like coronavirus (SL-CoV) evolution is critical to understanding the origins of SARS-CoV-2 and the molecular processes that could underlie future epidemic viruses. While genomic analyses suggest recombination was a factor in the emergence of SARS-CoV-2, few studies have quantified recombination rates among SL-CoVs. Here, we infer recombination rates of SL-CoVs from correlated substitutions in sequencing data using a coalescent model with recombination. Our computationally-efficient, non-phylogenetic method infers recombination parameters of both sampled sequences and the unsampled gene pools with which they recombine. We apply this approach to infer recombination parameters for a range of positive-sense RNA viruses. We then analyze a set of 191 SL-CoV sequences (including SARS-CoV-2) and find that ORF1ab and S genes frequently undergo recombination. We identify which SL-CoV sequence clusters have recombined with shared gene pools, and show that these pools have distinct structures and high recombination rates, with multiple recombination events occurring per synonymous substitution. We find that individual genes have recombined with different viral reservoirs. By decoupling contributions from mutation and recombination, we recover the phylogeny of non-recombined portions for many of these SL-CoVs, including the position of SARS-CoV-2 in this clonal phylogeny. Lastly, by analyzing >400,000 SARS-CoV-2 whole genome sequences, we show current diversity levels are insufficient to infer the within-population recombination rate of the virus since the pandemic began. Our work offers new methods for inferring recombination rates in RNA viruses with implications for understanding recombination in SARS-CoV-2 evolution and the structure of clonal relationships and gene pools shaping its origins.

Atypical Mutational Spectrum of SARS-CoV-2 Replicating in the Presence of Ribavirin

We report that ribavirin exerts an inhibitory and mutagenic activity on SARS-CoV-2-infecting Vero cells, with a therapeutic index higher than 10. Deep sequencing analysis of the mutant spectrum of SARS-CoV-2 replicating in the absence or presence of ribavirin indicated an increase in the number of mutations, but not in deletions, and modification of diversity indices, expected from a mutagenic activity. Notably, the major mutation types enhanced by replication in the presence of ribavirin were A→G and U→C transitions, a pattern which is opposite to the dominance of G→A and C→U transitions previously described for most RNA viruses. Implications of the inhibitory activity of ribavirin, and the atypical mutational bias produced on SARS-CoV-2, for the search for synergistic anti-COVID-19 lethal mutagen combinations are discussed.

Mutational Patterns Observed in SARS-CoV-2 Genomes Sampled From Successive Epochs Delimited by Major Public Health Events in Ontario, Canada: Genomic Surveillance Study

BACKGROUND

The emergence of SARS-CoV-2 variants with mutations associated with increased transmissibility and virulence is a public health concern in Ontario, Canada. Characterizing how the mutational patterns of the SARS-CoV-2 genome have changed over time can shed light on the driving factors, including selection for increased fitness and host immune response, that may contribute to the emergence of novel variants. Moreover, the study of SARS-CoV-2 in the microcosm of Ontario, Canada can reveal how different province-specific public health policies over time may be associated with observed mutational patterns as a model system.

OBJECTIVE

This study aimed to perform a comprehensive analysis of single base substitution (SBS) types, counts, and genomic locations observed in SARS-CoV-2 genomic sequences sampled in Ontario, Canada. Comparisons of mutational patterns were conducted between sequences sampled during 4 different epochs delimited by major public health events to track the evolution of the SARS-CoV-2 mutational landscape over 2 years.

METHODS

In total, 24,244 SARS-CoV-2 genomic sequences and associated metadata sampled in Ontario, Canada from January 1, 2020, to December 31, 2021, were retrieved from the Global Initiative on Sharing All Influenza Data database. Sequences were assigned to 4 epochs delimited by major public health events based on the sampling date. SBSs from each SARS-CoV-2 sequence were identified relative to the MN996528.1 reference genome. Catalogues of SBS types and counts were generated to estimate the impact of selection in each open reading frame, and identify mutation clusters. The estimation of mutational fitness over time was performed using the Augur pipeline.

RESULTS

The biases in SBS types and proportions observed support previous reports of host antiviral defense activity involving the SARS-CoV-2 genome. There was an increase in U>C substitutions associated with adenosine deaminase acting on RNA (ADAR) activity uniquely observed during Epoch 4. The burden of novel SBSs observed in SARS-CoV-2 genomic sequences was the greatest in Epoch 2 (median 5), followed by Epoch 3 (median 4). Clusters of SBSs were observed in the spike protein open reading frame, ORF1a, and ORF3a. The high proportion of nonsynonymous SBSs and increasing dN/dS metric (ratio of nonsynonymous to synonymous mutations in a given open reading frame) to above 1 in Epoch 4 indicate positive selection of the spike protein open reading frame.

CONCLUSIONS

Quantitative analysis of the mutational patterns of the SARS-CoV-2 genome in the microcosm of Ontario, Canada within early consecutive epochs of the pandemic tracked the mutational dynamics in the context of public health events that instigate significant shifts in selection and mutagenesis. Continued genomic surveillance of emergent variants will be useful for the design of public health policies in response to the evolving COVID-19 pandemic.

Divergent SARS-CoV-2 variant emerges in white-tailed deer with deer-to-human transmission

Wildlife reservoirs of broad-host-range viruses have the potential to enable evolution of viral variants that can emerge to infect humans. In North America, there is phylogenomic evidence of continual transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) from humans to white-tailed deer (Odocoileus virginianus) through unknown means, but no evidence of transmission from deer to humans. We carried out an observational surveillance study in Ontario, Canada during November and December 2021 (n = 300 deer) and identified a highly divergent lineage of SARS-CoV-2 in white-tailed deer (B.1.641). This lineage is one of the most divergent SARS-CoV-2 lineages identified so far, with 76 mutations (including 37 previously associated with non-human mammalian hosts). From a set of five complete and two partial deer-derived viral genomes we applied phylogenomic, recombination, selection and mutation spectrum analyses, which provided evidence for evolution and transmission in deer and a shared ancestry with mink-derived virus. Our analysis also revealed an epidemiologically linked human infection. Taken together, our findings provide evidence for sustained evolution of SARS-CoV-2 in white-tailed deer and of deer-to-human transmission.

The mutational spectrum of SARS-CoV-2 genomic and antigenomic RNA

The raw material for viral evolution is provided by intra-host mutations occurring during replication, transcription or post-transcription. Replication and transcription of Coronaviridae proceed through the synthesis of negative-sense 'antigenomes' acting as templates for positive-sense genomic and subgenomic RNA. Hence, mutations in the genomes of SARS-CoV-2 and other coronaviruses can occur during (and after) the synthesis of either negative-sense or positive-sense RNA, with potentially distinct patterns and consequences. We explored for the first time the mutational spectrum of SARS-CoV-2 (sub)genomic and anti(sub)genomic RNA. We use a high-quality deep sequencing dataset produced using a quantitative strand-aware sequencing method, controlled for artefacts and sequencing errors, and scrutinized for accurate detection of within-host diversity. The nucleotide differences between negative- and positive-sense strand consensus vary between patients and do not show dependence on age or sex. Similarities and differences in mutational patterns between within-host minor variants on the two RNA strands suggested strand-specific mutations or editing by host deaminases and oxidative damage. We observe generally neutral and slight negative selection on the negative strand, contrasting with purifying selection in ORF1a, ORF1b and S genes of the positive strand of the genome.

Advances in Next-Generation Coronavirus Vaccines in Response to Future Virus Evolution

Coronavirus disease (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has spread to more than 230 countries and territories worldwide since its outbreak in late 2019. In less than three years, infection by SARS-CoV-2 has resulted in over 600 million cases of COVID-19 and over 6.4 million deaths. Vaccines have been developed with unimaginable speed, and 11 have already been approved by the World Health Organization and given Emergency Use Listing. The administration of several first-generation SARS-CoV-2 vaccines has successfully decelerated the spread of COVID-19 but not stopped it completely. In the ongoing fight against viruses, genetic mutations frequently occur in the viral genome, resulting in a decrease in vaccine-induced antibody neutralization and widespread breakthrough infection. Facing the evolution and uncertainty of SARS-CoV-2 in the future, and the possibility of the spillover of other coronaviruses to humans, the need for vaccines with a broad spectrum of antiviral variants against multiple coronaviruses is recognized. It is imperative to develop a universal coronavirus or pan-coronavirus vaccine or drug to combat the ongoing COVID-19 pandemic as well as to prevent the next coronavirus pandemic. In this review, in addition to summarizing the protective effect of approved vaccines, we systematically summarize current work on the development of vaccines aimed at suppressing multiple SARS-CoV-2 variants of concern as well as multiple coronaviruses.

Retasking of canonical antiviral factors into proviral effectors

Under constant barrage by viruses, hosts have evolved a plethora of antiviral effectors and defense mechanisms. To survive, viruses must adapt to evade or subvert these defenses while still capturing cellular resources to fuel their replication cycles. Large-scale studies of the antiviral activities of cellular proteins and processes have shown that different viruses are controlled by distinct subsets of antiviral genes. The remaining antiviral genes are either ineffective in controlling infection, or in some cases, actually promote infection. In these cases, classically defined antiviral factors are retasked by viruses to enhance viral replication. This creates a more nuanced picture revealing the contextual nature of antiviral activity. The same protein can exert different effects on replication, depending on multiple factors, including the host, the target cells, and the specific virus infecting it. Here, we review numerous examples of viruses hijacking canonically antiviral proteins and retasking them for proviral purposes.

The roles of APOBEC-mediated RNA editing in SARS-CoV-2 mutations, replication and fitness

During COVID-19 pandemic, mutations of SARS-CoV-2 produce new strains that can be more infectious or evade vaccines. Viral RNA mutations can arise from misincorporation by RNA-polymerases and modification by host factors. Analysis of SARS-CoV-2 sequence from patients showed a strong bias toward C-to-U mutation, suggesting a potential mutational role by host APOBEC cytosine deaminases that possess broad anti-viral activity. We report the first experimental evidence demonstrating that APOBEC3A, APOBEC1, and APOBEC3G can edit on specific sites of SARS-CoV-2 RNA to produce C-to-U mutations. However, SARS-CoV-2 replication and viral progeny production in Caco-2 cells are not inhibited by the expression of these APOBECs. Instead, expression of wild-type APOBEC3 greatly promotes viral replication/propagation, suggesting that SARS-CoV-2 utilizes the APOBEC-mediated mutations for fitness and evolution. Unlike the random mutations, this study suggests the predictability of all possible viral genome mutations by these APOBECs based on the UC/AC motifs and the viral genomic RNA structure.

Phylogenetic and phylodynamic approaches to understanding and combating the early SARS-CoV-2 pandemic

Determining the transmissibility, prevalence and patterns of movement of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is central to our understanding of the impact of the pandemic and to the design of effective control strategies. Phylogenies (evolutionary trees) have provided key insights into the international spread of SARS-CoV-2 and enabled investigation of individual outbreaks and transmission chains in specific settings. Phylodynamic approaches combine evolutionary, demographic and epidemiological concepts and have helped track virus genetic changes, identify emerging variants and inform public health strategy. Here, we review and synthesize studies that illustrate how phylogenetic and phylodynamic techniques were applied during the first year of the pandemic, and summarize their contributions to our understanding of SARS-CoV-2 transmission and control.

Massive-scale genomic analysis reveals SARS-CoV-2 mutation characteristics and evolutionary trends

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic resulted in significant societal costs. Hence, an in-depth understanding of SARS-CoV-2 virus mutation and its evolution will help determine the direction of the COVID-19 pandemic. In this study, we identified 296,728 de novo mutations in more than 2,800,000 high-quality SARS-CoV-2 genomes. All possible factors affecting the mutation frequency of SARS-CoV-2 in human hosts were analyzed, including zinc finger antiviral proteins, sequence context, amino acid change, and translation efficiency. As a result, we proposed that when adenine (A) and tyrosine (T) bases are in the context of AM (M stands for adenine or cytosine) or TA motif, A or T base has lower mutation frequency. Furthermore, we hypothesized that translation efficiency can affect the mutation frequency of the third position of the codon by the selection, which explains why SARS-CoV-2 prefers AT3 codons usage. In addition, we found a host-specific asymmetric dinucleotide mutation frequency in the SARS-CoV-2 genome, which provides a new basis for determining the origin of the SARS-CoV-2. Finally, we summarize all possible factors affecting mutation frequency and provide insights into the mutation characteristics and evolutionary trends of SARS-CoV-2.

Heterologous immunity induced by 1st generation COVID-19 vaccines and its role in developing a pan-coronavirus vaccine

Severe acute respiratory syndrome virus-2 (SARS-CoV-2), the causative infectious agent of the COVID-19 pandemic, has led to multiple (4-6) waves of infections worldwide during the past two years. The development of vaccines against SARS-CoV-2 has led to successful mass immunizations worldwide, mitigating the worldwide mortality due the pandemic to a great extent. Yet the evolution of new variants highlights a need to develop a universal vaccine which can prevent infections from all virulent SARS-CoV-2. Most of the current first generation COVID-19 vaccines are based on the Spike protein from the original Wuhan-hu-1 virus strain. It is encouraging that they still protect from serious illnesses, hospitalizations and mortality against a number of mutated viral strains, to varying degrees. Understanding the mechanisms by which these vaccines provide heterologous protection against multiple highly mutated variants can reveal strategies to develop a universal vaccine. In addition, many unexposed individuals have been found to harbor T cells that are cross-reactive against SARS-CoV-2 antigens, with a possible protective role. In this review, we will discuss various aspects of natural or vaccine-induced heterologous (cross-reactive) adaptive immunity against SARS-CoV-2 and other coronaviruses, and their role in achieving the concept of a pan-coronavirus vaccine.

The SARS-CoV-2 differential genomic adaptation in response to varying UVindex reveals potential genomic resources for better COVID-19 diagnosis and prevention

Coronavirus disease 2019 (COVID-19) has been a pandemic disease reported in almost every country and causes life-threatening, severe respiratory symptoms. Recent studies showed that various environmental selection pressures challenge the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infectivity and, in response, the virus engenders new mutations, leading to the emergence of more virulent strains of WHO concern. Advance prediction of the forthcoming virulent SARS-CoV-2 strains in response to the principal environmental selection pressures like temperature and solar UV radiation is indispensable to overcome COVID-19. To discover the UV-solar radiation-driven genomic adaption of SARS-CoV-2, a curated dataset of 2,500 full-grade genomes from five different UVindex regions (25 countries) was subjected to in-depth downstream genome-wide analysis. The recurrent variants that best respond to UV-solar radiations were extracted and extensively annotated to determine their possible effects and impacts on gene functions. This study revealed 515 recurrent single nucleotide variants (rcntSNVs) as SARS-CoV-2 genomic responses to UV-solar radiation, of which 380 were found to be distinct. For all discovered rcntSNVs, 596 functional effects (rcntEffs) were detected, containing 290 missense, 194 synonymous, 81 regulatory, and 31 in the intergenic region. The highest counts of missense rcntSNVs in spike (27) and nucleocapsid (26) genes explain the SARS-CoV-2 genomic adjustment to escape immunity and prevent UV-induced DNA damage, respectively. Among all, the most commonly observed rcntEffs were four missenses (RdRp-Pro327Leu, N-Arg203Lys, N-Gly204Arg, and Spike-Asp614Gly) and one synonymous (ORF1ab-Phe924Phe) functional effects. The highest number of rcntSNVs found distinct and were uniquely attributed to the specific UVindex regions, proposing solar-UV radiation as one of the driving forces for SARS-CoV-2 differential genomic adaptation. The phylogenetic relationship indicated the high UVindex region populating SARS-CoV-2 as the recent progenitor of all included samples. Altogether, these results provide baseline genomic data that may need to be included for preparing UVindex region-specific future diagnostic and vaccine formulations.

The past, current and future epidemiological dynamic of SARS-CoV-2

SARS-CoV-2, the agent of the COVID-19 pandemic, emerged in late 2019 in China, and rapidly spread throughout the world to reach all continents. As the virus expanded in its novel human host, viral lineages diversified through the accumulation of around two mutations a month on average. Different viral lineages have replaced each other since the start of the pandemic, with the most successful Alpha, Delta and Omicron variants of concern (VoCs) sequentially sweeping through the world to reach high global prevalence. Neither Alpha nor Delta was characterized by strong immune escape, with their success coming mainly from their higher transmissibility. Omicron is far more prone to immune evasion and spread primarily due to its increased ability to (re-)infect hosts with prior immunity. As host immunity reaches high levels globally through vaccination and prior infection, the epidemic is expected to transition from a pandemic regime to an endemic one where seasonality and waning host immunization are anticipated to become the primary forces shaping future SARS-CoV-2 lineage dynamics. In this review, we consider a body of evidence on the origins, host tropism, epidemiology, genomic and immunogenetic evolution of SARS-CoV-2 including an assessment of other coronaviruses infecting humans. Considering what is known so far, we conclude by delineating scenarios for the future dynamic of SARS-CoV-2, ranging from the good-circulation of a fifth endemic 'common cold' coronavirus of potentially low virulence, the bad-a situation roughly comparable with seasonal flu, and the ugly-extensive diversification into serotypes with long-term high-level endemicity.

Emerging SARS-CoV-2 variants: Why, how, and what's next?

The emergence of the SARS-CoV-2 Omicron variant poses a striking threat to human society. More than 30 mutations in the Spike protein of the Omicron variant severely compromised the protective immunity elicited by either vaccination or prior infection. The persistent viral evolutionary trajectory generates Omicron-associated lineages, such as BA.1 and BA.2. Moreover, the virus recombination upon Delta and Omicron co-infections has been reported lately, although the impact remains to be assessed. This minireview summarizes the characteristics, evolution and mutation control, and immune evasion mechanisms of SARS-CoV-2 variants, which will be helpful for the in-depth understanding of the SARS-CoV-2 variants and policy-making related to COVID-19 pandemic control.

Transmission of SARS-CoV-2 from humans to animals and potential host adaptation

SARS-CoV-2, the causative agent of the COVID-19 pandemic, can infect a wide range of mammals. Since its spread in humans, secondary host jumps of SARS-CoV-2 from humans to multiple domestic and wild populations of mammals have been documented. Understanding the extent of adaptation to these animal hosts is critical for assessing the threat that the spillback of animal-adapted SARS-CoV-2 into humans poses. We compare the genomic landscapes of SARS-CoV-2 isolated from animal species to that in humans, profiling the mutational biases indicative of potentially different selective pressures in animals. We focus on viral genomes isolated from mink (Neovison vison) and white-tailed deer (Odocoileus virginianus) for which multiple independent outbreaks driven by onward animal-to-animal transmission have been reported. We identify five candidate mutations for animal-specific adaptation in mink (NSP9_G37E, Spike_F486L, Spike_N501T, Spike_Y453F, ORF3a_L219V), and one in deer (NSP3a_L1035F), though they appear to confer a minimal advantage for human-to-human transmission. No considerable changes to the mutation rate or evolutionary trajectory of SARS-CoV-2 has resulted from circulation in mink and deer thus far. Our findings suggest that minimal adaptation was required for onward transmission in mink and deer following human-to-animal spillover, highlighting the 'generalist' nature of SARS-CoV-2 as a mammalian pathogen.

Commentary on "Poor evidence for host-dependent regular RNA editing in the transcriptome of SARS-CoV-2"

Analysis of the SARS-CoV-2 transcriptome has revealed a background of low-frequency intra-host genetic changes with a strong bias towards transitions. A similar pattern is also observed when inter-host variability is considered. We and others have shown that the cellular RNA editing machinery based on ADAR and APOBEC host-deaminases could be involved in the onset of SARS-CoV-2 genetic variability. Our hypothesis is based both on similarities with other known forms of viral genome editing and on the excess of transition changes, which is difficult to explain with errors during viral replication. Zong et al. criticize our analysis on both conceptual and technical grounds. While ultimate proof of an involvement of host deaminases in viral RNA editing will depend on experimental validation, here, we address the criticism to suggest that viral RNA editing is the most reasonable explanation for the observed intra- and inter-host variability.

SARS-CoV-2 Point Mutation and Deletion Spectra and Their Association with Different Disease Outcomes

Mutant spectra of RNA viruses are important to understand viral pathogenesis and response to selective pressures. There is a need to characterize the complexity of mutant spectra in coronaviruses sampled from infected patients. In particular, the possible relationship between SARS-CoV-2 mutant spectrum complexity and disease associations has not been established. In the present study, we report an ultradeep sequencing (UDS) analysis of the mutant spectrum of amplicons from the nsp12 (polymerase)- and spike (S)-coding regions of 30 nasopharyngeal isolates (diagnostic samples) of SARS-CoV-2 of the first COVID-19 pandemic wave (Madrid, Spain, April 2020) classified according to the severity of ensuing COVID-19. Low-frequency mutations and deletions, counted relative to the consensus sequence of the corresponding isolate, were overwhelmingly abundant. We show that the average number of different point mutations, mutations per haplotype, and several diversity indices was significantly higher in SARS-CoV-2 isolated from patients who developed mild disease than in those associated with moderate or severe disease (exitus). No such bias was observed with RNA deletions. Location of amino acid substitutions in the three-dimensional structures of nsp12 (polymerase) and S suggest significant structural or functional effects. Thus, patients who develop mild symptoms may be a richer source of genetic variants of SARS-CoV-2 than patients with moderate or severe COVID-19. IMPORTANCE The study shows that mutant spectra of SARS-CoV-2 from diagnostic samples differ in point mutation abundance and complexity and that significantly larger values were observed in virus from patients who developed mild COVID-19 symptoms. Mutant spectrum complexity is not a uniform trait among isolates. The nature and location of low-frequency amino acid substitutions present in mutant spectra anticipate great potential for phenotypic diversification of SARS-CoV-2.

1H, 13C and 15N chemical shift assignment of the stem-loops 5b + c from the 5'-UTR of SARS-CoV-2

The ongoing pandemic of the respiratory disease COVID-19 is caused by the SARS-CoV-2 (SCoV2) virus. SCoV2 is a member of the Betacoronavirus genus. The 30 kb positive sense, single stranded RNA genome of SCoV2 features 5'- and 3'-genomic ends that are highly conserved among Betacoronaviruses. These genomic ends contain structured cis-acting RNA elements, which are involved in the regulation of viral replication and translation. Structural information about these potential antiviral drug targets supports the development of novel classes of therapeutics against COVID-19. The highly conserved branched stem-loop 5 (SL5) found within the 5'-untranslated region (5'-UTR) consists of a basal stem and three stem-loops, namely SL5a, SL5b and SL5c. Both, SL5a and SL5b feature a 5'-UUUCGU-3' hexaloop that is also found among Alphacoronaviruses. Here, we report the extensive ¹H, ¹³C and ¹⁵N resonance assignment of the 37 nucleotides (nts) long sequence spanning SL5b and SL5c (SL5b + c), as basis for further in-depth structural studies by solution NMR spectroscopy.

Statistical modeling of SARS-CoV-2 substitution processes: predicting the next variant

We build statistical models to describe the substitution process in the SARS-CoV-2 as a function of explanatory factors describing the sequence, its function, and more. These models serve two different purposes: first, to gain knowledge about the evolutionary biology of the virus; and second, to predict future mutations in the virus, in particular, non-synonymous amino acid substitutions creating new variants. We use tens of thousands of publicly available SARS-CoV-2 sequences and consider tens of thousands of candidate models. Through a careful validation process, we confirm that our chosen models are indeed able to predict new amino acid substitutions: candidates ranked high by our model are eight times more likely to occur than random amino acid changes. We also show that named variants were highly ranked by our models before their appearance, emphasizing the value of our models for identifying likely variants and potentially utilizing this knowledge in vaccine design and other aspects of the ongoing battle against COVID-19.

As the virus that causes COVID-19 continues to mutate and spread, new methods are needed to predict new potential variants. Here, the authors identify the best regression models for predicting likely mutation sites in the SARS-CoV-2 genome using a candidate set that considers sequence, gene location, and biological function.

Omicron: What Makes the Latest SARS-CoV-2 Variant of Concern So Concerning?

Emerging strains of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of the coronavirus disease 2019 (COVID-19) pandemic, that show increased transmission fitness and/or immune evasion are classified as "variants of concern" (VOCs). Recently, a SARS-CoV-2 variant first identified in November 2021 in South Africa has been recognized as a fifth VOC, termed "Omicron." What makes this VOC so alarming is the high number of changes, especially in the viral Spike protein, and accumulating evidence for increased transmission efficiency and escape from neutralizing antibodies. In an amazingly short time, the Omicron VOC has outcompeted the previously dominating Delta VOC. However, it seems that the Omicron VOC is overall less pathogenic than other SARS-CoV-2 VOCs. Here, we provide an overview of the mutations in the Omicron genome and the resulting changes in viral proteins compared to other SARS-CoV-2 strains and discuss their potential functional consequences.

SARS-CoV-2 Evolution: On the Sudden Appearance of the Omicron Variant

The Omicron variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spreads rapidly and harbors many mutations in the spike protein, but the origin of this virus variant remains unclear. We address the role of unusual virus evolution mechanisms such as hypermutation, out-of-frame reading, and recombination. Rather, regular Darwinian evolution, that is, the repeated selection of beneficial spike mutations, seems to have led to the appearance of the grossly altered spike protein of the Omicron variant.

Functions and consequences of AID/APOBEC-mediated DNA and RNA deamination

The AID/APOBEC polynucleotide cytidine deaminases have historically been classified as either DNA mutators or RNA editors based on their first identified nucleic acid substrate preference. DNA mutators can generate functional diversity at antibody genes but also cause genomic instability in cancer. RNA editors can generate informational diversity in the transcriptome of innate immune cells, and of cancer cells. Members of both classes can act as antiviral restriction factors. Recent structural work has illuminated differences and similarities between AID/APOBEC enzymes that can catalyse DNA mutation, RNA editing or both, suggesting that the strict functional classification of members of this family should be reconsidered. As many of these enzymes have been employed for targeted genome (or transcriptome) editing, a more holistic understanding will help improve the design of therapeutically relevant programmable base editors.

In this Perspective, Pecori et al. provide an overview of the AID/APOBEC cytidine deaminase family, discussing key structural features, how they contribute to viral and tumour evolution and how they can be harnessed for (potentially therapeutic) base-editing purposes.

Whole-Genome Sequencing of SARS-CoV-2 Infection in a Cluster of Immunocompromised Children in Indonesia

Background

Thus far, Indonesia has recorded over 4,000,000 confirmed COVID-19 cases and 144,000 fatalities; 12.8% of cases have been in children under 18 years. Whole-genome viral sequencing (WGS) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been demonstrated to help differentiate hospital-acquired infection from community-acquired coronavirus disease 2019 (COVID-19) infection. Our study highlighted the use of WGS to investigate the origin of infection among pediatric oncology patients in Jakarta. The aim of our study was to evaluate clinical and laboratory characteristics and also the efficacy of using WGS to confirm hospital-acquired COVID-19 infection in a cluster of immunocompromised children within a single ward of a tertiary hospital in metropolitan Jakarta based on quasispecies, viral load, and admission dates.

Method

Real-time reverse-transcription polymerase chain reaction (RT-PCR) from nasopharyngeal (NP) swabs was used to diagnose the patients and also guardians and healthcare workers (HCWs) in the ward, followed by WGS of RT-PCR positive cases to establish their phylogenetic relationships.

Result

Using WGS, we showed that SARS-CoV-2 transmission in a cluster of children with underlying malignancy was characterized by high similarity of whole virus genome, which suggests nosocomial transmission.

Inosine and its methyl derivatives: Occurrence, biogenesis, and function in RNA

Inosine is one of the most common post-transcriptional modifications. Since its discovery, it has been noted for its ability to contribute to non-Watson-Crick interactions within RNA. Rapidly accumulating evidence points to the widespread generation of inosine through hydrolytic deamination of adenosine to inosine by different classes of adenosine deaminases. Three naturally occurring methyl derivatives of inosine, i.e., 1-methylinosine, 2'-O-methylinosine and 1,2'-O-dimethylinosine are currently reported in RNA modification databases. These modifications are expected to lead to changes in the structure, folding, dynamics, stability and functions of RNA. The importance of the modifications is indicated by the strong conservation of the modifying enzymes across organisms. The structure, binding and catalytic mechanism of the adenosine deaminases have been well-studied, but the underlying mechanism of the catalytic reaction is not very clear yet. Here we extensively review the existing data on the occurrence, biogenesis and functions of inosine and its methyl derivatives in RNA. We also included the structural and thermodynamic aspects of these modifications in our review to provide a detailed and integrated discussion on the consequences of A-to-I editing in RNA and the contribution of different structural and thermodynamic studies in understanding its role in RNA. We also highlight the importance of further studies for a better understanding of the mechanisms of the different classes of deamination reactions. Further investigation of the structural and thermodynamic consequences and functions of these modifications in RNA should provide more useful information about their role in different diseases.

Impact of ADAR-induced editing of minor viral RNA populations on replication and transmission of SARS-CoV-2

Viral RNA may be edited by enzymes of the ADAR family that deaminate adenosine residues with ensuing A→G mutations. We found multiple A→G mutations in minor viral populations of the SARS-CoV-2 genome. A→G mutations accumulated in the receptor binding domain of the spike gene, which may cause structural changes by altering binding to the ACE2 receptor. Presence of A→G mutations in minor viral populations was associated with reduced viral load, implying that ADAR may limit viral replication. Analyses of >250,000 European samples from 2020 revealed that A→G mutations in SARS-CoV-2 RNA were inversely correlated with mortality as a reflection of incidence. ADAR may thus be important in providing new variants of SARS-CoV-2 with altered infectivity and transmissibility.

Adenosine deaminases acting on RNA (ADAR) are RNA-editing enzymes that may restrict viral infection. We have utilized deep sequencing to determine adenosine to guanine (A→G) mutations, signifying ADAR activity, in clinical samples retrieved from 93 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)–infected patients in the early phase of the COVID-19 pandemic. A→G mutations were detected in 0.035% (median) of RNA residues and were predominantly nonsynonymous. These mutations were rarely detected in the major viral population but were abundant in minor viral populations in which A→G was more prevalent than any other mutation (P < 0.001). The A→G substitutions accumulated in the spike protein gene at positions corresponding to amino acids 505 to 510 in the receptor binding motif and at amino acids 650 to 655. The frequency of A→G mutations in minor viral populations was significantly associated with low viral load (P < 0.001). We additionally analyzed A→G mutations in 288,247 SARS-CoV-2 major (consensus) sequences representing the dominant viral population. The A→G mutations observed in minor viral populations in the initial patient cohort were increasingly detected in European consensus sequences between March and June 2020 (P < 0.001) followed by a decline of these mutations in autumn and early winter (P < 0.001). We propose that ADAR-induced deamination of RNA is a significant source of mutated SARS-CoV-2 and hypothesize that the degree of RNA deamination may determine or reflect viral fitness and infectivity.