Potential APOBEC-mediated RNA editing of the genomes of SARS-CoV-2 and other coronaviruses and its impact on their longer term evolution

Enrichment of G-to-U Substitution in SARS-CoV-2 Functional Regions and Its Compensation via Concurrent Mutations

We surveyed single nucleotide variant (SNV) patterns from 5 903 647 complete SARS-CoV-2 genomes. Among 10 012 SNVs, APOBEC-mediated C-to-U (C > U) deamination was the most prevalent, followed by G > U and other RNA editing-related substitutions including (A > G, U > C, G > A). However, C > U mutations were less frequent in functional regions, for example, S protein, intrinsic disordered regions, and nonsynonymous mutations, where G > U were over-represented. Notably, G-loss substitutions rarely appeared together. Instead, G-gain mutations tended to more frequently co-occur with others, with a marked preference in the S protein, suggesting a compensatory mechanism for G loss in G > U mutations. The temporal patterns revealed C > U frequency declined until late 2021 then resurged in early 2022. Conversely, G > U steadily decreased, with a pronounced drop in January 2022, coinciding with reduced COVID-19 severity. Vaccinated individuals exhibited a slightly but significantly higher C > U frequency and a notably lower G > U frequency compared to the unvaccinated group. Additionally, cancer patients had higher G > U frequency than general patients during the same period. Interestingly, none of the C > U SNVs were uniquely identified in 2724 environmental samples. These findings suggest novel functional roles of G > U in COVID-19 symptoms, potentially linked to oxidative stress and reactive oxygen species, while C > U remains the dominant substitution, likely driven by host immune-mediated RNA editing.

The differential effect of SARS-CoV-2 NSP1 on mRNA translation and stability reveals new insights linking ribosome recruitment, codon usage, and virus evolution

The nonstructural protein 1 (NSP1) of SARS-CoV-2 blocks the messenger RNA (mRNA) entry channel of the 40S ribosomal subunit, causing inhibition of translation initiation and subsequent degradation of host mRNAs. However, target mRNA specificity and how viral mRNAs escape NSP1-mediated degradation have not been clarified to date. Here we found that NSP1 acts as a translational switch capable of blocking or enhancing translation depending on how preinitiation complex, 43S-PIC, is recruited to the mRNA, whereas NSP1-mediated mRNA degradation mostly depends on codon usage bias. Thus, fast-translating mRNAs with optimal codon usage for human cells that preferentially recruit 43S-PIC by threading showed a dramatic sensitivity to NSP1. Translation of SARS-CoV-2 mRNAs escapes NSP1-mediated inhibition by a proper combination of suboptimal codon usage and slotting-prone 5' UTR. Thus, the prevalence of nonoptimal codons found in SARS-CoV-2 and other coronavirus genomes is favored by the distinctive effect that NSP1 plays on translation and mRNA stability.

Intra-host variability of SARS-CoV-2: Patterns, causes and impact on COVID-19

Intra-host viral variability is related to pathogenicity, persistence, drug resistance, and the emergence of new clades. This work reviews the large amount of data on SARS-CoV-2 intra-host variability accumulated to date, addressing known and potential implications in COVID-19 and the emergence of VOCs and lineage-defining mutations. Topics covered include the distribution of intra-host polymorphisms across the genome, the corresponding mutational signatures, their patterns of emergence and extinction throughout infection, and the processes governing their abundance, frequency, and type (synonymous, nonsynonymous, indels, nonsense). Besides, evidence is reviewed that the virus can replicate and mutate in isolation at different anatomical compartments, which may imply that what we have learned from respiratory samples could be part of a broader picture.

C→U transition biases in SARS-CoV-2: still rampant 4 years from the start of the COVID-19 pandemic

The evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the pandemic and post-pandemic periods has been characterized by rapid adaptive changes that confer immune escape and enhanced human-to-human transmissibility. Sequence change is additionally marked by an excess number of C→U transitions suggested as being due to host-mediated genome editing. To investigate how these influence the evolutionary trajectory of SARS-CoV-2, 2,000 high-quality, coding complete genome sequences of SARS-CoV-2 variants collected pre-September 2020 and from each subsequently appearing alpha, delta, BA.1, BA.2, BA.5, XBB, EG, HK, and JN.1 lineages were downloaded from NCBI Virus in April 2024. C→U transitions were the most common substitution during the diversification of SARS-CoV-2 lineages over the 4-year observation period. A net loss of C bases and accumulation of U's occurred at a constant rate of approximately 0.2%-0.25%/decade. C→U transitions occurred in over a quarter of all sites with a C (26.5%; range 20.0%-37.2%) around five times more than observed for the other transitions (5.3%-6.8%). In contrast to an approximately random distribution of other transitions across the genome, most C→U substitutions occurred at statistically preferred sites in each lineage. However, only the most C→U polymorphic sites showed evidence for a preferred 5'U context previously associated with APOBEC 3A editing. There was a similarly weak preference for unpaired bases suggesting much less stringent targeting of RNA than mediated by A3 deaminases in DNA editing. Future functional studies are required to determine editing preferences, impacts on replication fitness in vivo of SARS-CoV-2 and other RNA viruses, and impact on host tropism.

IMPORTANCE

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the pandemic and post-pandemic periods has shown a remarkable capacity to adapt and evade human immune responses and increase its human-to-human transmissibility. The genome of SARS-CoV-2 is also increasingly scarred by the effects of multiple C→U mutations from host genome editing as a cellular defense mechanism akin to restriction factors for retroviruses. Through the analysis of large data sets of SARS-CoV-2 isolate sequences collected throughout the pandemic period and beyond, we show that C→U transitions have driven a base compositional change over time amounting to a net loss of C bases and accumulation of U's at a rate of approximately 0.2%-0.25%/decade. Most C→U substitutions occurred in the absence of the preferred upstream-base context or targeting of unpaired RNA bases previously associated with the host RNA editing protein, APOBEC 3A. The analyses provide a series of testable hypotheses that can be experimentally investigated in the future.

MDA5 Is a Major Determinant of Developing Symptoms in Critically Ill COVID-19 Patients

Apart from the skin and mucosal immune barrier, the first line of defense of the human immune system includes MDA5 (ifih1 gene) which acts as a cellular sensor protein for certain viruses including SARS-CoV-2. Upon binding with viral RNA, MDA5 activates cell-intrinsic innate immunity, humoral responses, and MAVS (mitochondrial antiviral signaling). MAVS signaling induces type I and III interferon (IFN) expressions that further induce ISGs (interferon stimulatory genes) expressions to initiate human cell-mediated immune responses and attenuate viral replication. SARS-CoV-2 counteracts by producing NSP1, NSP2, NSP3, NSP5, NSP7, NSP12, ORF3A, ORF9, N, and M protein and directs anti-MDA5 antibody production presumably to antagonize IFN signaling. Furthermore, COVID-19 resembles several diseases that carry anti-MDA5 antibodies and the current COVID-19 vaccines induced anti-MDA5 phenotypes in healthy individuals. GWAS (genome-wide association studies) identified several polymorphisms (SNPs) in the ifih1-ifn pathway genes including rs1990760 in ifih1 that are strongly associated with COVID-19, and the associated risk allele is correlated with reduced IFN production. The genetic association of SNPs in ifih1 and ifih1-ifn pathway genes reinforces the molecular findings of the critical roles of MDA5 in sensing SARS-CoV-2 and subsequently the IFN responses to inhibit viral replication and host immune evasion. Thus, MDA5 or its pathway genes could be targeted for therapeutic development of COVID-19.

Analysis of the abundance and diversity of RNA secondary structure elements in RNA viruses using the RNAsselem Python package

Recent advancements in experimental and computational methods for RNA secondary structure detection have revealed the crucial role of RNA structural elements in diverse molecular processes within living cells. It has been demonstrated that the secondary structure of the entire viral genome is often responsible for performing crucial functions in the viral life cycle and also influences virus evolution. To investigate the role of viral RNA secondary structure, alongside experimental techniques, the use of bioinformatics tools is important for analyzing various secondary structure patterns, including hairpin loops, internal loops, multifurcations, external loops, bulges, stems, and pseudoknots. Here, we have introduced a Python package for analyzing RNA secondary structure elements in viral genomes, which includes the recognition of common secondary structure patterns, the generation of descriptive statistics for these structural elements, and the provision of their basic properties. We applied the developed package to analyze the secondary structures of complete viral genomes collected from the literature, aiming to gain insights into viral function and evolution. Both the package and the collection of secondary structures of viral genomes are available at http://github.com/KazanovLab/RNAsselem .

Apolipoprotein-B mRNA-editing complex 3B could be a new potential therapeutic target in endometriosis

This study investigated the correlation of Apolipoprotein-B mRNA-editing complex 3B (APOBEC3B) expression with hypoxia inducible factor 1α (HIF-1α), Kirsten rat sarcoma virus (KRAS) and phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) in endometriosis patients, and the inhibitory effects of APOBEC3B knockdown in a human endometriotic cell line. Here, APOBEC3B, HIF-1α, KRAS, and PIK3CA were examined in patients with and without endometriosis using reverse transcription polymerase chain reaction (RT-PCR). The apoptosis, cell proliferation, invasion, migration, and biological function of APOBEC3B knockdown were explored in 12Z immortalized human endometriotic cell line. We observed APOBEC3B, HIF-1α, KRAS and PIK3CA expressions were significantly higher in endometriosis patients (p < 0.001, p < 0.001, p = 0.029, p = 0.001). Knockdown of APOBEC3B increased apoptosis, which was 28.03% and 22.27% higher than in mock and control siRNA samples, respectively. APOBEC3B knockdown also decreased PIK3CA expression and increased Caspase 8 expression, suggesting a potential role in the regulation of apoptosis. Furthermore, knockdown of APOBEC3B significantly inhibited cell proliferation, invasion, and migration compared to mock and control siRNA. (Cell proliferation: mock: p < 0.001 and control siRNA: p = 0.049. Cell invasion: mock: p < 0.001 and control siRNA: p = 0.029. Cell migration: mock: p = 0.004, and control siRNA: p = 0.014). In conclusion, this study suggests that APOBEC3B may be a new potential therapeutic target for endometriosis.

The Emergence and Evolution of SARS-CoV-2

The origin of SARS-CoV-2 has evoked heated debate and strong accusations, yet seemingly little resolution. I review the scientific evidence on the origin of SARS-CoV-2 and its subsequent spread through the human population. The available data clearly point to a natural zoonotic emergence within, or closely linked to, the Huanan Seafood Wholesale Market in Wuhan. There is no direct evidence linking the emergence of SARS-CoV-2 to laboratory work conducted at the Wuhan Institute of Virology. The subsequent global spread of SARS-CoV-2 was characterized by a gradual adaptation to humans, with dual increases in transmissibility and virulence until the emergence of the Omicron variant. Of note has been the frequent transmission of SARS-CoV-2 from humans to other animals, marking it as a strongly host generalist virus. Unless lessons from the origin of SARS-CoV-2 are learned, it is inevitable that more zoonotic events leading to more epidemics and pandemics will plague human populations.

The nuts and bolts of recombination in the generation of SARS-CoV-2 variants; from XA to XBB

Since the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), new variants with enhanced transmissibility and pathogenicity have surfaced. The World Health Organization has designated five such variants-Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2), and Omicron (B.1.1.529)-as variants of concern. Each variant exhibits distinct characteristics, with many displaying a combination of point mutations and insertions/deletions (indels). These genetic alterations, including mutations, recombinations, and rearrangements, contribute to the emergence of new strains that may exhibit modified phenotypes. However, identifying recombinant forms can be challenging due to their resemblance to other lineages. It is critical to monitor the evolution of new recombinant variants, particularly in light of the potential for vaccine-resistant strains and their accelerated propagation. Recombination has played a pivotal role in the development of certain SARS-CoV-2 variants, such as XA, XD, XF, XE, and XBB, among others. This report delves into the significance of recombination in the evolution of SARS-CoV-2 variants, especially Omicron sublineages, underscoring the necessity for continuous surveillance of the SARS-CoV-2 genome to identify newly emerged recombinant variants.

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.

RNA editing in host lncRNAs as potential modulator in SARS-CoV-2 variants-host immune response dynamics

Both host and viral RNA editing plays a crucial role in host's response to infection, yet our understanding of host RNA editing remains limited. In this study of in-house generated RNA sequencing (RNA-seq) data of 211 hospitalized COVID-19 patients with PreVOC, Delta, and Omicron variants, we observed a significant differential editing frequency and patterns in long non-coding RNAs (lncRNAs), with Delta group displaying lower RNA editing compared to PreVOC/Omicron patients. Notably, multiple transcripts of UGDH-AS1 and NEAT1 exhibited high editing frequencies. Expression of ADAR1/APOBEC3A/APOBEC3G and differential abundance of repeats were possible modulators of differential editing across patient groups. We observed a shift in crucial infection-related pathways wherein the pathways were downregulated in Delta compared to PreVOC and Omicron. Our genomics-based evidence suggests that lncRNA editing influences stability, miRNA binding, and expression of both lncRNA and target genes. Overall, the study highlights the role of lncRNAs and how editing within host lncRNAs modulates the disease severity.

Vaccination Shapes Within-Host SARS-CoV-2 Diversity of Omicron BA.2.2 Breakthrough Infection

BACKGROUND

Low-frequency intrahost single-nucleotide variants of SARS-CoV-2 have been recognized as predictive indicators of selection. However, the impact of vaccination on the intrahost evolution of SARS-CoV-2 remains uncertain at present.

METHODS

We investigated the genetic variation of SARS-CoV-2 in individuals who were unvaccinated, partially vaccinated, or fully vaccinated during Shanghai's Omicron BA.2.2 wave. We substantiated the connection between particular amino acid substitutions and immune-mediated selection through a pseudovirus neutralization assay or by cross-verification with the human leukocyte antigen-associated T-cell epitopes.

RESULTS

In contrast to those with immunologic naivety or partial vaccination, participants who were fully vaccinated had intrahost variant spectra characterized by reduced diversity. Nevertheless, the distribution of mutations in the fully vaccinated group was enriched in the spike protein. The distribution of intrahost single-nucleotide variants in individuals who were immunocompetent did not demonstrate notable signs of positive selection, in contrast to the observed adaptation in 2 participants who were immunocompromised who had an extended period of viral shedding.

CONCLUSIONS

In SARS-CoV-2 infections, vaccine-induced immunity was associated with decreased diversity of within-host variant spectra, with milder inflammatory pathophysiology. The enrichment of mutations in the spike protein gene indicates selection pressure exerted by vaccination on the evolution of SARS-CoV-2.

Most accurate mutations in SARS-CoV-2 genomes identified in Uzbek patients show novel amino acid changes

PURPOSE

The rapid changes in the coronavirus genomes created new strains after the first variation was found in Wuhan in 2019. SARS-CoV-2 genotypes should periodically undergo whole genome sequencing to control it because it has been extremely helpful in combating the virus. Many diagnoses, treatments, and vaccinations have been developed against it based on genome sequencing. With its practical implications, this study aimed to determine changes in the delta variant of SARS-CoV-2 widespread in Uzbekistan during the pandemic by genome sequencing, thereby providing crucial insights for developing effective control strategies that can be directly applied in the field.

DESIGN

We meticulously generated 17 high-quality whole-genome sequence data from 48 SARS-CoV-2 genotypes of COVID-19 patients who tested positive by PCR in Tashkent, Uzbekistan. Our rigorous approach, which includes stringent quality control measures and multiple rounds of verification, ensures the accuracy and reliability of our findings.

METHODS

Our study employed a unique combination of genome sequencing and bioinformatics web tools to analyze amino acid (AA) changes in the virus genomes. This approach allowed us to understand the genetic changes in the delta variant of SARS-CoV-2 widespread in Uzbekistan during the pandemic.

RESULTS

Our study revealed significant nucleotide polymorphisms, including non-synonymous (missense) and synonymous mutations in the coding regions of the sequenced sample genomes. These findings, categorized by phylogenetic analysis into the G clade (or GK sub-clade), contribute to our understanding of the delta variant of SARS-CoV-2 widespread in Uzbekistan during the pandemic. A total of 134 mutations were identified, consisting of 65 shared and 69 unique mutations. These nucleotide changes, including one frameshift mutation, one conservative and disruptive insertion-deletion, four upstream region mutations, four downstream region mutations, 39 synonymous mutations, and 84 missense mutations, are crucial in the ongoing battle against the virus.

CONCLUSION

The comprehensive whole-genome sequencing data presented in this study aids in tracing the origins and sources of circulating SARS-CoV-2 variants and analyzing emerging variations within Uzbekistan and globally. The genome sequencing of SARS-CoV-2 from samples collected in Uzbekistan in late 2021, during the peak of the pandemic's second wave nationwide, is detailed here. Following acquiring these sequences, research efforts have focused on developing DNA and plant-based edible vaccines utilizing prevalent SARS-CoV-2 strains in Uzbekistan, which are currently undergoing clinical trials.

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.

The RdRp genotyping of SARS-CoV-2 isolated from patients with different clinical spectrum of COVID-19

BACKGROUND

The evolution of SARS-CoV-2 has been observed from the very beginning of the fight against COVID-19, some mutations are indicators of potentially dangerous variants of the virus. However, there is no clear association between the genetic variants of SARS-CoV-2 and the severity of COVID-19. We aimed to analyze the genetic variability of RdRp in correlation with different courses of COVID-19.

RESULTS

The prospective study included 77 samples of SARS-CoV-2 isolated from outpatients (1st degree of severity) and hospitalized patients (2nd, 3rd and 4th degree of severity). The retrospective analyses included 15,898,266 cases of SARS-CoV-2 genome sequences deposited in the GISAID repository. Single-nucleotide variants were identified based on the four sequenced amplified fragments of SARS-CoV-2. The analysis of the results was performed using appropriate statistical methods, with p < 0.05, considered statistically significant. Additionally, logistic regression analysis was performed to predict the strongest determinants of the observed relationships. The number of mutations was positively correlated with the severity of the COVID-19, and older male patients. We detected four mutations that significantly increased the risk of hospitalization of COVID-19 patients (14676C > T, 14697C > T, 15096 T > C, and 15279C > T), while the 15240C > T mutation was common among strains isolated from outpatients. The selected mutations were searched worldwide in the GISAID database, their presence was correlated with the severity of COVID-19.

CONCLUSION

Identified mutations have the potential to be used to assess the increased risk of hospitalization in COVID-19 positive patients. Experimental studies and extensive epidemiological data are needed to investigate the association between individual mutations and the severity of COVID-19.

A simple method for estimating time-irreversible nucleotide substitution rates in the SARS-CoV-2 genome

SARS-CoV-2 is the cause of the current worldwide pandemic of severe acute respiratory syndrome. The change of nucleotide composition of the SARS-CoV-2 genome is crucial for understanding the spread and transmission dynamics of the virus because viral nucleotide sequences are essential in identifying viral strains. Recent studies have shown that cytosine (C) to uracil (U) substitutions are overrepresented in SARS-CoV-2 genome sequences. These asymmetric substitutions between C and U indicate that traditional time-reversible substitution models cannot be applied to the evolution of SARS-CoV-2 sequences. Thus, we develop a new time-irreversible model of nucleotide substitutions to estimate the substitution rates in SARS-CoV-2 genomes. We investigated the number of nucleotide substitutions among the 7862 genomic sequences of SARS-CoV-2 registered in the Global Initiative on Sharing All Influenza Data (GISAID) that have been sampled from all over the world. Using the new method, the substitution rates in SARS-CoV-2 genomes were estimated. The C-to-U substitution rates of SARS-CoV-2 were estimated to be 1.95 × 10-3 ± 4.88 × 10-4 per site per year, compared with 1.48 × 10-4 ± 7.42 × 10-5 per site per year for all other types of substitutions.

The Intricate Interplay between APOBEC3 Proteins and DNA Tumour Viruses

APOBEC3 proteins are cytidine deaminases that play a crucial role in the innate immune response against viruses, including DNA viruses. Their main mechanism for restricting viral replication is the deamination of cytosine to uracil in viral DNA during replication. This process leads to hypermutation of the viral genome, resulting in loss of viral fitness and, in many cases, inactivation of the virus. APOBEC3 proteins inhibit the replication of a number of DNA tumour viruses, including herpesviruses, papillomaviruses and hepadnaviruses. Different APOBEC3s restrict the replication of different virus families in different ways and this restriction is not limited to one APOBEC3. Infection with DNA viruses often leads to the development and progression of cancer. APOBEC3 mutational signatures have been detected in various cancers, indicating the importance of APOBEC3s in carcinogenesis. Inhibition of DNA viruses by APOBEC3 proteins appears to play a dual role in this process. On the one hand, it is an essential component of the innate immune response to viral infections, and, on the other hand, it contributes to the pathogenesis of persistent viral infections and the progression of cancer. The current review examines the complex interplay between APOBEC3 proteins and DNA viruses and sheds light on the mechanisms of action, viral countermeasures and the impact on carcinogenesis. Deciphering the current issues in the interaction of APOBEC/DNA viruses should enable the development of new targeted cancer therapies.

Resistance-associated mutations to the anti-SARS-CoV-2 agent nirmatrelvir: Selection not induction

Mutations associated with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) resistance to antiprotease nirmatrelvir were reported. We aimed to detect them in SARS-CoV-2 genomes and quasispecies retrieved in our institute before drug availability in January 2022 and to analyze the impact of mutations on protease (3CLpro) structure. We sought for 38 3CLpro nirmatrelvir resistance mutations in a set of 62 673 SARS-CoV-2 genomes obtained in our institute from respiratory samples collected between 2020 and 2023 and for these mutations in SARS-CoV-2 quasispecies for 90 samples collected in 2020, using Python. SARS-CoV-2 protease with major mutation E166V was generated with Swiss Pdb Viewer and Molegro Molecular Viewer. We detected 22 (58%) of the resistance-associated mutations in 417 (0.67%) of the genomes analyzed; 325 (78%) of these genomes had been obtained from samples collected in 2020-2021. APOBEC signatures were found for 12/22 mutations. We also detected among viral quasispecies from 90 samples some minority reads harboring any of 15 nirmatrelvir resistance mutations, including E166V. Also, we predicted that E166V has a very limited effect on 3CLpro structure but may prevent drug attachment. Thus, we evidenced that mutations associated with nirmatrelvir resistance pre-existed in SARS-CoV-2 before drug availability. These findings further warrant SARS-CoV-2 genomic surveillance and SARS-CoV-2 quasispecies characterization.

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.

Anelloviruses versus human immunity: how do we control these viruses?

One continuous companion and one of the major players in the human blood virome are members of the Anelloviridae family. Anelloviruses are probably found in all humans, infection occurs early in life and the composition (anellome) is thought to remain stable and personal during adulthood. The stable anellome implies a great balance between the host immune system and the virus. However, the lack of a robust culturing system hampers direct investigation of interactions between virus and host cells. Other techniques, however, including next generation sequencing, AnelloScan-antibody tests, evolution selection pressure analysis, and virus protein structures, do provide new insights into the interactions between anelloviruses and the host immune system. This review aims at providing an overview of the current knowledge on the immune mechanisms acting on anelloviruses and the countering viral mechanisms allowing immune evasion.

The Impact of First-Time SARS-CoV-2 Infection on Human Anelloviruses

Members of the Anelloviridae family dominate the blood virome, emerging early in life. The anellome, representing the variety of anelloviruses within an individual, stabilizes by adulthood. Despite their supposedly commensal nature, elevated anellovirus concentrations under immunosuppressive treatment indicate an equilibrium controlled by immunity. Here, we investigated whether anelloviruses are sensitive to the immune activation that accompanies a secondary infection. As a model, we investigated 19 health care workers (HCWs) with initial SARS-CoV-2 infection, with blood sampling performed pre and post infection every 4 weeks in a 3-month-follow-up during the early 2020 COVID-19 pandemic. A concurrently followed control group (n = 27) remained SARS-CoV-2-negative. Serum anellovirus loads were measured using qPCR. A significant decrease in anellovirus load was found in the first weeks after SARS-CoV-2 infection, whereas anellovirus concentrations remained stable in the uninfected control group. A restored anellovirus load was seen approximately 10 weeks after SARS-CoV-2 infection. For five subjects, an in-time anellome analysis via Illumina sequencing could be performed. In three of the five HCWs, the anellome visibly changed during SARS-CoV-2 infection and returned to baseline in two of these cases. In conclusion, anellovirus loads in blood can temporarily decrease upon an acute secondary infection.

Is it time to switch to a formulation other than the live attenuated poliovirus vaccine to prevent poliomyelitis?

The polioviruses (PVs) are mainly transmitted by direct contact with an infected person through the fecal-oral route and respiratory secretions (or more rarely via contaminated water or food) and have a primary tropism for the gut. After their replication in the gut, in rare cases (far less than 1% of the infected individuals), PVs can spread to the central nervous system leading to flaccid paralysis, which can result in respiratory paralysis and death. By the middle of the 20th century, every year the wild polioviruses (WPVs) are supposed to have killed or paralyzed over half a million people. The introduction of the oral poliovirus vaccines (OPVs) through mass vaccination campaigns (combined with better application of hygiene measures), was a success story which enabled the World Health Organization (WHO) to set the global eradication of poliomyelitis as an objective. However this strategy of viral eradication has its limits as the majority of poliomyelitis cases today arise in individuals infected with circulating vaccine-derived polioviruses (cVDPVs) which regain pathogenicity following reversion or recombination. In recent years (between January 2018 and May 2023), the WHO recorded 8.8 times more cases of polio which were linked to the attenuated OPV vaccines (3,442 polio cases after reversion or recombination events) than cases linked to a WPV (390 cases). Recent knowledge of the evolution of RNA viruses and the exchange of genetic material among biological entities of the intestinal microbiota, call for a reassessment of the polio eradication vaccine strategies.

Defining Distinct RNA-Protein Interactomes of SARS-CoV-2 Genomic and Subgenomic RNAs

Host RNA binding proteins recognize viral RNA and play key roles in virus replication and antiviral mechanisms. SARS-CoV-2 generates a series of tiered subgenomic RNAs (sgRNAs), each encoding distinct viral protein(s) that regulate different aspects of viral replication. Here, for the first time, we demonstrate the successful isolation of SARS-CoV-2 genomic RNA and three distinct sgRNAs (N, S, and ORF8) from a single population of infected cells and characterize their protein interactomes. Over 500 protein interactors (including 260 previously unknown) were identified as associated with one or more target RNA. These included protein interactors unique to a single RNA pool and others present in multiple pools, highlighting our ability to discriminate between distinct viral RNA interactomes despite high sequence similarity. Individual interactomes indicated viral associations with cell response pathways, including regulation of cytoplasmic ribonucleoprotein granules and posttranscriptional gene silencing. We tested the significance of three protein interactors in these pathways (APOBEC3F, PPP1CC, and MSI2) using siRNA knockdowns, with several knockdowns affecting viral gene expression, most consistently PPP1CC. This study describes a new technology for high-resolution studies of SARS-CoV-2 RNA regulation and reveals a wealth of new viral RNA-associated host factors of potential functional significance to infection.

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.

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.

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.

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.

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 roles of nucleic acid editing in adaptation of zoonotic viruses to humans

Following spillover, viruses must adapt to new selection pressures exerted by antiviral responses in their new hosts. In mammals, cellular defense mechanisms often include viral nucleic acid editing pathways mediated through protein families apolipoprotein-B mRNA-editing complex (APOBEC) and Adenosine Deaminase Acting on ribonucleic acid (ADAR). APOBECs induce C→U transitions in viral genomes; the APOBEC locus is highly polymorphic with variable numbers of APOBEC3 paralogs and target preferences in humans and other mammals. APOBEC3 paralogs have shaped the evolutionary history of human immunodeficiency virus, with compelling bioinformatic evidence also for its mutagenic impact on monkeypox virus and severe acute respiratory syndrome coronavirus 2. ADAR-1 induces adenose-to-inosine (A→I) substitutions in double-stranded ribonucleic acid (RNA); its role in virus adaptation is less clear, as are epigenetic modifications to viral genomes, such as methylation. Nucleic acid editing restricts evolutionary space in which viruses can explore and may restrict viral-host range.

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.

Defining distinct RNA-protein interactomes of SARS-CoV-2 genomic and subgenomic RNAs

Host RNA binding proteins recognize viral RNA and play key roles in virus replication and antiviral defense mechanisms. SARS-CoV-2 generates a series of tiered subgenomic RNAs (sgRNAs), each encoding distinct viral protein(s) that regulate different aspects of viral replication. Here, for the first time, we demonstrate the successful isolation of SARS-CoV-2 genomic RNA and three distinct sgRNAs (N, S, and ORF8) from a single population of infected cells and characterize their protein interactomes. Over 500 protein interactors (including 260 previously unknown) were identified as associated with one or more target RNA at either of two time points. These included protein interactors unique to a single RNA pool and others present in multiple pools, highlighting our ability to discriminate between distinct viral RNA interactomes despite high sequence similarity. The interactomes indicated viral associations with cell response pathways including regulation of cytoplasmic ribonucleoprotein granules and posttranscriptional gene silencing. We validated the significance of five protein interactors predicted to exhibit antiviral activity (APOBEC3F, TRIM71, PPP1CC, LIN28B, and MSI2) using siRNA knockdowns, with each knockdown yielding increases in viral production. This study describes new technology for studying SARS-CoV-2 and reveals a wealth of new viral RNA-associated host factors of potential functional significance to infection.

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.

Evolution of the SARS-CoV-2 Mutational Spectrum

SARS-CoV-2 evolves rapidly in part because of its high mutation rate. Here, we examine whether this mutational process itself has changed during viral evolution. To do this, we quantify the relative rates of different types of single-nucleotide mutations at 4-fold degenerate sites in the viral genome across millions of human SARS-CoV-2 sequences. We find clear shifts in the relative rates of several types of mutations during SARS-CoV-2 evolution. The most striking trend is a roughly 2-fold decrease in the relative rate of G→T mutations in Omicron versus early clades, as was recently noted by Ruis et al. (2022. Mutational spectra distinguish SARS-CoV-2 replication niches. bioRxiv, doi:10.1101/2022.09.27.509649). There is also a decrease in the relative rate of C→T mutations in Delta, and other subtle changes in the mutation spectrum along the phylogeny. We speculate that these changes in the mutation spectrum could arise from viral mutations that affect genome replication, packaging, and antagonization of host innate-immune factors, although environmental factors could also play a role. Interestingly, the mutation spectrum of Omicron is more similar than that of earlier SARS-CoV-2 clades to the spectrum that shaped the long-term evolution of sarbecoviruses. Overall, our work shows that the mutation process is itself a dynamic variable during SARS-CoV-2 evolution and suggests that human SARS-CoV-2 may be trending toward a mutation spectrum more similar to that of other animal sarbecoviruses.

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.

The P323L substitution in the SARS-CoV-2 polymerase (NSP12) confers a selective advantage during infection

BACKGROUND

The mutational landscape of SARS-CoV-2 varies at the dominant viral genome sequence and minor genomic variant population. During the COVID-19 pandemic, an early substitution in the genome was the D614G change in the spike protein, associated with an increase in transmissibility. Genomes with D614G are accompanied by a P323L substitution in the viral polymerase (NSP12). However, P323L is not thought to be under strong selective pressure.

RESULTS

Investigation of P323L/D614G substitutions in the population shows rapid emergence during the containment phase and early surge phase during the first wave. These substitutions emerge from minor genomic variants which become dominant viral genome sequence. This is investigated in vivo and in vitro using SARS-CoV-2 with P323 and D614 in the dominant genome sequence and L323 and G614 in the minor variant population. During infection, there is rapid selection of L323 into the dominant viral genome sequence but not G614. Reverse genetics is used to create two viruses (either P323 or L323) with the same genetic background. L323 shows greater abundance of viral RNA and proteins and a smaller plaque morphology than P323.

CONCLUSIONS

These data suggest that P323L is an important contribution in the emergence of variants with transmission advantages. Sequence analysis of viral populations suggests it may be possible to predict the emergence of a new variant based on tracking the frequency of minor variant genomes. The ability to predict an emerging variant of SARS-CoV-2 in the global landscape may aid in the evaluation of medical countermeasures and non-pharmaceutical interventions.

Oligonucleotide usage in coronavirus genomes mimics that in exon regions in host genomes

BACKGROUND

Viruses use various host factors for their growth, and efficient growth requires efficient use of these factors. Our previous study revealed that the occurrence frequency of oligonucleotides in the influenza virus genome is distinctly different among derived hosts, and the frequency tends to adapt to the host cells in which they grow. We aimed to study the adaptation mechanisms of a zoonotic virus to host cells.

METHODS

Herein, we compared the frequency of oligonucleotides in the genome of alpha- and betacoronavirus with those in the genomes of humans and bats, which are typical hosts of the viruses.

RESULTS

By comparing the oligonucleotide frequency in coronaviruses and their host genomes, we found a statistically tested positive correlation between the frequency of coronaviruses and that of the exon regions of the host from which the virus is derived. To examine the characteristics of early-stage changes in the viral genome, which are assumed to accompany the host change from non-humans to humans, we compared the oligonucleotide frequency between severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) at the beginning of the pandemic and the prevalent variants thereafter, and found changes towards the frequency of the host exon regions.

CONCLUSIONS

In alpha- and betacoronaviruses, the genome oligonucleotide frequency is thought to change in response to the cellular environment in which the virus is replicating, and actually the frequency has approached the frequency in exon regions in the host.

Perceptions into causes and consequences of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) variants

Severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) variants are emerging worldwide and pathogenicity varies widely from no symptoms to death. The SARS‐CoV‐2 is evolving as lineages like Alpha, Beta, Gamma, Epsilon, Iota, Delta, and Omicron in the course of time. The main reasons for such viral evolution are (a) the imperfect nature of SARS‐CoV‐2 RNA polymerase, and viral exonuclease mediated proofreading functions resulting in the generation of mutations in viral genomes; (b) fusions of the 5′ leader sequence to unexpected 3′ sites, and transcription regulatory sequences (TRSs) in subgenomic RNAs (sgRNAs), which result in the generation of structural variants and novel open reading frames; (c) these viruses are combated by the host type I interferons (IFNs). In such a process IFNs upregulate viral RNA editing APOBEC3G/F and ADAR1 genes, which induce mutations in viral genomes. These factors play important roles in causing viral evolution and the emergence of more efficient SARS‐CoV‐2 genomes, which escape the host immune defense system, and vaccine‐elicited antibodies and impede therapeutic strategies. The main challenges we now face are how to control future SARS‐CoV‐2 evolution, the elimination of their deleterious side effects, and the onset of new diseases as aftermaths of SARS‐CoV‐2 infections. Preventive measures like (a) the development of broadly neutralizing antibodies and novel vaccines, therapies based on genomics and proteomics data will help in avoiding, and/or minimizing SARS‐CoV‐2 infections; (b) targeted therapies, application of patient‐based precision medicine methodology can help in achieving the goal and avoiding unwanted deleterious side effects and the onset of SARS‐CoV‐2 infections mediated several diseases in future.

Role of error catastrophe in transmission ability of virus

The role played by error catastrophe is explicitly taken into account in a mathematical formulation to analyze COVID-19 data. The idea is to combine the mathematical genetics formalism of the error catastrophe of mutations in virus gene loci with the standard model of epidemics, which lacks the explicit incorporation of the effect of mutation on the spreading of viruses. We apply this formalism to the case of SARS-CoV-2 virus. We assume the universality of the error catastrophe in the process of analyzing the data. This means that some basic parameter to describe the error catastrophe is independent of which group (country or city) we deal with. Concretely, we analyze Omicron variant data from South Africa and then analyze cases from Japan using the same value of the basic parameter derived in the South Africa analysis. The excellent fit between the two sets of data, one from South Africa and the other from Japan, using the common values of genetic parameters, justifies our assumption of the universality of these parameters.

A Long-Running Arms Race between APOBEC1 Genes and Retroviruses in Tetrapods

Activation-induced cytidine deaminase/apolipoprotein B mRNA editing catalytic polypeptide-like (AID/APOBEC) proteins are cytosine deaminases implicated in diverse biological functions. APOBEC1 (A1) proteins have long been thought to regulate lipid metabolism, whereas the evolutionary significance of A1 proteins in antiviral defense remains largely obscure. Endogenous retroviruses (ERVs) document past retroviral infections and are ubiquitous within the vertebrate genomes. Here, we identify the A1 gene repertoire, characterize the A1-mediated mutation footprints in ERVs, and interrogate the evolutionary arms race between A1 genes and ERVs across vertebrate species. We find that A1 genes are widely present in tetrapods, recurrently amplified and lost in certain lineages, suggesting that A1 genes might have originated during the early evolution of tetrapods. A1-mediated mutation footprints can be detected in ERVs across tetrapods. Moreover, A1 genes appear to have experienced episodic positive selection in many tetrapod lineages. Taken together, we propose that a long-running arms race between A1 genes and retroviruses might have persisted throughout the evolutionary course of tetrapods. IMPORTANCE APOBEC3 (A3) genes have been thought to function in defense against retroviruses, whereas the evolutionary significance of A1 proteins in antiviral defense remains largely obscure. In this study, we identify the A1 gene repertoire, characterize the A1-mediated mutation footprints in endogenous retroviruses (ERVs), and explore the evolutionary arms race between A1 genes and ERVs across vertebrate species. We found A1 proteins originated during the early evolution of tetrapods, and detected the footprints of A1-induced hypermutations in retroviral fossils. A1 genes appear to have experienced pervasive positive selection in tetrapods. Our study indicates a long-running arms race between A1 genes and retroviruses taking place throughout the evolutionary course of tetrapods.

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.

Convergent Evolution in SARS-CoV-2 Spike Creates a Variant Soup from Which New COVID-19 Waves Emerge

The first 2 years of the COVID-19 pandemic were mainly characterized by recurrent mutations of SARS-CoV-2 Spike protein at residues K417, L452, E484, N501 and P681 emerging independently across different variants of concern (Alpha, Beta, Gamma, and Delta). Such homoplasy is a marker of convergent evolution. Since Spring 2022 and the third year of the pandemic, with the advent of Omicron and its sublineages, convergent evolution has led to the observation of different lineages acquiring an additional group of mutations at different amino acid residues, namely R346, K444, N450, N460, F486, F490, Q493, and S494. Mutations at these residues have become increasingly prevalent during Summer and Autumn 2022, with combinations showing increased fitness. The most likely reason for this convergence is the selective pressure exerted by previous infection- or vaccine-elicited immunity. Such accelerated evolution has caused failure of all anti-Spike monoclonal antibodies, including bebtelovimab and cilgavimab. While we are learning how fast coronaviruses can mutate and recombine, we should reconsider opportunities for economically sustainable escape-proof combination therapies, and refocus antibody-mediated therapeutic efforts on polyclonal preparations that are less likely to allow for viral immune escape.

Control of Human Anelloviruses by Cytosine to Uracil Genome Editing

Anelloviruses are the most common viruses infecting humans. Every human carries a nonpathogenic personal anellovirus virome (anellome), yet it is unknown which mechanisms contribute to its stability. Here, we assessed the dynamics and impact of a host antiviral defense mechanism-cytidine deaminase activity leading to C to U editing in anelloviruses-on the stability of the anellome. We investigated anellome sequence data obtained from serum samples collected every 6 months from two healthy subjects followed for more than 30 years. The subjects were infected by a total of 64 anellovirus lineages. Minus-stranded C to U editing was observed in lineages belonging to the Alpha-, Beta-, and Gammatorquevirus genera. The edited genomes were present within virus particles, therefore editing must have occurred at the late stages of the virus life cycle. Editing was favored by 5'-TC contexts in the virus genome, indicating that apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like, catalytic subunit 3 or A3 (APOBEC3) proteins are involved. Within a lineage, mutational dynamics varied over time and few fixations of mutations were detected, indicating that C to U editing is a dead end for a virus genome. We detected an editing coldspot in the GC-rich regions, suggesting that the GC-rich region is crucial for genome packaging, since only packaged virus particles were included in the analysis. Finally, we noticed a lineage-specific reduced concentration after an editing event, yet no clearance. In conclusion, cytidine deaminase activity does not clear anelloviruses, nor does it play a major role in virus evolution, but it does contribute to the stability of the anellome. IMPORTANCE Despite significant attention on anellovirus research, the interaction between the anellovirus virome and the human host remains unknown. We show the dynamics of APOBEC3-mediated cytidine deaminase activity on anelloviruses during a 30-year period of chronic infection and postulate that this antiviral mechanism controls anelloviruses. These results expand our knowledge of anellovirus-host interactions, which may be important for the design of gene therapies.

Characterization of SARS-CoV-2 Mutational Signatures from 1.5+ Million Raw Sequencing Samples

We present a large-scale analysis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) substitutions, considering 1,585,456 high-quality raw sequencing samples, aimed at investigating the existence and quantifying the effect of mutational processes causing mutations in SARS-CoV-2 genomes when interacting with the human host. As a result, we confirmed the presence of three well-differentiated mutational processes likely ruled by reactive oxygen species (ROS), apolipoprotein B editing complex (APOBEC), and adenosine deaminase acting on RNA (ADAR). We then evaluated the activity of these mutational processes in different continental groups, showing that some samples from Africa present a significantly higher number of substitutions, most likely due to higher APOBEC activity. We finally analyzed the activity of mutational processes across different SARS-CoV-2 variants, and we found a significantly lower number of mutations attributable to APOBEC activity in samples assigned to the Omicron variant.

Molecular insights into onset of autoimmunity in SARS-CoV-2 infected patients

Some of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infected patients are facing long-term devastating effects like induction of autoimmune diseases. Here, I discuss molecular mechanisms and risk factors involved in the induction of autoimmune diseases after SARS-CoV-2 infections. Transcript editing genes were upregulated during SARS-CoV-2 infections, which might have edited host gene transcripts and paved the way for autoantigens generation and presented as nonself to generate autoantibodies followed by auto immunogenicity after SARS-CoV-2 infections. Therefore, some SARS-CoV-2 patients acquire autoimmunity. The transient and/or innocuous autoimmune response in some SARS-CoV-2 infected patients may be due to a lack of repeated production of autoantibodies to host autoantigens and/or viral antigens, which are needed to boost autoimmune response. In the future, SARS-CoV-2 mediated autoimmune disease onset will be a challenging task. Therefore, possible preventive measures and strategies to minimize and/or preclude such SARS-CoV-2 mediated autoimmune diseases have been presented in this commentary.

Cellular miR-6741-5p as a Prognostic Biomarker Predicting Length of Hospital Stay among COVID-19 Patients

Wide variability exists with host response to SARS-CoV-2 infection among individuals. Circulatory micro RNAs (miRNAs) are being recognized as promising biomarkers for complex traits, including viral pathogenesis. We hypothesized that circulatory miRNAs at 48 h post hospitalization may predict the length of stay (LOS) and prognosis of COVID-19 patients. Plasma miRNA levels were compared between three groups: (i) healthy volunteers (C); (ii) COVID-19 patients treated with remdesivir (an antiviral) plus dexamethasone (a glucocorticoid) (with or without baricitinib, a Janus kinase inhibitor) on the day of hospitalization (I); and COVID-19 patients at 48 h post treatment (T). Results showed that circulatory miR-6741-5p expression levels were significantly different between groups C and I (p < 0.0000001); I and T (p < 0.0000001); and C and T (p = 0.001). Our ANOVA model estimated that all patients with less than 12.42 Log2 CPM had a short LOS, or a good prognosis, whereas all patients with over 12.42 Log2 CPM had a long LOS, or a poor prognosis. In sum, we show that circulatory miR-6741-5p may serve as a prognostic biomarker effectively predicting mortality risk and LOS of hospitalized COVID-19 patients.

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.

The coevolution between APOBEC3 and retrotransposons in primates

Retrotransposons are genetic elements with the ability to replicate in the genome using reverse transcriptase: they have been associated with the development of different biological structures, such as the Central Nervous System (CNS), and their high mutagenic potential has been linked to various diseases, including cancer and neurological disorders. Throughout evolution and over time, Primates and Homo had to cope with infections from viruses and bacteria, and also with endogenous retroelements. Therefore, host genomes have evolved numerous methods to counteract the activity of endogenous and exogenous pathogens, and the APOBEC3 family of mutators is a prime example of a defensive mechanism in this context.In most Primates, there are seven members of the APOBEC3 family of deaminase proteins: among their functions, there is the ability to inhibit the mobilization of retrotransposons and the functionality of viruses. The evolution of the APOBEC3 proteins found in Primates is correlated with the expansion of two major families of retrotransposons, i.e. ERV and LINE-1.In this review, we will discuss how the rapid expansion of the APOBEC3 family is linked to the evolution of retrotransposons, highlighting the strong evolutionary arms race that characterized the history of APOBEC3s and endogenous retroelements in Primates. Moreover, the possible role of this relationship will be assessed in the context of embryonic development and brain-associated diseases.

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.