Extensive C->U transition biases in the genomes of a wide range of mammalian RNA viruses; potential associations with transcriptional mutations, damage- or host-mediated editing of viral RNA

Mapping of susceptibility loci for Ebola virus pathogenesis in mice

Ebola virus (EBOV), a major global health concern, causes severe, often fatal EBOV disease (EVD) in humans. Host genetic variation plays a critical role, yet the identity of host susceptibility loci in mammals remains unknown. Using genetic reference populations, we generate an F2 mapping cohort to identify host susceptibility loci that regulate EVD. While disease-resistant mice display minimal pathogenesis, susceptible mice display severe liver pathology consistent with EVD-like disease and transcriptional signatures associated with inflammatory and liver metabolic processes. A significant quantitative trait locus (QTL) for virus RNA load in blood is identified in chromosome (chr)8, and a severe clinical disease and mortality QTL is mapped to chr7, which includes the Trim5 locus. Using knockout mice, we validate the Trim5 locus as one potential driver of liver failure and mortality after infection. The identification of susceptibility loci provides insight into molecular genetic mechanisms regulating EVD progression and severity, potentially informing therapeutics and vaccination strategies.

Properties and Mechanisms of Deletions, Insertions, and Substitutions in the Evolutionary History of SARS-CoV-2

SARS-CoV-2 has accumulated many mutations since its emergence in late 2019. Nucleotide substitutions leading to amino acid replacements constitute the primary material for natural selection. Insertions, deletions, and substitutions appear to be critical for coronavirus's macro- and microevolution. Understanding the molecular mechanisms of mutations in the mutational hotspots (positions, loci with recurrent mutations, and nucleotide context) is important for disentangling roles of mutagenesis and selection. In the SARS-CoV-2 genome, deletions and insertions are frequently associated with repetitive sequences, whereas C>U substitutions are often surrounded by nucleotides resembling the APOBEC mutable motifs. We describe various approaches to mutation spectra analyses, including the context features of RNAs that are likely to be involved in the generation of recurrent mutations. We also discuss the interplay between mutations and natural selection as a complex evolutionary trend. The substantial variability and complexity of pipelines for the reconstruction of mutations and the huge number of genomic sequences are major problems for the analyses of mutations in the SARS-CoV-2 genome. As a solution, we advocate for the development of a centralized database of predicted mutations, which needs to be updated on a regular basis.

Entropic Dynamics of Mutations in SARS-CoV-2 Genomic Sequences

In this paper, we investigate a certain class of mutations in genomic sequences by studying the evolution of the entropy and relative entropy associated with the base frequencies of a given genomic sequence. Even if the method is, in principle, applicable to every sequence which varies randomly, the case of SARS-CoV-2 RNA genome is particularly interesting to analyze, due to the richness of the available sequence database containing more than a million sequences. Our model is able to track known features of the mutation dynamics like the Cytosine-Thymine bias, but also to reveal new features of the virus mutation dynamics. We show that these new findings can be studied using an approach that combines the mean field approximation of a Markov dynamics within a stochastic thermodynamics framework.

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.

Optimal entropic properties of SARS-CoV-2 RNA sequences

The reaction of the scientific community against the COVID-19 pandemic has generated a huge (approx. 106 entries) dataset of genome sequences collected worldwide and spanning a relatively short time window. These unprecedented conditions together with the certain identification of the reference viral genome sequence allow for an original statistical study of mutations in the virus genome. In this paper, we compute the Shannon entropy of every sequence in the dataset as well as the relative entropy and the mutual information between the reference sequence and the mutated ones. These functions, originally developed in information theory, measure the information content of a sequence and allows us to study the random character of mutation mechanism in terms of its entropy and information gain or loss. We show that this approach allows us to set in new format known features of the SARS-CoV-2 mutation mechanism like the CT bias, but also to discover new optimal entropic properties of the mutation process in the sense that the virus mutation mechanism track closely theoretically computable lower bounds for the entropy decrease and the information transfer.

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.

The role of post-transcriptional regulation in SARS-CoV-2 infection and pathogenicity

The COVID-19 pandemic caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) virus has had a significant impact on global social and economic stability. To combat this, researchers have turned to omics approaches, particularly epitranscriptomics, to limit infection and develop effective therapeutic strategies. Multi-omics can provide the host response dynamics during multiple disease phases to reveal the molecular and cellular landscapes. Epitranscriptomics focuses on the mechanisms of gene transcription in cells and tissues and the relationship between genetic material and epigenetic regulation. This review highlights the role of post-transcriptional regulation in SARS-CoV-2, which affect various processes such as virus infection, replication, immunogenicity, and pathogenicity. The review also explains the formation mechanism of post-transcriptional modifications and how they can be regulated to combat viral infection and pathogenicity.

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.

Evolutionary potential of the monkeypox genome arising from interactions with human APOBEC3 enzymes

APOBEC3, an enzyme subfamily that plays a role in virus restriction by generating mutations at particular DNA motifs or mutational 'hotspots', can drive viral mutagenesis with host-specific preferential hotspot mutations contributing to pathogen variation. While previous analysis of viral genomes from the 2022 Mpox (formerly Monkeypox) disease outbreak has shown a high frequency of C>T mutations at TC motifs, suggesting recent mutations are human APOBEC3-mediated, how emerging monkeypox virus (MPXV) strains will evolve as a consequence of APOBEC3-mediated mutations remains unknown. By measuring hotspot under-representation, depletion at synonymous sites, and a combination of the two, we analyzed APOBEC3-driven evolution in human poxvirus genomes, finding varying hotspot under-representation patterns. While the native poxvirus molluscum contagiosum exhibits a signature consistent with extensive coevolution with human APOBEC3, including depletion of TC hotspots, variola virus shows an intermediate effect consistent with ongoing evolution at the time of eradication. MPXV, likely the result of recent zoonosis, showed many genes with more TC hotspots than expected by chance (over-representation) and fewer GC hotspots than expected (under-representation). These results suggest the MPXV genome: (1) may have evolved in a host with a particular APOBEC GC hotspot preference, (2) has inverted terminal repeat (ITR) regions-which may be exposed to APOBEC3 for longer during viral replication-and longer genes likely to evolve faster, and therefore (3) has a heightened potential for future human APOBEC3-meditated evolution as the virus spreads in the human population. Our predictions of MPXV mutational potential can both help guide future vaccine development and identification of putative drug targets and add urgency to the task of containing human Mpox disease transmission and uncovering the ecology of the virus in its reservoir host.

COVID-19 annual update: a narrative review

Three and a half years after the pandemic outbreak, now that WHO has formally declared that the emergency is over, COVID-19 is still a significant global issue. Here, we focus on recent developments in genetic and genomic research on COVID-19, and we give an outlook on state-of-the-art therapeutical approaches, as the pandemic is gradually transitioning to an endemic situation. The sequencing and characterization of rare alleles in different populations has made it possible to identify numerous genes that affect either susceptibility to COVID-19 or the severity of the disease. These findings provide a beginning to new avenues and pan-ethnic therapeutic approaches, as well as to potential genetic screening protocols. The causative virus, SARS-CoV-2, is still in the spotlight, but novel threatening virus could appear anywhere at any time. Therefore, continued vigilance and further research is warranted. We also note emphatically that to prevent future pandemics and other world-wide health crises, it is imperative to capitalize on what we have learnt from COVID-19: specifically, regarding its origins, the world's response, and insufficient preparedness. This requires unprecedented international collaboration and timely data sharing for the coordination of effective response and the rapid implementation of containment measures.

Increased interregional virus exchange and nucleotide diversity outline the expansion of chikungunya virus in Brazil

The emergence and reemergence of mosquito-borne diseases in Brazil such as yellow fever, zika, chikungunya, and dengue have had serious impacts on public health. Concerns have been raised due to the rapid dissemination of the chikungunya virus across the country since its first detection in 2014 in Northeast Brazil. In this work, we carried out on-site training activities in genomic surveillance in partnership with the National Network of Public Health Laboratories that have led to the generation of 422 chikungunya virus genomes from 12 Brazilian states over the past two years (2021-2022), a period that has seen more than 312 thousand chikungunya fever cases reported in the country. These genomes increased the amount of available data and allowed a more comprehensive characterization of the dispersal dynamics of the chikungunya virus East-Central-South-African lineage in Brazil. Tree branching patterns revealed the emergence and expansion of two distinct subclades. Phylogeographic analysis indicated that the northeast region has been the leading hub of virus spread towards other regions. Increased frequency of C > T transitions among the new genomes suggested that host restriction factors from the immune system such as ADAR and AID/APOBEC deaminases might be driving the genetic diversity of the chikungunya virus in Brazil.

Evolutionary potential of the monkeypox genome arising from interactions with human APOBEC3 enzymes

APOBEC3, an enzyme subfamily that plays a role in virus restriction by generating mutations at particular DNA motifs or mutational "hotspots," can drive viral mutagenesis with host-specific preferential hotspot mutations contributing to pathogen variation. While previous analysis of viral genomes from the 2022 Mpox (formerly Monkeypox) disease outbreak has shown a high frequency of C>T mutations at T C motifs, suggesting recent mutations are human APOBEC3-mediated, how emerging monkeypox virus (MPXV) strains will evolve as a consequence of APOBEC3-mediated mutations remains unknown. By measuring hotspot under-representation, depletion at synonymous sites, and a combination of the two, we analyzed APOBEC3-driven evolution in human poxvirus genomes, finding varying hotspot under-representation patterns. While the native poxvirus molluscum contagiosum exhibits a signature consistent with extensive coevolution with human APOBEC3, including depletion of T C hotspots, variola virus shows an intermediate effect consistent with ongoing evolution at the time of eradication. MPXV, likely the result of recent zoonosis, showed many genes with more T C hotspots than expected by chance (over-representation) and fewer G C hotspots than expected (under-representation). These results suggest the MPXV genome: 1) may have evolved in a host with a particular APOBEC G C hotspot preference, 2) has inverted terminal repeat (ITR) regions -which may be exposed to APOBEC3 for longer during viral replication- and longer genes likely to evolve faster, and therefore 3) has a heightened potential for future human APOBEC3-meditated evolution as the virus spreads in the human population. Our predictions of MPXV mutational potential can both help guide future vaccine development and identification of putative drug targets and add urgency to the task of containing human Mpox disease transmission and uncovering the ecology of the virus in its reservoir host.

The Mutational Landscape of SARS-CoV-2

Mutation research is crucial for detecting and treating SARS-CoV-2 and developing vaccines. Using over 5,300,000 sequences from SARS-CoV-2 genomes and custom Python programs, we analyzed the mutational landscape of SARS-CoV-2. Although almost every nucleotide in the SARS-CoV-2 genome has mutated at some time, the substantial differences in the frequency and regularity of mutations warrant further examination. C>U mutations are the most common. They are found in the largest number of variants, pangolin lineages, and countries, which indicates that they are a driving force behind the evolution of SARS-CoV-2. Not all SARS-CoV-2 genes have mutated in the same way. Fewer non-synonymous single nucleotide variations are found in genes that encode proteins with a critical role in virus replication than in genes with ancillary roles. Some genes, such as spike (S) and nucleocapsid (N), show more non-synonymous mutations than others. Although the prevalence of mutations in the target regions of COVID-19 diagnostic RT-qPCR tests is generally low, in some cases, such as for some primers that bind to the N gene, it is significant. Therefore, ongoing monitoring of SARS-CoV-2 mutations is crucial. The SARS-CoV-2 Mutation Portal provides access to a database of SARS-CoV-2 mutations.

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.

The evolution of SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused millions of deaths and substantial morbidity worldwide. Intense scientific effort to understand the biology of SARS-CoV-2 has resulted in daunting numbers of genomic sequences. We witnessed evolutionary events that could mostly be inferred indirectly before, such as the emergence of variants with distinct phenotypes, for example transmissibility, severity and immune evasion. This Review explores the mechanisms that generate genetic variation in SARS-CoV-2, underlying the within-host and population-level processes that underpin these events. We examine the selective forces that likely drove the evolution of higher transmissibility and, in some cases, higher severity during the first year of the pandemic and the role of antigenic evolution during the second and third years, together with the implications of immune escape and reinfections, and the increasing evidence for and potential relevance of recombination. In order to understand how major lineages, such as variants of concern (VOCs), are generated, we contrast the evidence for the chronic infection model underlying the emergence of VOCs with the possibility of an animal reservoir playing a role in SARS-CoV-2 evolution, and conclude that the former is more likely. We evaluate uncertainties and outline scenarios for the possible future evolutionary trajectories of SARS-CoV-2.

Increased interregional virus exchange and nucleotide diversity outline the expansion of the chikungunya virus ECSA lineage in Brazil

The emergence and reemergence of mosquito-borne diseases in Brazil such as Yellow Fever, Zika, Chikungunya, and Dengue have had serious impacts on public health. Concerns have been raised due to the rapid dissemination of the chikungunya virus (CHIKV) across the country since its first detection in 2014 in Northeast Brazil. Faced with this scenario, on-site training activities in genomic surveillance carried out in partnership with the National Network of Public Health Laboratories have led to the generation of 422 CHIKV genomes from 12 Brazilian states over the past two years (2021-2022), a period that has seen more than 312 thousand chikungunya fever cases reported in the country. These new genomes increased the amount of available data and allowed a more comprehensive characterization of the dispersion dynamics of the CHIKV East-Central-South-African (ECSA) lineage in Brazil. Tree branching patterns revealed the emergence and expansion of two distinct subclades. Phylogeographic analysis indicated that the northeast region has been the leading hub of virus spread towards other regions. Increased frequency of C>T transitions among the new genomes suggested that host restriction factors from the immune system such as ADAR and AID/APOBEC deaminases might be driving CHIKV ECSA lineage genetic diversity in Brazil.

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.

C-to-U RNA deamination is the driving force accelerating SARS-CoV-2 evolution

Understanding the molecular mechanism underlying the rampant mutation of SARS-CoV-2 would help us control the COVID-19 pandemic. The APOBEC-mediated C-to-U deamination is a major mutation type in the SARS-CoV-2 genome. However, it is unclear whether the novel mutation rate u is higher for C-to-U than for other mutation types, and what the detailed driving force is. By analyzing the time course SARS-CoV-2 global population data, we found that C-to-U has the highest novel mutation rate u among all mutation types and that this u is still increasing with time (du/dt > 0). Novel C-to-U events, rather than other mutation types, have a preference over particular genomic regions. A less local RNA structure is correlated with a high novel C-to-U mutation rate. A cascade model nicely explains the du/dt > 0 for C-to-U deamination. In SARS-CoV-2, the RNA structure serves as the molecular basis of the extremely high and continuously accelerating C-to-U deamination rate. This mechanism is the driving force of the mutation, adaptation, and evolution of SARS-CoV-2. Our findings help us understand the dynamic evolution of the virus mutation rate.

Examples of Structural Motifs in Viral Genomes and Approaches for RNA Structure Characterization

The relationship between conserved structural motifs and their biological function in the virus replication cycle is the interest of many researchers around the world. RNA structure is closely related to RNA function. Therefore, technological progress in high-throughput approaches for RNA structure analysis and the development of new ones are very important. In this mini review, we discuss a few perspectives on the structural elements of viral genomes and some methods used for RNA structure prediction and characterization. Based on the recent literature, we describe several examples of studies concerning the viral genomes, especially severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and influenza A virus (IAV). Herein, we emphasize that a better understanding of viral genome architecture allows for the discovery of the structure-function relationship, and as a result, the discovery of new potential antiviral therapeutics.

Prediction of Recurrent Mutations in SARS-CoV-2 Using Artificial Neural Networks

Predicting SARS-CoV-2 mutations is difficult, but predicting recurrent mutations driven by the host, such as those caused by host deaminases, is feasible. We used machine learning to predict which positions from the SARS-CoV-2 genome will hold a recurrent mutation and which mutations will be the most recurrent. We used data from April 2021 that we separated into three sets: a training set, a validation set, and an independent test set. For the test set, we obtained a specificity value of 0.69, a sensitivity value of 0.79, and an Area Under the Curve (AUC) of 0.8, showing that the prediction of recurrent SARS-CoV-2 mutations is feasible. Subsequently, we compared our predictions with updated data from January 2022, showing that some of the false positives in our prediction model become true positives later on. The most important variables detected by the model's Shapley Additive exPlanation (SHAP) are the nucleotide that mutates and RNA reactivity. This is consistent with the SARS-CoV-2 mutational bias pattern and the preference of some host deaminases for specific sequences and RNA secondary structures. We extend our investigation by analyzing the mutations from the variants of concern Alpha, Beta, Delta, Gamma, and Omicron. Finally, we analyzed amino acid changes by looking at the predicted recurrent mutations in the M-pro and spike proteins.

The molecular epidemiology of multiple zoonotic origins of SARS-CoV-2

Understanding the circumstances that lead to pandemics is important for their prevention. We analyzed the genomic diversity of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) early in the coronavirus disease 2019 (COVID-19) pandemic. We show that SARS-CoV-2 genomic diversity before February 2020 likely comprised only two distinct viral lineages, denoted "A" and "B." Phylodynamic rooting methods, coupled with epidemic simulations, reveal that these lineages were the result of at least two separate cross-species transmission events into humans. The first zoonotic transmission likely involved lineage B viruses around 18 November 2019 (23 October to 8 December), and the separate introduction of lineage A likely occurred within weeks of this event. These findings indicate that it is unlikely that SARS-CoV-2 circulated widely in humans before November 2019 and define the narrow window between when SARS-CoV-2 first jumped into humans and when the first cases of COVID-19 were reported. As with other coronaviruses, SARS-CoV-2 emergence likely resulted from multiple zoonotic events.

The complexity landscape of viral genomes

BACKGROUND

Viruses are among the shortest yet highly abundant species that harbor minimal instructions to infect cells, adapt, multiply, and exist. However, with the current substantial availability of viral genome sequences, the scientific repertory lacks a complexity landscape that automatically enlights viral genomes' organization, relation, and fundamental characteristics.

RESULTS

This work provides a comprehensive landscape of the viral genome's complexity (or quantity of information), identifying the most redundant and complex groups regarding their genome sequence while providing their distribution and characteristics at a large and local scale. Moreover, we identify and quantify inverted repeats abundance in viral genomes. For this purpose, we measure the sequence complexity of each available viral genome using data compression, demonstrating that adequate data compressors can efficiently quantify the complexity of viral genome sequences, including subsequences better represented by algorithmic sources (e.g., inverted repeats). Using a state-of-the-art genomic compressor on an extensive viral genomes database, we show that double-stranded DNA viruses are, on average, the most redundant viruses while single-stranded DNA viruses are the least. Contrarily, double-stranded RNA viruses show a lower redundancy relative to single-stranded RNA. Furthermore, we extend the ability of data compressors to quantify local complexity (or information content) in viral genomes using complexity profiles, unprecedently providing a direct complexity analysis of human herpesviruses. We also conceive a features-based classification methodology that can accurately distinguish viral genomes at different taxonomic levels without direct comparisons between sequences. This methodology combines data compression with simple measures such as GC-content percentage and sequence length, followed by machine learning classifiers.

CONCLUSIONS

This article presents methodologies and findings that are highly relevant for understanding the patterns of similarity and singularity between viral groups, opening new frontiers for studying viral genomes' organization while depicting the complexity trends and classification components of these genomes at different taxonomic levels. The whole study is supported by an extensive website (https://asilab.github.io/canvas/) for comprehending the viral genome characterization using dynamic and interactive approaches.

Rampant C-to-U deamination accounts for the intrinsically high mutation rate in SARS-CoV-2 spike gene

The high mutation rate of SARS-CoV-2 largely complicates our control of the pandemic. In particular, it is currently unclear why the spike (S) gene has an extraordinarily high mutation rate among all SARS-CoV-2 genes. By analyzing the occurrence of fixed synonymous mutations between SARS-CoV-2 and RaTG13, and profiling the DAF (derived allele frequency) of polymorphic synonymous sites among millions of worldwide SARS-CoV-2 strains, we found that both fixed and polymorphic mutations show higher mutation rates in the S gene than other genes. The majority of mutations are C-to-T, representing the APOBEC-mediated C-to-U deamination instead of the previously proposed A-to-I deamination. Both in silico and in vivo evidence indicated that the S gene is more likely to be single-stranded compared to other SARS-CoV-2 genes, agreeing with the APOBEC preference of ssRNA. We conclude that the single-stranded property of the S gene makes it a favorable target for C-to-U deamination, leading to its excessively high mutation rate compared to other non-S genes. In conclusion, APOBEC, rather than ADAR, is the "editor-in-chief" of SARS-CoV-2 RNAs. This work helps us to understand the molecular mechanism underlying the mutation and evolution of SARS-CoV-2, and we believe it will contribute to the control of the pandemic.

Variation in synonymous nucleotide composition among genomes of sarbecoviruses and consequences for the origin of COVID-19

The subgenus Sarbecovirus includes two human viruses, SARS-CoV and SARS-CoV-2, respectively responsible for the SARS epidemic and COVID-19 pandemic, as well as many bat viruses and two pangolin viruses. Here, the synonymous nucleotide composition (SNC) of Sarbecovirus genomes was analysed by examining third codon-positions, dinucleotides, and degenerate codons. The results show evidence for the eight following groups: (i) SARS-CoV related coronaviruses (SCoVrC including many bat viruses from China), (ii) SARS-CoV-2 related coronaviruses (SCoV2rC; including five bat viruses from Cambodia, Thailand and Yunnan), (iii) pangolin sarbecoviruses, (iv) three bat sarbecoviruses showing evidence of recombination between SCoVrC and SCoV2rC genomes, (v) two highly divergent bat sarbecoviruses from Yunnan, (vi) the bat sarbecovirus from Japan, (vii) the bat sarbecovirus from Bulgaria, and (viii) the bat sarbecovirus from Kenya. All these groups can be diagnosed by specific nucleotide compositional features except the one concerned by recombination between SCoVrC and SCoV2rC. In particular, SCoV2rC genomes have less cytosines and more uracils at third codon-positions than other sarbecoviruses, whereas the genomes of pangolin sarbecoviruses show more adenines at third codon-positions. I suggest that taxonomic differences in the imbalanced nucleotide pools available in host cells during viral replication can explain the eight groups of SNC here detected among Sarbecovirus genomes. A related effect due to hibernating bats and their latitudinal distribution is also discussed. I conclude that the two independent host switches from Rhinolophus bats to pangolins resulted in convergent mutational constraints and that SARS-CoV-2 emerged directly from a horseshoe bat sarbecovirus.

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.

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.

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.