Evidence for host-dependent RNA editing in the transcriptome of SARS-CoV-2

SARS-CoV-2 point mutations are over-represented in terminal loops of RNA stem-loop structures that can be resolved by Nsp13 helicase in a unique manner with respect to nucleotide dependence

To improve health outcomes for COVID-19 (coronavirus disease 2019) patients, the factors that influence coronavirus genome variation need to be ascertained. The SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) genome is rich in predicted RNA secondary structures, particularly stem-loops (SLs) formed by intramolecular base pairing within palindromic sequences. We analyzed the NCBI Virus collection of SARS-CoV-2 genome sequences from COVID-19 individuals to map variants relative to SL structural elements. Point mutations in the SARS-CoV-2 genome, with a C-to-U transition bias, were over-represented in unpaired nucleotides and, more specifically, within the terminal loops of RNA SL structures. As the sole helicase encoded by SARS-CoV-2, Nsp13 may operate in the timely resolution of secondary RNA structures to facilitate SARS-CoV-2 RNA copying or processing. We characterized Nsp13 to resolve SARS-CoV-2 sequence-derived unimolecular RNA SL substrates and determined that it does so in a functionally cooperative manner. In addition to ATP, Nsp13 resolves the unimolecular RNA SL structure in the absence of nucleotide, in contrast to the strict ATP requirement for a bimolecular RNA forked duplex. We suggest a model in which a series of binary and ternary complex interactions of Nsp13 with nucleotide and/or RNA SL pose mechanistic implications for RNA SL resolution.

Unveiling RNA Editing by ADAR and APOBEC Protein Gene Families

RNA editing is a crucial post-transcriptional modification that alters the transcriptome and proteome and affects many cellular processes, including splicing, microRNA specificity, stability of RNA molecules, and protein structure. Enzymes from the adenosine deaminase acting on RNA (ADAR) and apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) protein families mediate RNA editing and can alter a variety of non-coding and coding RNAs, including all regions of mRNA molecules, leading to tumor development and progression. This review provides novel insights into the potential use of RNA editing parameters, such as editing levels, expression of ADAR and APOBEC genes, and specifically edited genes, as biomarkers for cancer progression, distinguishing it from previous studies that focused on isolated aspects of RNA editing mechanisms. The methodological section offers clues to accelerate high-throughput analysis of RNA or DNA sequencing data for the identification of RNA editing events.

Within-Host Fitness and Antigenicity Shift Are Key Factors Influencing the Prevalence of Within-Host Variations in the SARS-CoV-2 S Gene

Within-host evolution plays a critical role in shaping the diversity of SARS-CoV-2. However, understanding the primary factors contributing to the prevalence of intra-host single nucleotide variants (iSNVs) in the viral population remains elusive. Here, we conducted a comprehensive analysis of over 556,000 SARS-CoV-2 sequencing data and prevalence data of different SARS-CoV-2 S protein amino acid mutations to elucidate key factors influencing the prevalence of iSNVs in the SARS-CoV-2 S gene. Within-host diversity analysis revealed the presence of mutational hotspots within the S gene, mainly located in NTD, RBD, TM, and CT domains. Additionally, we generated a single amino acid resolution selection status map of the S protein. We observed a significant variance in within-host fitness among iSNVs in the S protein. The majority of iSNVs exhibited low to no within-host fitness and displayed low alternate allele frequency (AAF), suggesting that they will be eliminated due to the narrow transmission bottleneck of SARS-CoV-2. Notably, iSNVs with moderate AAFs (0.06-0.12) were found to be more prevalent than those with high AAFs. Furthermore, iSNVs with the potential to alter antigenicity were more prevalent. These findings underscore the significance of within-host fitness and antigenicity shift as two key factors influencing the prevalence of iSNVs in the SARS-CoV-2 S gene.

APOBEC3-Related Editing and Non-Editing Determinants of HIV-1 and HTLV-1 Restriction

The apolipoprotein B mRNA editing enzyme catalytic polypeptide-like 3 (APOBEC3/A3) family of cytosine deaminases serves as a key innate immune barrier against invading retroviruses and endogenous retroelements. The A3 family's restriction activity against these parasites primarily arises from their ability to catalyze cytosine-to-uracil conversions, resulting in genome editing and the accumulation of lethal mutations in viral genomes. Additionally, non-editing mechanisms, including deaminase-independent pathways, such as blocking viral reverse transcription, have been proposed as antiviral strategies employed by A3 family proteins. Although viral factors can influence infection progression, the determinants that govern A3-mediated restriction are critical in shaping retroviral infection outcomes. This review examines the interactions between retroviruses, specifically human immunodeficiency virus type 1 and human T-cell leukemia virus type 1, and A3 proteins to better understand how editing and non-editing activities contribute to the trajectory of these retroviral infections.

An Expanding Universe of Mutational Signatures and Its Rapid Evolution in Single-Stranded RNA Viruses

The study of mutational processes in somatic genomes has gained recent momentum, uncovering a wide array of endogenous and exogenous factors associated with somatic changes. However, the overall landscape of mutational processes in germline mutations across the tree of life and associated evolutionary driving forces are rather unclear. In this study, we analyzed mutational processes in single-stranded RNA (ssRNA) viruses which are known to jump between different hosts with divergent exogenous environments. We found that mutational spectra in different ssRNA viruses differ significantly and are mainly associated with their genetic divergence. Surprisingly, host environments contribute much less significantly to the mutational spectrum, challenging the prevailing view that the exogenous cellular environment is a major determinant of the mutational spectrum in viruses. To dissect the evolutionary forces shaping viral spectra, we selected two important scenarios, namely the inter-host evolution between different viral strains as well as the intra-host evolution. In both scenarios, we found mutational spectra change significantly through space and time, strongly correlating with levels of natural selection. Combining the mutations across all ssRNA viruses, we identified a suite of mutational signatures with varying degrees of similarity to somatic signatures in humans, indicating universal and divergent mutational processes across the tree of life. Taken together, we unraveled an unprecedented dynamic landscape of mutational processes in ssRNA viruses, pinpointing important evolutionary forces shaping fast evolution of mutational spectra in different species.

The mutation rate of SARS-CoV-2 is highly variable between sites and is influenced by sequence context, genomic region, and RNA structure

RNA viruses like SARS-CoV-2 have a high mutation rate, which contributes to their rapid evolution. The rate of mutations depends on the mutation type (e.g., A→C, A→G, etc.) and can vary between sites in the viral genome. Understanding this variation can shed light on the mutational processes at play, and is crucial for quantitative modeling of viral evolution. Using the millions of available SARS-CoV-2 full-genome sequences, we estimate rates of synonymous mutations for all 12 possible nucleotide mutation types and examine how much these rates vary between sites. We find a surprisingly high level of variability and several striking patterns: the rates of four mutation types suddenly increase at one of two gene boundaries; the rates of most mutation types strongly depend on a site's local sequence context, with up to 56-fold differences between contexts; consistent with a previous study, the rates of some mutation types are lower at sites engaged in RNA secondary structure. A simple log-linear model of these features explains ~15-60% of the fold-variation of mutation rates between sites, depending on mutation type; more complex models only modestly improve predictive power out of sample. We estimate the fitness effect of each mutation based on the number of times it actually occurs versus the number of times it is expected to occur based on the model. We identify several small regions of the genome where synonymous or noncoding mutations occur much less often than expected, indicative of strong purifying selection on the RNA sequence that is independent of protein sequence. Overall, this work expands our basic understanding of SARS-CoV-2's evolution by characterizing the virus's mutation process at the level of individual sites and uncovering several striking mutational patterns that arise from unknown mechanisms.

SARS-CoV-2 introductions to the island of Ireland: a phylogenetic and geospatiotemporal study of infection dynamics

BACKGROUND

Ireland's COVID-19 response combined extensive SARS-CoV-2 testing to estimate incidence, with whole genome sequencing (WGS) for genome surveillance. As an island with two political jurisdictions-Northern Ireland (NI) and Republic of Ireland (RoI)-and access to detailed passenger travel data, Ireland provides a unique setting to study virus introductions and evaluate public health measures. Using a substantial Irish genomic dataset alongside global data from GISAID, this study aimed to trace the introduction and spread of SARS-CoV-2 across the island.

METHODS

We recursively searched for 29,518 SARS-CoV-2 genome sequences collected in Ireland from March 2020 to June 2022 within the global SARS-CoV-2 phylogenetic tree and identified clusters based on shared last common non-Irish ancestors. A maximum parsimony approach was used to assign a likely country of origin to each cluster. The geographic locations and collection dates of the samples in each introduction cluster were used to map the spread of the virus across Ireland. Downsampling was used to model the impact of varying levels of sequencing and normalisation for population permitted comparison between jurisdictions.

RESULTS

Six periods spanning the early introductions and the emergence of Alpha, Delta, and Omicron variants were studied in detail. Among 4439 SARS-CoV-2 introductions to Ireland, 2535 originated in England, with additional cases largely from the rest of Great Britain, United States of America, and Northwestern Europe. Introduction clusters ranged in size from a single to thousands of cases. Introductions were concentrated in the densely populated Dublin and Belfast areas, with many clusters spreading islandwide. Genetic phylogeny was able to effectively trace localised transmission patterns. Introduction rates were similar in NI and RoI for most variants, except for Delta, which was more frequently introduced to NI.

CONCLUSIONS

Tracking individual introduction events enables detailed modelling of virus spread patterns and clearer assessment of the effectiveness of control measures. Stricter travel restrictions in RoI likely reduced Delta introductions but not infection rates, which were similar across jurisdictions. Local and global sequencing levels influence the information available from phylogenomic analyses and we describe an approach to assess the ability of a chosen WGS level to detect virus introductions.

Past, Present, and Future of RNA Modifications in Infectious Disease Research

In early 2024, the National Academies of Sciences, Engineering, and Medicine (NASEM) released a roadmap for the future of research into mapping ribonucleic acid (RNA) modifications, which underscored the importance of better defining these diverse chemical changes to the RNA macromolecule. As nearly all mature RNA molecules harbor some form of modification, we must understand RNA modifications to fully appreciate the functionality of RNA. The NASEM report calls for massive mobilization of resources and investment akin to the transformative Human Genome Project of the early 1990s. Like the Human Genome Project, a concerted effort in improving our ability to assess every single modification on every single RNA molecule in an organism will change the way we approach biological questions, accelerate technological advance, and improve our understanding of the molecular world. Consequently, we are also at the start of a revolution in defining the impact of RNA modifications in the context of host-microbe and even microbe-microbe interactions. In this perspective, we briefly introduce RNA modifications to the infection biologist, highlight key aspects of the NASEM report and exciting examples of RNA modifications contributing to host and pathogen biology, and finally postulate where infectious disease research may benefit from this exciting new endeavor in globally mapping RNA modifications.

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.

Viral N protein hijacks deaminase-containing RNA granules to enhance SARS-CoV-2 mutagenesis

Host cell-encoded deaminases act as antiviral restriction factors to impair viral replication and production through introducing mutations in the viral genome. We sought to understand whether deaminases are involved in SARS-CoV-2 mutation and replication, and how the viral factors interact with deaminases to trigger these processes. Here, we show that APOBEC and ADAR deaminases act as the driving forces for SARS-CoV-2 mutagenesis, thereby blocking viral infection and production. Mechanistically, SARS-CoV-2 nucleocapsid (N) protein, which is responsible for packaging viral genomic RNA, interacts with host deaminases and co-localizes with them at stress granules to facilitate viral RNA mutagenesis. N proteins from several coronaviruses interact with host deaminases at RNA granules in a manner dependent on its F17 residue, suggesting a conserved role in modulation of viral mutagenesis in other coronaviruses. Furthermore, mutant N protein bearing a F17A substitution cannot localize to deaminase-containing RNA granules and leads to reduced mutagenesis of viral RNA, providing support for its function in enhancing deaminase-dependent viral RNA editing. Our study thus provides further insight into virus-host cell interactions mediating SARS-CoV-2 evolution.

SARS-CoV-2 disrupts host gene networks: Unveiling key hub genes as potential therapeutic targets for COVID-19 management

PURPOSE

Although the end of COVID-19 as a public health emergency was declared on May 2023, still new cases of the infection are reported and the risk remains of new variants emerging that may cause new surges in cases and deaths. While clinical symptoms have been rapidly defined worldwide, the basic body responses and pathogenetic mechanisms acting in patients with SARS-CoV-2 infection over time until recovery or death require further investigation. The understanding of the molecular mechanisms underlying the development and course of the disease is essential in designing effective preventive and therapeutic approaches, and ultimately reducing mortality and disease spreading.

METHODS

The current investigation aimed to identify the key genes engaged in SARS-CoV-2 infection. To achieve this goal high-throughput RNA sequencing of peripheral blood samples collected from healthy donors and COVID-19 patients was performed. The resulting sequence data were processed using a wide range of bioinformatics tools to obtain detailed modifications within five transcriptomic phenomena: expression of genes and long non-coding RNAs, alternative splicing, allel-specific expression and circRNA production. The in silico procedure was completed with a functional analysis of the identified alterations.

RESULTS

The transcriptomic analysis revealed that SARS-CoV-2 has a significant impact on multiple genes encoding ribosomal proteins (RPs). Results show that these genes differ not only in terms of expression but also manifest biases in alternative splicing and ASE ratios. The integrated functional analysis exposed that RPs mostly affected pathways and processes related to infection-COVID-19 and NOD-like receptor signaling pathway, SARS-CoV-2-host interactions and response to the virus. Furthermore, our results linked the multiple intronic ASE variants and exonic circular RNA differentiations with SARS-CoV-2 infection, suggesting that these molecular events play a crucial role in mRNA maturation and transcription during COVID-19 disease.

CONCLUSIONS

By elucidating the genetic mechanisms induced by the virus, the current research provides significant information that can be employed to create new targeted therapeutic strategies for future research and treatment related to COVID-19. Moreover, the findings highlight potentially promising therapeutic biomarkers for early risk assessment of critically ill patients.

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 .

MethylCallR : a comprehensive analysis framework for Illumina Methylation Beadchip

DNA methylation is a molecular process that mediates gene-environment interactions. Epigenome-wide association studies (EWAS) using the Illumina Human Methylation BeadChip are powerful tools for quantifying the relationship between DNA methylation and phenotypes. Recently, the Illumina Methylation EPICv2 BeadChip (EPICv2) was released, which includes new features, such as duplicated probes and changed probe names. Several published algorithms have been updated to address these features in EPICv2. However, appropriate EPICv2 preprocessing and integration with previous microarray versions remain complex. Therefore, MethylCallR, an open-source R package designed to provide standard procedures for performing EWAS using Illumina methylation microarrays including EPICv2, was developed. MethylCallR can be used to control duplicated probes in EPICv2, by using pre-set data implemented in MethylCallR or new customized data. MethylCallR includes a straightforward conversion function between different types of Illumina Human Methylation BeadChips. Using MethylCallR, potential outlier sample detection and statistical power estimation were conducted and used to select meaningful probes. Publicly available data was analyzed using MethylCallR and the findings were compared to that of a previous study.

Addressing pandemic-wide systematic errors in the SARS-CoV-2 phylogeny

The SARS-CoV-2 genome occupies a unique place in infection biology - it is the most highly sequenced genome on earth (making up over 20% of public sequencing datasets) with fine scale information on sampling date and geography, and has been subject to unprecedented intense analysis. As a result, these phylogenetic data are an incredibly valuable resource for science and public health. However, the vast majority of the data was sequenced by tiling amplicons across the full genome, with amplicon schemes that changed over the pandemic as mutations in the viral genome interacted with primer binding sites. In combination with the disparate set of genome assembly workflows and lack of consistent quality control (QC) processes, the current genomes have many systematic errors that have evolved with the virus and amplicon schemes. These errors have significant impacts on the phylogeny, and therefore over the last few years, many thousands of hours of researchers time has been spent in "eyeballing" trees, looking for artefacts, and then patching the tree. Given the huge value of this dataset, we therefore set out to reprocess the complete set of public raw sequence data in a rigorous amplicon-aware manner, and build a cleaner phylogeny. Here we provide a global tree of 4,471,579 samples, built from a consistently assembled set of high quality consensus sequences from all available public data as of June 2024, viewable at https://viridian.taxonium.org. Each genome was constructed using a novel assembly tool called Viridian (https://github.com/iqbal-lab-org/viridian), developed specifically to process amplicon sequence data, eliminating artefactual errors and mask the genome at low quality positions. We provide simulation and empirical validation of the methodology, and quantify the improvement in the phylogeny. We hope the tree, consensus sequences and Viridian will be a valuable resource for researchers.

Mass spectrometry-based metabolomics reveals metabolism of molnupiravir may lead to metabolic disorders and hepatotoxicity

Molnupiravir (MO) is a pyrimidine nucleoside anti-SARS-CoV-2 drug. MO treatment could cause mild liver injury. However, the underlying mechanism of MO-induced liver injury and the metabolic pathway of MO in vivo are unclear. In this study, metabolomics analysis and molecular biology methods were used to explore these issues. Through metabolomics analysis, it was found that the homeostasis of pyrimidine, purine, lysophosphatidylcholine (LPC), and amino acids in mice was destroyed after MO treatment. A total of 80 changed metabolites were detected. Among these changed metabolites, 4-ethylphenyl sulfate, dihydrouracil, and LPC 20:0 was related to the elevation of alkaline phosphatase (ALP), interleukin-6 (IL6), and nuclear factor kappa-B (NF-κB). The levels of 4-ethylphenyl sulfate, dihydrouracil, and LPC 20:0 in plasma were positively correlated with their levels in the liver, suggesting that these metabolites were associated with MO-induced liver injury. MO treatment could increase NHC and cytidine levels, activate cytidine deaminase (CDA), and increase LPC levels. CDA and LPC could increase the mRNA expression level of toll-like receptor (TLR). The current study indicated that the elevation of hepatic TLR may be an important reason for MO leading to the liver injury.

SARS-CoV-2 within-host population expansion, diversification and adaptation in zoo tigers, lions and hyenas

SARS-CoV-2 rapidly adapts to new hosts following cross-species transmission; this is highly relevant as novel within-host variants have emerged following infection of susceptible wild and domestic animal species. Furthermore, SARS-CoV-2 transmission from animals (e.g., white-tailed deer, mink, domestic cats, and others) back to humans has also been observed, documenting the potential of novel animal-derived variants to infect humans. We investigated SARS-CoV-2 evolution and host-specific adaptation during an outbreak in Amur tigers (Panthera tigris altaica), African lions (Panthera leo), and spotted hyenas (Crocuta crocuta) at Denver Zoo in late 2021. SARS-CoV-2 genomes from longitudinal samples collected from 16 individuals were evaluated for within-host variation and genomic signatures of selection. The outbreak was likely initiated by a single spillover of a rare Delta sublineage subsequently transmitted from tigers to lions to hyenas. Within-host virus populations rapidly expanded and diversified. We detected signatures of purifying and positive selection, including strong positive selection in hyenas and in the nucleocapsid (N) gene in all animals. Four candidate species-specific adaptive mutations were identified: N A254V in lions and hyenas, and ORF1a E1724D, spike T274I, and N P326L in hyenas. These results reveal accelerated SARS-CoV-2 adaptation following host shifts in three non-domestic species in daily contact with humans.

Genomic epidemiology reveals the variation and transmission properties of SARS-CoV-2 in a single-source community outbreak

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused the coronavirus disease 2019 (COVID-19) pandemic, which is still a global public health concern. During March 2022, a rapid and confined single-source outbreak of SARS-CoV-2 was identified in a community in Nanjing municipal city. Overall, 95 individuals had laboratory-confirmed SARS-CoV-2 infection. The whole genomes of 61 viral samples were obtained, which were all members of the BA.2.2 lineage and clearly demonstrated the presence of one large clade, and all the infections could be traced back to the original index case. The most distant sequence from the index case presented a difference of 4 SNPs, and 118 intrahost single-nucleotide variants (iSNVs) at 74 genomic sites were identified. Some minor iSNVs can be transmitted and subsequently rapidly fixed in the viral population. The minor iSNVs transmission resulted in at least two nucleotide substitutions among all seven SNPs identified in the outbreak, generating genetically diverse populations. We estimated the overall transmission bottleneck size to be 3 using 11 convincing donor-recipient transmission pairs. Our study provides new insights into genomic epidemiology and viral transmission, revealing how iSNVs become fixed in local clusters, followed by viral transmission across the community, which contributes to population diversity.

Host-dependent C-to-U RNA editing in SARS-CoV-2 creates novel viral genes with optimized expressibility

Rampant C-to-U RNA editing drives the mutation and evolution of SARS-CoV-2. While much attention has been paid to missense mutations, the C-to-U events leading to AUG and thus creating novel ORFs were uninvestigated. By utilizing the public time-course mutation data from the worldwide SARS-CoV-2 population, we systematically identified the "AUG-gain mutations" caused by C-to-U RNA editing. Synonymous mutations were of special focus. A total of 58 synonymous C-to-U sites are able to create out-of-frame AUG in coding sequence (CDS). These 58 synonymous sites showed significantly higher allele frequency (AF) and increasing rate (dAF/dt) than other C-to-U synonymous sites in the SARS-CoV-2 population, suggesting that these 58 AUG-gain events conferred additional benefits to the virus and are subjected to positive selection. The 58 predicted new ORFs created by AUG-gain events showed the following advantages compared to random expectation: they have longer lengths, higher codon adaptation index (CAI), higher Kozak scores, and higher tRNA adaptation index (tAI). The 58 putatively novel ORFs have high expressibility and are very likely to be functional, providing an explanation for the positive selection on the 58 AUG-gain mutations. Our study proposed a possible mechanism of the emergence of de novo genes in SARS-CoV-2. This idea should be helpful in studying the mutation and evolution of SARS-CoV-2.

Host response to Aplysia Abyssovirus 1 in nervous system and gill

The California sea hare (Aplysia californica) is a model for age associated cognitive decline. Recent researched identified a novel nidovirus, Aplysia Abyssovirus 1, with broad tropism enriched in the Aplysia nervous system. This virus is ubiquitous in wild and maricultured, young and old animals without obvious pathology. Here we re-evaluated gene expression data from several previous studies to investigate differential expression in the nervous system and gill in response to virus and aging as well as the mutational spectrum observed in the viral sequences obtained from these datasets. Viral load and age were highly correlated, indicating persistent infection. Upregulated genes in response to virus were enriched for immune genes and signatures of ER and proteostatic stress, while downregulated genes were enriched for mitochondrial metabolism. Differential expression with respect to age suggested increased iron accumulation and decreased glycolysis, fatty acid metabolism, and proteasome function. Interaction of gene expression trends associated with viral infection and aging suggest that viral infection likely plays a role in aging in the Aplysia nervous system. Mutation analysis of viral RNA identified signatures suggesting ADAR and AID/APOBEC like deaminase act as part of Aplysia anti-viral defense.

Genomic Landscape and Regulation of RNA Editing in Pekin Ducks Susceptible to Duck Hepatitis A Virus Genotype 3 Infection

RNA editing is increasingly recognized as a post-transcriptional modification that directly affects viral infection by regulating RNA stability and recoding proteins. the duck hepatitis A virus genotype 3 (DHAV-3) infection is seriously detrimental to the Asian duck industry. However, the landscape and roles of RNA editing in the susceptibility and resistance of Pekin ducks to DHAV-3 remain unclear. Here, we profiled dynamic RNA editing events in liver tissue and investigated their potential functions during DHAV-3 infection in Pekin ducks. We identified 11,067 informative RNA editing sites in liver tissue from DHAV-3-susceptible and -resistant ducklings at three time points during virus infection. Differential RNA editing sites (DRESs) between S and R ducks were dynamically changed during infection, which were enriched in genes associated with vesicle-mediated transport and immune-related pathways. Moreover, we predicted and experimentally verified that RNA editing events in 3'-UTR could result in loss or gain of miRNA-mRNA interactions, thereby changing the expression of target genes. We also found a few DRESs in coding sequences (CDSs) that altered the amino acid sequences of several proteins that were vital for viral infection. Taken together, these data suggest that dynamic RNA editing has significant potential to tune physiological processes in response to virus infection in Pekin ducks, thus contributing to host differential susceptibility to DHAV-3.

Impact of Vaccination on Intra-Host Genetic Diversity of Patients Infected with SARS-CoV-2 Gamma Lineage

The high transmissibility, rapid evolution, and immune escape of SARS-CoV-2 variants can influence the course of infection and, in turn, morbidity and mortality in COVID-19, posing a challenge in controlling transmission rates and contributing to the emergence and spread of new variants. Understanding the factors that shape viral genetic variation is essential for comprehending the evolution and transmission of SARS-CoV-2, especially in vaccinated individuals where immune response plays a role in the progression and spread of this disease. In this context, we evaluated the impact of immunity induced by the CoronaVac vaccine (Butantan/Sinovac) on intra-host genetic diversity, analyzing 118 whole-genome sequences of SARS-CoV-2 from unvaccinated and vaccinated patients infected with the Gamma variant. Vaccination with CoronaVac favors negative selection at the intra-host level in different genomic regions. It prevents greater genetic diversity of SARS-CoV-2, reinforcing the importance of vaccination in reducing the emergence of new mutations and virus transmission.

RNA editing in response to COVID-19 vaccines: unveiling dynamic epigenetic regulation of host immunity

Background

COVID-19 vaccines are crucial for reducing the threat and burden of the pandemic on global public health, yet the epigenetic, especially RNA editing in response to the vaccines remains unelucidated.

Results

Our current study performed an epitranscriptomic analysis of RNA-Seq data of 260 blood samples from 102 healthy and SARS-CoV-2 naïve individuals receiving different doses of the COVID-19 vaccine and revealed dynamic, transcriptome-wide adenosine to inosine (A-to-I) RNA editing changes in response to COVID-19 vaccines (RNA editing in response to COVID-19 vaccines). 5592 differential RNA editing (DRE) sites in 1820 genes were identified, with most of them showing up-regulated RNA editing and correlated with increased expression of edited genes. These deferentially edited genes were primarily involved in immune- and virus-related gene functions and pathways. Differential ADAR expression probably contributed to RNA editing in response to COVID-19 vaccines. One of the most significant DRE in RNA editing in response to COVID-19 vaccines was in apolipoprotein L6 (APOL6) 3' UTR, which positively correlated with its up-regulated expression. In addition, recoded key antiviral and immune-related proteins such as IFI30 and GBP1 recoded by missense editing was observed as an essential component of RNA editing in response to COVID-19 vaccines. Furthermore, both RNA editing in response to COVID-19 vaccines and its functions dynamically depended on the number of vaccine doses.

Conclusion

Our results thus underscored the potential impact of blood RNA editing in response to COVID-19 vaccines on the host's molecular immune system.

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.

Refining SARS-CoV-2 intra-host variation by leveraging large-scale sequencing data

Understanding viral genome evolution during host infection is crucial for grasping viral diversity and evolution. Analyzing intra-host single nucleotide variants (iSNVs) offers insights into new lineage emergence, which is important for predicting and mitigating future viral threats. Despite next-generation sequencing's potential, challenges persist, notably sequencing artifacts leading to false iSNVs. We developed a workflow to enhance iSNV detection in large NGS libraries, using over 130 000 SARS-CoV-2 libraries to distinguish mutations from errors. Our approach integrates bioinformatics protocols, stringent quality control, and dimensionality reduction to tackle batch effects and improve mutation detection reliability. Additionally, we pioneer the application of the PHATE visualization approach to genomic data and introduce a methodology that quantifies how related groups of data points are represented within a two-dimensional space, enhancing clustering structure explanation based on genetic similarities. This workflow advances accurate intra-host mutation detection, facilitating a deeper understanding of viral diversity and evolution.

RNA editing regulates host immune response and T cell homeostasis in SARS-CoV-2 infection

Adenosine to inosine (A-to-I) RNA editing by ADAR1 has been implicated in maintaining self-tolerance, preventing autoimmunity, and mediating antiviral immunity. Foreign viral double-stranded RNA triggers rapid interferon response and activates ADAR1 in the host immune system. Emerging data points to a role of ADAR1 A-to-I editing in the inflammatory response associated with severe COVID-19 disease. We identify A-to-I editing events within human whole transcriptome data from SARS-CoV-2 infected individuals, non-infected individuals, and individuals with other viral illnesses from nasopharyngeal swabs. High levels of RNA editing in host cells are associated with low SARS-CoV-2 viral load (p = 9.27 E-06), suggesting an inhibitory effect of ADAR1 on viral infection. Additionally, we find differentially expressed genes associated with RNA-modifications and interferon response. Single cell RNA-sequencing analysis of SARS-CoV-2 infected nasopharyngeal swabs reveals that cytotoxic CD8 T cells upregulate ADAR1 in COVID-19 positive samples (p = 0.0269). We further reveal ADAR1 expression increases with CD4 and CD8 T cell activation, and knockdown of ADAR1 leads to apoptosis and aberrant IL-2 secretion. Together, our data suggests A-to-I RNA editing is required to maintain healthy homeostasis of activated T cells to combat SARS-CoV-2 infection.

Studying the intersection of nucleoside modifications and SARS-CoV-2 RNA-dependent RNA transcription using an in vitro reconstituted system

There is growing recognition that viral RNA genomes possess enzymatically incorporated modified nucleosides. These small chemical changes are analogous to epigenomic modifications in DNA and have the potential to be similarly important modulators of viral transcription and evolution. However, the molecular level consequences of individual sites of modification remain to be broadly explored. Here we describe an in vitro assay to examine the impact of nucleoside modifications on the rate and fidelity of SARS-CoV-2 RNA transcription. Establishing the role of modified nucleotides in SARS-CoV-2 is of interest both for advancing fundamental knowledge of RNA modifications in viruses, and because modulating the modification-landscape of SARS-CoV-2 may represent a therapeutic strategy to interfere with viral RNA replication. Our approach can be used to assess the influence both of modifications present in a template RNA, as well nucleotide analog inhibitors. These methods provide a reproducible guide for generating active SARS-CoV-2 replication/transcription complexes capable of establishing how RNA modifications influence the pre-steady state rate constants of nucleotide addition by RNA-dependent RNA polymerases.

Adenosine modifications impede SARS-CoV-2 RNA-dependent RNA transcription

SARS-CoV-2, the causative virus of the COVID-19 pandemic, follows SARS and MERS as recent zoonotic coronaviruses causing severe respiratory illness and death in humans. The recurrent impact of zoonotic coronaviruses demands a better understanding of their fundamental molecular biochemistry. Nucleoside modifications, which modulate many steps of the RNA life cycle, have been found in SARS-CoV-2 RNA, although whether they confer a pro- or antiviral effect is unknown. Regardless, the viral RNA-dependent RNA polymerase will encounter these modifications as it transcribes through the viral genomic RNA. We investigated the functional consequences of nucleoside modification on the pre-steady state kinetics of SARS-CoV-2 RNA-dependent RNA transcription using an in vitro reconstituted transcription system with modified RNA templates. Our findings show that N 6-methyladenosine and 2'-O-methyladenosine modifications slow the rate of viral transcription at magnitudes specific to each modification, which has the potential to impact SARS-CoV-2 genome maintenance.

Potential Role of APOBEC3 Family Proteins in SARS-CoV-2 Replication

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has acquired multiple mutations since its emergence. Analyses of the SARS-CoV-2 genomes from infected patients exhibit a bias toward C-to-U mutations, which are suggested to be caused by the apolipoprotein B mRNA editing enzyme polypeptide-like 3 (APOBEC3, A3) cytosine deaminase proteins. However, the role of A3 enzymes in SARS-CoV-2 replication remains unclear. To address this question, we investigated the effect of A3 family proteins on SARS-CoV-2 replication in the myeloid leukemia cell line THP-1 lacking A3A to A3G genes. The Wuhan, BA.1, and BA.5 variants had comparable viral replication in parent and A3A-to-A3G-null THP-1 cells stably expressing angiotensin-converting enzyme 2 (ACE2) protein. On the other hand, the replication and infectivity of these variants were abolished in A3A-to-A3G-null THP-1-ACE2 cells in a series of passage experiments over 20 days. In contrast to previous reports, we observed no evidence of A3-induced SARS-CoV-2 mutagenesis in the passage experiments. Furthermore, our analysis of a large number of publicly available SARS-CoV-2 genomes did not reveal conclusive evidence for A3-induced mutagenesis. Our studies suggest that A3 family proteins can positively contribute to SARS-CoV-2 replication; however, this effect is deaminase-independent.

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.

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

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

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.

Proteolytic cleavage and inactivation of the TRMT1 tRNA modification enzyme by SARS-CoV-2 main protease

Nonstructural protein 5 (Nsp5) is the main protease of SARS-CoV-2 that cleaves viral polyproteins into individual polypeptides necessary for viral replication. Here, we show that Nsp5 binds and cleaves human tRNA methyltransferase 1 (TRMT1), a host enzyme required for a prevalent post-transcriptional modification in tRNAs. Human cells infected with SARS-CoV-2 exhibit a decrease in TRMT1 protein levels and TRMT1-catalyzed tRNA modifications, consistent with TRMT1 cleavage and inactivation by Nsp5. Nsp5 cleaves TRMT1 at a specific position that matches the consensus sequence of SARS-CoV-2 polyprotein cleavage sites, and a single mutation within the sequence inhibits Nsp5-dependent proteolysis of TRMT1. The TRMT1 cleavage fragments exhibit altered RNA binding activity and are unable to rescue tRNA modification in TRMT1-deficient human cells. Compared to wild-type human cells, TRMT1-deficient human cells infected with SARS-CoV-2 exhibit reduced levels of intracellular viral RNA. These findings provide evidence that Nsp5-dependent cleavage of TRMT1 and perturbation of tRNA modification patterns contribute to the cellular pathogenesis of SARS-CoV-2 infection.

Targeted accurate RNA consensus sequencing (tARC-seq) reveals mechanisms of replication error affecting SARS-CoV-2 divergence

RNA viruses, like SARS-CoV-2, depend on their RNA-dependent RNA polymerases (RdRp) for replication, which is error prone. Monitoring replication errors is crucial for understanding the virus's evolution. Current methods lack the precision to detect rare de novo RNA mutations, particularly in low-input samples such as those from patients. Here we introduce a targeted accurate RNA consensus sequencing method (tARC-seq) to accurately determine the mutation frequency and types in SARS-CoV-2, both in cell culture and clinical samples. Our findings show an average of 2.68 × 10-5 de novo errors per cycle with a C > T bias that cannot be solely attributed to APOBEC editing. We identified hotspots and cold spots throughout the genome, correlating with high or low GC content, and pinpointed transcription regulatory sites as regions more susceptible to errors. tARC-seq captured template switching events including insertions, deletions and complex mutations. These insights shed light on the genetic diversity generation and evolutionary dynamics of SARS-CoV-2.

Ocular A-to-I RNA editing signatures associated with SARS-CoV-2 infection

Ophthalmic manifestations have recently been observed in acute and post-acute complications of COVID-19 caused by SARS-CoV-2 infection. Our precious study has shown that host RNA editing is linked to RNA viral infection, yet ocular adenosine to inosine (A-to-I) RNA editing during SARS-CoV-2 infection remains uninvestigated in COVID-19. Herein we used an epitranscriptomic pipeline to analyze 37 samples and investigate A-to-I editing associated with SARS-CoV-2 infection, in five ocular tissue types including the conjunctiva, limbus, cornea, sclera, and retinal organoids. Our results revealed dramatically altered A-to-I RNA editing across the five ocular tissues. Notably, the transcriptome-wide average level of RNA editing was increased in the cornea but generally decreased in the other four ocular tissues. Functional enrichment analysis showed that differential RNA editing (DRE) was mainly in genes related to ubiquitin-dependent protein catabolic process, transcriptional regulation, and RNA splicing. In addition to tissue-specific RNA editing found in each tissue, common RNA editing was observed across different tissues, especially in the innate antiviral immune gene MAVS and the E3 ubiquitin-protein ligase MDM2. Analysis in retinal organoids further revealed highly dynamic RNA editing alterations over time during SARS-CoV-2 infection. Our study thus suggested the potential role played by RNA editing in ophthalmic manifestations of COVID-19, and highlighted its potential transcriptome impact, especially on innate immunity.

Molnupiravir increases SARS-CoV-2 genome diversity and complexity: A case-control cohort study

Molnupiravir, an oral direct-acting antiviral effective in vitro against SARS-CoV-2, has been largely employed during the COVID-19 pandemic, since December 2021. After marketing and widespread usage, a progressive increase in SARS-CoV-2 lineages characterized by a higher transition/transversion ratio, a characteristic signature of molnupiravir action, appeared in the Global Initiative on Sharing All Influenza Data (GISAID) and International Nucleotide Sequence Database Collaboration (INSDC) databases. Here, we assessed the drug effects by SARS-CoV-2 whole-genome sequencing on 38 molnupiravir-treated persistently positive COVID-19 outpatients tested before and after treatment. Seventeen tixagevimab/cilgavimab-treated outpatients served as controls. Mutational analyses confirmed that SARS-CoV-2 exhibits an increased transition/transversion ratio seven days after initiation of molnupiravir. Moreover we observed an increased G->A ratio compared to controls, which was not related to apolipoprotein B mRNAediting enzyme, catalytic polypeptide-like (APOBEC) activity. In addition, we demonstrated for the first time an increased diversity and complexity of the viral quasispecies.

High frequency of transition to transversion ratio in the stem region of RNA secondary structure of untranslated region of SARS-CoV-2

Introduction

The propensity of nucleotide bases to form pairs, causes folding and the formation of secondary structure in the RNA. Therefore, purine (R): pyrimidine (Y) base-pairing is vital to maintain uniform lateral dimension in RNA secondary structure. Transversions or base substitutions between R and Y bases, are more detrimental to the stability of RNA secondary structure, than transitions derived from substitutions between A and G or C and T. The study of transversion and transition base substitutions is important to understand evolutionary mechanisms of RNA secondary structure in the 5'  and 3'  untranslated (UTR) regions of SARS-CoV-2. In this work, we carried out comparative analysis of transition and transversion base substitutions in the stem and loop regions of RNA secondary structure of SARS-CoV-2.

Methods

We have considered the experimentally determined and well documented stem and loop regions of 5' and 3' UTR regions of SARS-CoV-2 for base substitution analysis. The secondary structure comprising of stem and loop regions were visualized using the RNAfold web server. The GISAID repository was used to extract base sequence alignment of the UTR regions. Python scripts were developed for comparative analysis of transversion and transition frequencies in the stem and the loop regions.

Results

The results of base substitution analysis revealed a higher transition (ti) to transversion (tv) ratio (ti/tv) in the stem region of UTR of RNA secondary structure of SARS-CoV-2 reported during the early stage of the pandemic. The higher ti/tv ratio in the stem region suggested the influence of secondary structure in selecting the pattern of base substitutions. This differential pattern of ti/tv values between stem and loop regions was not observed among the Delta and Omicron variants that dominated the later stage of the pandemic. It is noteworthy that the ti/tv values in the stem and loop regions were similar among the later dominant Delta and Omicron variant strains which is to be investigated to understand the rapid evolution and global adaptation of SARS-CoV-2.

Conclusion

Our findings implicate the lower frequency of transversions than the transitions in the stem regions of UTRs of SARS-CoV-2. The RNA secondary structures are associated with replication, translation, and packaging, further investigations are needed to understand these base substitutions across different variants of SARS-CoV-2.

PlantC2U: deep learning of cross-species sequence landscapes predicts plastid C-to-U RNA editing in plants

In plants, C-to-U RNA editing mainly occurs in plastid and mitochondrial transcripts, which contributes to a complex transcriptional regulatory network. More evidence reveals that RNA editing plays critical roles in plant growth and development. However, accurate detection of RNA editing sites using transcriptome sequencing data alone is still challenging. In the present study, we develop PlantC2U, which is a convolutional neural network, to predict plastid C-to-U RNA editing based on the genomic sequence. PlantC2U achieves >95% sensitivity and 99% specificity, which outperforms the PREPACT tool, random forests, and support vector machines. PlantC2U not only further checks RNA editing sites from transcriptome data to reduce possible false positives, but also assesses the effect of different mutations on C-to-U RNA editing based on the flanking sequences. Moreover, we found the patterns of tissue-specific RNA editing in the mangrove plant Kandelia obovata, and observed reduced C-to-U RNA editing rates in the cold stress response of K. obovata, suggesting their potential regulatory roles in plant stress adaptation. In addition, we present RNAeditDB, available online at https://jasonxu.shinyapps.io/RNAeditDB/. Together, PlantC2U and RNAeditDB will help researchers explore the RNA editing events in plants and thus will be of broad utility for the plant research community.

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.

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.

Proteolytic cleavage and inactivation of the TRMT1 tRNA modification enzyme by SARS-CoV-2 main protease

Nonstructural protein 5 (Nsp5) is the main protease of SARS-CoV-2 that cleaves viral polyproteins into individual polypeptides necessary for viral replication. Here, we show that Nsp5 binds and cleaves human tRNA methyltransferase 1 (TRMT1), a host enzyme required for a prevalent post-transcriptional modification in tRNAs. Human cells infected with SARS-CoV-2 exhibit a decrease in TRMT1 protein levels and TRMT1-catalyzed tRNA modifications, consistent with TRMT1 cleavage and inactivation by Nsp5. Nsp5 cleaves TRMT1 at a specific position that matches the consensus sequence of SARS-CoV-2 polyprotein cleavage sites, and a single mutation within the sequence inhibits Nsp5-dependent proteolysis of TRMT1. The TRMT1 cleavage fragments exhibit altered RNA binding activity and are unable to rescue tRNA modification in TRMT1-deficient human cells. Compared to wildtype human cells, TRMT1-deficient human cells infected with SARS-CoV-2 exhibit reduced levels of intracellular viral RNA. These findings provide evidence that Nsp5-dependent cleavage of TRMT1 and perturbation of tRNA modification patterns contribute to the cellular pathogenesis of SARS-CoV-2 infection.

Liver cancer-specific mutations in functional domains of ADAR2 lead to the elevation of coding and non-coding RNA editing in multiple tumor-related genes

Mutation is the major cause of phenotypic innovations. Apart from DNA mutations, the alteration on RNA such as the ADAR-mediated A-to-I RNA editing could also shape the phenotype. These two layers of variations have not been systematically combined to study their collective roles in cancers. We collected the high-quality transcriptomes of ten hepatocellular carcinoma (HCC) and the matched control samples. We systematically identified HCC-specific mutations in the exonic regions and profiled the A-to-I RNA editome in each sample. All ten HCC samples had mutations in the CDS of ADAR2 gene (dsRNA-binding domain or catalytic domain). The consequence of these mutations converged to the elevation of ADAR2 efficiency as reflected by the global increase of RNA editing levels in HCC. The up-regulated editing sites (UES) were enriched in the CDS and UTR of oncogenes and tumor suppressor genes (TSG), indicating the possible roles of these target genes in HCC oncogenesis. We present the mutation-ADAR2-UES-oncogene/TSG-HCC axis that explains how mutations at different layers would finally lead to abnormal phenotype. In the light of central dogma, our work provides novel insights into how to fully take advantage of the transcriptome data to decipher the consequence of mutations.

In silico and experimental approaches for validating RNA editing events in transcriptomes

As a typical RNA virus, SARS-CoV-2 is subjected to RNA editing in host cells. While some researchers believe that a traditional variant calling pipeline retrieves all true-positive RNA editing events from the transcriptome, others argue that conventional methods identify many false-positive sites. Here, I describe several additional in silico and experimental approaches to validate the authenticity of RNA editing in SARS-CoV-2. These approaches include requiring strand-specific sequencing, analysis of hyperedited reads, linkage analysis, orthogonal methods like mass spectrometry, and the use of ADAR-deficient host cells. These findings may improve future analyses on the identification of RNA editing, especially in RNA viruses.

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

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

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

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

An orthology-based methodology as a complementary approach to retrieve evolutionarily conserved A-to-I RNA editing sites

Adar-mediated adenosine-to-inosine (A-to-I) mRNA editing is a conserved mechanism that exerts diverse regulatory functions during the development, evolution, and adaptation of metazoans. The accurate detection of RNA editing sites helps us understand their biological significance. In this work, with an improved genome assembly of honeybee (Apis mellifera), we used a new orthology-based methodology to complement the traditional pipeline of (de novo) RNA editing detection. Compared to the outcome of traditional pipeline, we retrieved many novel editing sites in CDS that are deeply conserved between honeybee and other distantly related insects. The newly retrieved sites were missed by the traditional de novo identification due to the stringent criteria for controlling false-positive rate. Caste-specific editing sites are identified, including an Ile>Met auto-recoding site in Adar. This recoding was even conserved between honeybee and bumblebee, suggesting its putative regulatory role in shaping the phenotypic plasticity of eusocial Hymenoptera. In summary, we proposed a complementary approach to the traditional pipeline and retrieved several previously unnoticed CDS editing sites. From both technical and biological aspects, our works facilitate future researches on finding the functional editing sites and advance our understanding on the connection between RNA editing and the great phenotypic diversity of organisms.

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

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

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

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

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.