ADARs: viruses and innate immunity

Prediction of Prospective Mutational Landscape of SARS-CoV-2 Spike ssRNA and Evolutionary Basis of Its Host Interaction

Booster doses are crucial against severe COVID-19, as rapid virus mutations and variant emergence prolong the pandemic crisis. The virus's quick evolution, short generation-time, and adaptive changes impact virulence and evolvability, helping predictions about variant of concerns' (VOCs') landscapes. Here, in this study, we used a new computational algorithm, to predict the mutational pattern in SARS-CoV-2 ssRNA, proteomics, structural identification, mutation stability, and functional correlation, as well as immune escape mechanisms. Interestingly, the sequence diversity of SARS Coronavirus-2 has demonstrated a predominance of G- > A and C- > U substitutions. The best validation statistics are explored here in seven homologous models of the expected mutant SARS-CoV-2 spike ssRNA and employed for hACE2 and IgG interactions. The interactome profile of SARS-CoV-2 spike with hACE2 and IgG revealed a strong correlation between phylogeny and divergence time. The systematic adaptation of SARS-CoV-2 spike ssRNA influences infectivity and immune escape. Data suggest higher propensity of Adenine rich sequence promotes MHC system avoidance, preferred by A-rich codons. Phylogenetic data revealed the evolution of SARS-CoV-2 lineages' epidemiology. Our findings may unveil processes governing the genesis of immune-resistant variants, prompting a critical reassessment of the coronavirus mutation rate and exploration of hypotheses beyond mechanical aspects.

Inosine-Induced Base Pairing Diversity during Reverse Transcription

A-to-I editing catalyzed by adenosine deaminase acting on RNAs impacts numerous physiological and biochemical processes that are essential for cellular functions and is a big contributor to the infectivity of certain RNA viruses. The outcome of this deamination leads to changes in the eukaryotic transcriptome functionally resembling A-G transitions since inosine preferentially pairs with cytosine. Moreover, hyper-editing or multiple A to G transitions in clusters were detected in measles virus. Inosine modifications either directly on viral RNA or on cellular RNA can have antiviral or pro-viral repercussions. While many of the significant roles of inosine in cellular RNAs are well understood, the effects of hyper-editing of A to I on viral polymerase activity during RNA replication remain elusive. Moreover, biological strategies such as molecular cloning and RNA-seq for transcriptomic interrogation rely on RT-polymerase chain reaction with little to no emphasis placed on the first step, reverse transcription, which may reshape the sequencing results when hypermodification is present. In this study, we systematically explore the influence of inosine modification, varying the number and position of inosines, on decoding outcomes using three different reverse transcriptases (RTs) followed by standard Sanger sequencing. We find that inosine alone or in clusters can differentially affect the RT activity. To gain structural insights into the accommodation of inosine in the polymerase site of HIV-1 reverse transcriptase (HIV-1-RT) and how this structural context affects the base pairing rules for inosine, we performed molecular dynamics simulations of the HIV-1-RT. The simulations highlight the importance of the protein-nucleotide interaction as a critical factor in deciphering the base pairing behavior of inosine clusters. This effort sets the groundwork for decrypting the physiological significance of inosine and linking the fidelity of reverse transcriptase and the possible diverse transcription outcomes of cellular RNAs and/or viral RNAs where hyper-edited inosines are present in the transcripts.

The epitranscriptome of high-grade gliomas: a promising therapeutic target with implications from the tumor microenvironment to endogenous retroviruses

Glioblastoma (GBM) comprises 45.6% of all primary malignant brain cancers and is one of the most common and aggressive intracranial tumors in adults. Intratumoral heterogeneity with a wide range of proteomic, genetic, and epigenetic dysregulation contributes to treatment resistance and poor prognosis, thus demanding novel therapeutic approaches. To date, numerous clinical trials have been developed to target the proteome and epigenome of high-grade gliomas with promising results. However, studying RNA modifications, or RNA epitranscriptomics, is a new frontier within neuro-oncology. RNA epitranscriptomics was discovered in the 1970s, but in the last decade, the extent of modification of mRNA and various non-coding RNAs has emerged and been implicated in transposable element activation and many other oncogenic processes within the tumor microenvironment. This review provides background information and discusses the therapeutic potential of agents modulating epitranscriptomics in high-grade gliomas. A particular emphasis will be placed on how combination therapies that include immune agents targeting hERV-mediated viral mimicry could improve the treatment of GBM.

Hepatitis D infection induces IFN-β-mediated NK cell activation and TRAIL-dependent cytotoxicity

Background and aims

The co-infection of hepatitis B (HBV) patients with the hepatitis D virus (HDV) causes the most severe form of viral hepatitis and thus drastically worsens the course of the disease. Therapy options for HBV/HDV patients are still limited. Here, we investigated the potential of natural killer (NK) cells that are crucial drivers of the innate immune response against viruses to target HDV-infected hepatocytes.

Methods

We established in vitro co-culture models using HDV-infected hepatoma cell lines and human peripheral blood NK cells. We determined NK cell activation by flow cytometry, transcriptome analysis, bead-based cytokine immunoassays, and NK cell-mediated effects on T cells by flow cytometry. We validated the mechanisms using CRISPR/Cas9-mediated gene deletions. Moreover, we assessed the frequencies and phenotype of NK cells in peripheral blood of HBV and HDV superinfected patients.

Results

Upon co-culture with HDV-infected hepatic cell lines, NK cells upregulated activation markers, interferon-stimulated genes (ISGs) including the death receptor ligand tumor necrosis factor-related apoptosis-inducing ligand (TRAIL), produced interferon (IFN)-γ and eliminated HDV-infected cells via the TRAIL-TRAIL-R2 axis. We identified IFN-β released by HDV-infected cells as an important enhancer of NK cell activity. In line with our in vitro data, we observed activation of peripheral blood NK cells from HBV/HDV co-infected, but not HBV mono-infected patients.

Conclusion

Our data demonstrate NK cell activation in HDV infection and their potential to eliminate HDV-infected hepatoma cells via the TRAIL/TRAIL-R2 axis which implies a high relevance of NK cells for the design of novel anti-viral therapies.

The regulation of antiviral innate immunity through non-m6A RNA modifications

The post-transcriptional RNA modifications impact the dynamic regulation of gene expression in diverse biological and physiological processes. Host RNA modifications play an indispensable role in regulating innate immune responses against virus infection in mammals. Meanwhile, the viral RNAs can be deposited with RNA modifications to interfere with the host immune responses. The N6-methyladenosine (m6A) has boosted the recent emergence of RNA epigenetics, due to its high abundance and a transcriptome-wide widespread distribution in mammalian cells, proven to impact antiviral innate immunity. However, the other types of RNA modifications are also involved in regulating antiviral responses, and the functional roles of these non-m6A RNA modifications have not been comprehensively summarized. In this Review, we conclude the regulatory roles of 2'-O-methylation (Nm), 5-methylcytidine (m5C), adenosine-inosine editing (A-to-I editing), pseudouridine (Ψ), N1-methyladenosine (m1A), N7-methylguanosine (m7G), N6,2'-O-dimethyladenosine (m6Am), and N4-acetylcytidine (ac4C) in antiviral innate immunity. We provide a systematic introduction to the biogenesis and functions of these non-m6A RNA modifications in viral RNA, host RNA, and during virus-host interactions, emphasizing the biological functions of RNA modification regulators in antiviral responses. Furthermore, we discussed the recent research progress in the development of antiviral drugs through non-m6A RNA modifications. Collectively, this Review conveys knowledge and inspiration to researchers in multiple disciplines, highlighting the challenges and future directions in RNA epitranscriptome, immunology, and virology.

Adenosine deaminase acting on RNA 1 (ADAR1) as crucial regulators in cardiovascular diseases: structures, pathogenesis, and potential therapeutic approach

Cardiovascular diseases (CVDs) are a group of diseases that have a major impact on global health and are the leading cause of death. A large number of chemical base modifications in ribonucleic acid (RNA) are associated with cardiovascular diseases. A variety of ribonucleic acid modifications exist in cells, among which adenosine deaminase-dependent modification is one of the most common ribonucleic acid modifications. Adenosine deaminase acting on ribonucleic acid 1 (Adenosine deaminase acting on RNA 1) is a widely expressed double-stranded ribonucleic acid adenosine deaminase that forms inosine (A-to-I) by catalyzing the deamination of adenosine at specific sites of the target ribonucleic acid. In this review, we provide a comprehensive overview of the structure of Adenosine deaminase acting on RNA 1 and summarize the regulatory mechanisms of ADAR1-mediated ribonucleic acid editing in cardiovascular diseases, indicating Adenosine deaminase acting on RNA 1 as a promising therapeutic target in cardiovascular diseases.

SARS-CoV-2-associated organs failure and inflammation: a focus on the role of cellular and viral microRNAs

SARS-CoV-2 has been responsible for the recent pandemic all over the world, which has caused many complications. One of the hallmarks of SARS-CoV-2 infection is an induced immune dysregulation, in some cases resulting in cytokine storm syndrome, acute respiratory distress syndrome and many organs such as lungs, brain, and heart that are affected during the SARS-CoV-2 infection. Several physiological parameters are altered as a result of infection and cytokine storm. Among them, microRNAs (miRNAs) might reflect this poor condition since they play a significant role in immune cellular performance including inflammatory responses. Both host and viral-encoded miRNAs are crucial for the successful infection of SARS-CoV-2. For instance, dysregulation of miRNAs that modulate multiple genes expressed in COVID-19 patients with comorbidities (e.g., type 2 diabetes, and cerebrovascular disorders) could affect the severity of the disease. Therefore, altered expression levels of circulating miRNAs might be helpful to diagnose this illness and forecast whether a COVID-19 patient could develop a severe state of the disease. Moreover, a number of miRNAs could inhibit the expression of proteins, such as ACE2, TMPRSS2, spike, and Nsp12, involved in the life cycle of SARS-CoV-2. Accordingly, miRNAs represent potential biomarkers and therapeutic targets for this devastating viral disease. In the current study, we investigated modifications in miRNA expression and their influence on COVID-19 disease recovery, which may be employed as a therapy strategy to minimize COVID-19-related disorders.

Enzymes in the time of COVID-19: An overview about the effects in the human body, enzyme market, and perspectives for new drugs

The rising pandemic caused by a coronavirus, resulted in a scientific quest to discover some effective treatments against its etiologic agent, the severe acute respiratory syndrome‐coronavirus 2 (SARS‐CoV‐2). This research represented a significant scientific landmark and resulted in many medical advances. However, efforts to understand the viral mechanism of action and how the human body machinery is subverted during the infection are still ongoing. Herein, we contributed to this field with this compilation of the roles of both viral and human enzymes in the context of SARS‐CoV‐2 infection. In this sense, this overview reports that proteases are vital for the infection to take place: from SARS‐CoV‐2 perspective, the main protease (M^pro) and papain‐like protease (PL^pro) are highlighted; from the human body, angiotensin‐converting enzyme‐2, transmembrane serine protease‐2, and cathepsins (CatB/L) are pointed out. In addition, the influence of the virus on other enzymes is reported as the JAK/STAT pathway and the levels of lipase, enzymes from the cholesterol metabolism pathway, amylase, aspartate aminotransferase, alanine aminotransferase, lactate dehydrogenase, and glyceraldehyde 3‐phosphate dehydrogenase are also be disturbed in SARS‐CoV‐2 infection. Finally, this paper discusses the importance of detailed enzymatic studies for future treatments against SARS‐CoV‐2, and how some issues related to the syndrome treatment can create opportunities in the biotechnological market of enzymes and the development of new drugs.

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

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

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

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

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

Detection of A-to-I RNA Editing in SARS-COV-2

ADAR1-mediated deamination of adenosines in long double-stranded RNAs plays an important role in modulating the innate immune response. However, recent investigations based on metatranscriptomic samples of COVID-19 patients and SARS-COV-2-infected Vero cells have recovered contrasting findings. Using RNAseq data from time course experiments of infected human cell lines and transcriptome data from Vero cells and clinical samples, we prove that A-to-G changes observed in SARS-COV-2 genomes represent genuine RNA editing events, likely mediated by ADAR1. While the A-to-I editing rate is generally low, changes are distributed along the entire viral genome, are overrepresented in exonic regions, and are (in the majority of cases) nonsynonymous. The impact of RNA editing on virus-host interactions could be relevant to identify potential targets for therapeutic interventions.

Host-dependent editing of SARS-CoV-2 in COVID-19 patients

A common trait among RNA viruses is their high capability to acquire genetic variability due to viral and host mechanisms. Next-generation sequencing (NGS) analysis enables the deep study of the viral quasispecies in samples from infected individuals. In this study, the viral quasispecies complexity and single nucleotide polymorphisms of the SARS-CoV-2 spike gene of coronavirus disease 2019 (COVID-19) patients with mild or severe disease were investigated using next-generation sequencing (Illumina platform). SARS-CoV-2 spike variability was higher in patients with long-lasting infection. Most substitutions found were present at frequencies lower than 1%, and had an A → G or T → C pattern, consistent with variants caused by adenosine deaminase acting on RNA-1 (ADAR1). ADAR1 affected a small fraction of replicating genomes, but produced multiple, mainly non-synonymous mutations. ADAR1 editing during replication rather than the RNA-dependent RNA polymerase (nsp12) was the predominant mechanism generating SARS-CoV-2 genetic variability. However, the mutations produced are not fixed in the infected human population, suggesting that ADAR1 may have an antiviral role, whereas nsp12-induced mutations occurring in patients with high viremia and persistent infection are the main source of new SARS-CoV-2 variants.

ADAR Editing in Viruses: An Evolutionary Force to Reckon with

Adenosine Deaminases that Act on RNA (ADARs) are RNA editing enzymes that play a dynamic and nuanced role in regulating transcriptome and proteome diversity. This editing can be highly selective, affecting a specific site within a transcript, or nonselective, resulting in hyperediting. ADAR editing is important for regulating neural functions and autoimmunity, and has a key role in the innate immune response to viral infections, where editing can have a range of pro- or antiviral effects and can contribute to viral evolution. Here we examine the role of ADAR editing across a broad range of viral groups. We propose that the effect of ADAR editing on viral replication, whether pro- or antiviral, is better viewed as an axis rather than a binary, and that the specific position of a given virus on this axis is highly dependent on virus- and host-specific factors, and can change over the course of infection. However, more research needs to be devoted to understanding these dynamic factors and how they affect virus-ADAR interactions and viral evolution. Another area that warrants significant attention is the effect of virus-ADAR interactions on host-ADAR interactions, particularly in light of the crucial role of ADAR in regulating neural functions. Answering these questions will be essential to developing our understanding of the relationship between ADAR editing and viral infection. In turn, this will further our understanding of the effects of viruses such as SARS-CoV-2, as well as many others, and thereby influence our approach to treating these deadly diseases.

Strategies for controlling the innate immune activity of conventional and self-amplifying mRNA therapeutics: Getting the message across

The recent approval of messenger RNA (mRNA)-based vaccines to combat the SARS-CoV-2 pandemic highlights the potential of both conventional mRNA and self-amplifying mRNA (saRNA) as a flexible immunotherapy platform to treat infectious diseases. Besides the antigen it encodes, mRNA itself has an immune-stimulating activity that can contribute to vaccine efficacy. This self-adjuvant effect, however, will interfere with mRNA translation and may influence the desired therapeutic outcome. To further exploit its potential as a versatile therapeutic platform, it will be crucial to control mRNA's innate immune-stimulating properties. In this regard, we describe the mechanisms behind the innate immune recognition of mRNA and provide an extensive overview of strategies to control its innate immune-stimulating activity. These strategies range from modifications to the mRNA backbone itself, optimization of production and purification processes to the combination with innate immune inhibitors. Furthermore, we discuss the delicate balance of the self-adjuvant effect in mRNA vaccination strategies, which can be both beneficial and detrimental to the therapeutic outcome.

3' untranslated regions of Marburg and Ebola virus mRNAs possess negative regulators of translation that are modulated by ADAR1 editing

The filovirus family includes deadly pathogens such as Ebola virus (EBOV) and Marburg virus (MARV). A substantial portion of filovirus genomes encode 5′ and 3′ untranslated regions (UTRs) of viral mRNAs. Select viral genomic RNA sequences corresponding to 3′ UTRs are prone to editing by adenosine deaminase acting on RNA 1 (ADAR1). A reporter mRNA approach, in which different 5′ or 3′ UTRs were inserted into luciferase-encoding mRNAs, demonstrates that MARV 3′ UTRs yield different levels of reporter gene expression, suggesting modulation of translation. The modulation occurs in cells unable to produce microRNAs (miRNAs) and can be recapitulated in a MARV minigenome assay. Deletion mutants identified negative regulatory regions at the ends of the MARV nucleoprotein (NP) and large protein (L) 3′ UTRs. Apparent ADAR1 editing mutants were previously identified within the MARV NP 3′ UTR. Introduction of these changes into the MARV nucleoprotein (NP) 3′ UTR or deletion of the region targeted for editing enhances translation, as indicated by reporter assays and polysome analysis. In addition, the parental NP 3′ UTR, but not the edited or deletion mutant NP 3′ UTRs, induces a type I interferon (IFN) response upon transfection into cells. Because some EBOV isolates from the West Africa outbreak exhibited ADAR1 editing of the viral protein of 40 kDa (VP40) 3′ UTR, VP40 3′ UTRs with parental and edited sequences were similarly assayed. The EBOV VP40 3′ UTR edits also enhanced translation, but neither the wild-type nor the edited 3′ UTRs induced IFN. These findings implicate filoviral mRNA 3′ UTRs as negative regulators of translation that can be inactivated by innate immune responses that induce ADAR1.

IMPORTANCE UTRs comprise a large percentage of filovirus genomes and are apparent targets of editing by ADAR1, an enzyme with pro- and antiviral activities. However, the functional significance of the UTRs and ADAR1 editing has been uncertain. This study demonstrates that MARV and EBOV 3′ UTRs can modulate translation, in some cases negatively. ADAR1 editing or deletion of select regions within the translation suppressing 3′ UTRs relieves the negative effects of the UTRs. These data indicate that filovirus 3′ UTRs contain translation regulatory elements that are modulated by activation of ADAR1, suggesting a complex interplay between filovirus gene expression and innate immunity.

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

The rapid evolution of RNA viruses has been long considered to result from a combination of high copying error frequencies during RNA replication, short generation times and the consequent extensive fixation of neutral or adaptive changes over short periods. While both the identities and sites of mutations are typically modelled as being random, recent investigations of sequence diversity of SARS coronavirus 2 (SARS-CoV-2) have identified a preponderance of C->U transitions, proposed to be driven by an APOBEC-like RNA editing process. The current study investigated whether this phenomenon could be observed in datasets of other RNA viruses. Using a 5% divergence filter to infer directionality, 18 from 36 datasets of aligned coding region sequences from a diverse range of mammalian RNA viruses (including Picornaviridae, Flaviviridae, Matonaviridae, Caliciviridae and Coronaviridae) showed a >2-fold base composition normalised excess of C->U transitions compared to U->C (range 2.1x-7.5x), with a consistently observed favoured 5' U upstream context. The presence of genome scale RNA secondary structure (GORS) was the only other genomic or structural parameter significantly associated with C->U/U->C transition asymmetries by multivariable analysis (ANOVA), potentially reflecting RNA structure dependence of sites targeted for C->U mutations. Using the association index metric, C->U changes were specifically over-represented at phylogenetically uninformative sites, potentially paralleling extensive homoplasy of this transition reported in SARS-CoV-2. Although mechanisms remain to be functionally characterised, excess C->U substitutions accounted for 11-14% of standing sequence variability of structured viruses and may therefore represent a potent driver of their sequence diversification and longer-term evolution.

Evolutionary Comparative Analyses of DNA-Editing Enzymes of the Immune System: From 5-Dimensional Description of Protein Structures to Immunological Insights and Applications to Protein Engineering

The immune system is unique among all biological sub-systems in its usage of DNA-editing enzymes to introduce targeted gene mutations and double-strand DNA breaks to diversify antigen receptor genes and combat viral infections. These processes, initiated by specific DNA-editing enzymes, often result in mistargeted induction of genome lesions that initiate and drive cancers. Like other molecules involved in human health and disease, the DNA-editing enzymes of the immune system have been intensively studied in humans and mice, with little attention paid (< 1% of published studies) to the same enzymes in evolutionarily distant species. Here, we present a systematic review of the literature on the characterization of one such DNA-editing enzyme, activation-induced cytidine deaminase (AID), from an evolutionary comparative perspective. The central thesis of this review is that although the evolutionary comparative approach represents a minuscule fraction of published works on this and other DNA-editing enzymes, this approach has made significant impacts across the fields of structural biology, immunology, and cancer research. Using AID as an example, we highlight the value of the evolutionary comparative approach in discoveries already made, and in the context of emerging directions in immunology and protein engineering. We introduce the concept of 5-dimensional (5D) description of protein structures, a more nuanced view of a structure that is made possible by evolutionary comparative studies. In this higher dimensional view of a protein's structure, the classical 3-dimensional (3D) structure is integrated in the context of real-time conformations and evolutionary time shifts (4^th dimension) and the relevance of these dynamics to its biological function (5^th dimension).

The Evolution of Severe Acute Respiratory Syndrome Coronavirus-2 during Pandemic and Adaptation to the Host

Severe Acute Respiratory Syndrome Coronavirus-2 is a zoonotic virus with a possible origin in bats and potential transmission to humans through an intermediate host. When zoonotic viruses jump to a new host, they undergo both mutational and natural selective pressures that result in non-synonymous and synonymous adaptive changes, necessary for efficient replication and rapid spread of diseases in new host species. The nucleotide composition and codon usage pattern of SARS-CoV-2 indicate the presence of a highly conserved, gene-specific codon usage bias. The codon usage pattern of SARS-CoV-2 is mostly antagonistic to human and bat codon usage. SARS-CoV-2 codon usage bias is mainly shaped by the natural selection, while mutational pressure plays a minor role. The time-series analysis of SARS-CoV-2 genome indicates that the virus is slowly evolving. Virus isolates from later stages of the outbreak have more biased codon usage and nucleotide composition than virus isolates from early stages of the outbreak.

A mark of disease: how mRNA modifications shape genetic and acquired pathologies

RNA modifications have recently emerged as a widespread and complex facet of gene expression regulation. Counting more than 170 distinct chemical modifications with far-reaching implications for RNA fate, they are collectively referred to as the epitranscriptome. These modifications can occur in all RNA species, including messenger RNAs (mRNAs) and noncoding RNAs (ncRNAs). In mRNAs the deposition, removal, and recognition of chemical marks by writers, erasers and readers influence their structure, localization, stability, and translation. In turn, this modulates key molecular and cellular processes such as RNA metabolism, cell cycle, apoptosis, and others. Unsurprisingly, given their relevance for cellular and organismal functions, alterations of epitranscriptomic marks have been observed in a broad range of human diseases, including cancer, neurological and metabolic disorders. Here, we will review the major types of mRNA modifications and editing processes in conjunction with the enzymes involved in their metabolism and describe their impact on human diseases. We present the current knowledge in an updated catalog. We will also discuss the emerging evidence on the crosstalk of epitranscriptomic marks and what this interplay could imply for the dynamics of mRNA modifications. Understanding how this complex regulatory layer can affect the course of human pathologies will ultimately lead to its exploitation toward novel epitranscriptomic therapeutic strategies.

The role of N6-methyladenosine modification in the life cycle and disease pathogenesis of hepatitis B and C viruses

N6-methyladenosine (m6A) is the most prevalent modification of mammalian cellular RNAs. m6A methylation is linked to epigenetic regulation of several aspects of gene expression, including RNA stability, splicing, nuclear export, RNA folding, and translational activity. m6A modification is reversibly catalyzed by methyltransferases (m6A writers) and demethylases (m6A erasers), and the dynamics of m6A-modified RNA are regulated by m6A-binding proteins (m6A readers). Recently, several studies have shown that m6A methylation sites have been identified in hepatitis B virus (HBV) transcripts and the hepatitis C virus (HCV) RNA genome. Here, we review the role of m6A modification in HBV/HCV replication and its contribution to liver disease pathogenesis. A better understanding of the functions of m6A methylation in the life cycles of HBV and HCV is required to establish the role of these modifications in liver diseases associated with these viral infections.

Biased Mutation and Selection in RNA Viruses

RNA viruses are responsible for some of the worst pandemics known to mankind, including outbreaks of Influenza, Ebola, and COVID-19. One major challenge in tackling RNA viruses is the fact they are extremely genetically diverse. Nevertheless, they share common features that include their dependence on host cells for replication, and high mutation rates. We set out to search for shared evolutionary characteristics that may aid in gaining a broader understanding of RNA virus evolution, and constructed a phylogeny-based data set spanning thousands of sequences from diverse single-stranded RNA viruses of animals. Strikingly, we found that the vast majority of these viruses have a skewed nucleotide composition, manifested as adenine rich (A-rich) coding sequences. In order to test whether A-richness is driven by selection or by biased mutation processes, we harnessed the effects of incomplete purifying selection at the tips of virus phylogenies. Our results revealed consistent mutational biases toward U rather than A in genomes of all viruses. In +ssRNA viruses, we found that this bias is compensated by selection against U and selection for A, which leads to A-rich genomes. In -ssRNA viruses, the genomic mutational bias toward U on the negative strand manifests as A-rich coding sequences, on the positive strand. We investigated possible reasons for the advantage of A-rich sequences including weakened RNA secondary structures, codon usage bias, and selection for a particular amino acid composition, and conclude that host immune pressures may have led to similar biases in coding sequence composition across very divergent RNA viruses.

Adenosine-to-inosine RNA editing in neurological development and disease

Adenosine-to-inosine (A-to-I) editing is one of the most prevalent post-transcriptional RNA modifications in metazoan. This reaction is catalysed by enzymes called adenosine deaminases acting on RNA (ADARs). RNA editing is involved in the regulation of protein function and gene expression. The numerous A-to-I editing sites have been identified in both coding and non-coding RNA transcripts. These editing sites are also found in various genes expressed in the central nervous system (CNS) and play an important role in neurological development and brain function. Aberrant regulation of RNA editing has been associated with the pathogenesis of neurological and psychiatric disorders, suggesting the physiological significance of RNA editing in the CNS. In this review, we discuss the current knowledge of editing on neurological disease and development.

MicroRNAs and SARS-CoV-2 life cycle, pathogenesis, and mutations: biomarkers or therapeutic agents?

To date, proposed therapies and antiviral drugs have been failed to cure coronavirus disease 2019 (COVID-19) patients. However, at least two drug companies have applied for emergency use authorization with the United States Food and Drug Administration for their coronavirus vaccine candidates and several other vaccines are in various stages of development to determine safety and efficacy. Recently, some studies have shown the role of different human and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) microRNAs (miRNAs) in the pathophysiology of COVID-19. miRNAs are non-coding single-stranded RNAs, which are involved in several physiological and pathological conditions, such as cell proliferation, differentiation, and metabolism. They act as negative regulators of protein synthesis through binding to the 3' untranslated region (3' UTR) of the complementary target mRNA, leading to mRNA degradation or inhibition. The databases of Google Scholar, Scopus, PubMed, and Web of Science were searched for literature regarding the importance of miRNAs in the SARS-CoV-2 life cycle, pathogenesis, and genomic mutations. Furthermore, promising miRNAs as a biomarker or antiviral agent in COVID-19 therapy are reviewed.

Drivers of within-host genetic diversity in acute infections of viruses

Genetic diversity is the fuel of evolution and facilitates adaptation to novel environments. However, our understanding of what drives differences in the genetic diversity during the early stages of viral infection is somewhat limited. Here, we use ultra-deep sequencing to interrogate 43 clinical samples taken from early infections of the human-infecting viruses HIV, RSV and CMV. Hundreds to thousands of virus templates were sequenced per sample, allowing us to reveal dramatic differences in within-host genetic diversity among virus populations. We found that increased diversity was mostly driven by presence of multiple divergent genotypes in HIV and CMV samples, which we suggest reflect multiple transmitted/founder viruses. Conversely, we detected an abundance of low frequency hyper-edited genomes in RSV samples, presumably reflecting defective virus genomes (DVGs). We suggest that RSV is characterized by higher levels of cellular co-infection, which allow for complementation and hence elevated levels of DVGs.

The few days or weeks following infection with a virus, termed acute infection, are critical for virus establishment. Here we sought to characterize what leads to differences in the genetic diversity of different viruses sampled during acute infection. We performed ultra-deep sequencing of hundreds to thousands viral genomes from forty-three samples spanning three pathogenic human viruses: HIV, RSV and CMV. We found major differences in the genetic diversity of these different viruses, and in different patients infected with the same virus. We investigated the factors responsible for these differences. We found that the DNA virus CMV was less diverse, most likely since it has a lower mutation rate than the RNA viruses HIV and RSV. We also found that the samples with the highest genetic diversity, which included one CMV sample and two HIV samples, bore evidence for multiple genotype infection. In other words, patients from whom these samples were taken were infected with two different “strains” of the virus. Finally, we also found evidence that viral genomes of HIV, and in particular RSV, are edited by the innate immune system of the host, leading to the presence of defective virus genomes.

Species-Specific Evolution of Ebola Virus during Replication in Human and Bat Cells

Ebola virus (EBOV) causes a severe, often fatal disease in humans and nonhuman primates. Within the past decade, EBOV has caused two large and difficult-to-control outbreaks, one of which recently ended in the Democratic Republic of the Congo. Bats are the likely reservoir of EBOV, but little is known of their relationship with the virus. We perform serial passages of EBOV in human and bat cells and use circular sequencing to compare the short-term evolution of the virus. Virus populations passaged in bat cells have sequence markers indicative of host RNA editing enzyme activity, including evidence for ADAR editing of the EBOV glycoprotein. Multiple regions in the EBOV genome appear to have undergone adaptive evolution when passaged in bat and human cells. Individual mutated viruses are rescued and characterized. Our results provide insight into the host species-specific evolution of EBOV and highlight the adaptive flexibility of the virus.

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

The COVID-19 outbreak has become a global health risk, and understanding the response of the host to the SARS-CoV-2 virus will help to combat the disease. RNA editing by host deaminases is an innate restriction process to counter virus infection, but it is not yet known whether this process operates against coronaviruses. Here, we analyze RNA sequences from bronchoalveolar lavage fluids obtained from coronavirus-infected patients. We identify nucleotide changes that may be signatures of RNA editing: adenosine-to-inosine changes from ADAR deaminases and cytosine-to-uracil changes from APOBEC deaminases. Mutational analysis of genomes from different strains of Coronaviridae from human hosts reveals mutational patterns consistent with those observed in the transcriptomic data. However, the reduced ADAR signature in these data raises the possibility that ADARs might be more effective than APOBECs in restricting viral propagation. Our results thus suggest that both APOBECs and ADARs are involved in coronavirus genome editing, a process that may shape the fate of both virus and patient.

Relative Contributions of the cGAS-STING and TLR3 Signaling Pathways to Attenuation of Herpes Simplex Virus 1 Replication

The innate immune response is crucial for defense against viral infections. Cells recognize virus infection through pattern recognition receptors and induce type I interferons as well as proinflammatory cytokines to orchestrate an innate immune response. Herpes simplex virus 1 (HSV-1) triggers both the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) and Toll-like receptor 3 (TLR3) pathways. It is well known that TLR3 uses the adaptor protein Toll/interleukin-1 receptor (IL-1R) domain-containing adaptor-inducing beta interferon (TRIF) for signaling, but we recently reported that STING signaling also requires TRIF. Because STING directly binds to TRIF, we identified the STING-interacting domain of TRIF and generated STING-noninteracting mutants of human and mouse TRIFs. The mutant TRIFs were unable to support STING signaling, although they were fully functional in the TLR3 pathway. These mutants were used to assess the relative contributions of the TLR3 and STING pathways to the attenuation of HSV-1 replication in mouse and human cell lines. For this purpose, the mouse L929 and NB41A3 cell lines and the human HT1080 and HeLa-M cell lines, in which both the TLR3 and the STING pathways are operational, were used. The TRIF gene was disrupted in these lines by CRISPR/Cas9, before reconstituting them with mutant and wild-type TRIF expression vectors. Infection of the reconstituted cells with HSV-1 revealed that both the cGAS-STING and the TLR3 signaling pathways are required for the attenuation of virus replication, but their relative contributions in attenuating HSV-1 replication were found to be different in mouse versus human cell lines. Thus, our study suggests that the relative contributions of the cGAS-STING and the TLR3 pathways in the attenuation of viral infection may be species specific.IMPORTANCE The magnitude of fatal infections caused by all different viruses in human and animal populations justifies a better understanding of the host innate immune response process that attenuates virus replication. In particular, the relative contributions of different signaling pathways which are responsible for the generation of the innate immune response are still largely unknown. In this study, we used STING-noninteracting TRIF mutants to decipher the relative contributions of the TLR3 and cGAS-STING signaling pathways to the attenuation of HSV-1 infection. We show that the relative contributions of the two pathways to the attenuation of viral infection are different in mouse versus human cell lines. Together, our results provide new insights into the relative contributions of two different signaling pathways in the attenuation of viral infection and may lead to the development of new antiviral strategies aimed at blocking viral infection at very early stages.

Therapeutic targeting of the PI4K2A/PKR lysosome network is critical for misfolded protein clearance and survival in cancer cells

The role of RNA-dependent protein kinase R (PKR) and its association with misfolded protein expression in cancer cells are unclear. Herein we report that PKR regulates misfolded protein clearance by preventing it release through exosomes and promoting lysosomal degradation of misfolded prion proteins in cancer cells. We demonstrated that PKR contributes to the lysosome function and regulates misfolded prion protein clearance. We hypothesized that PKR-associated lysosome function is critical for cancer but not normal cell survival, representing an effective approach for highly targeted cancer therapy. In screening a compound library, we identified two PKR-associated compounds 1 and 2 (Pac 1 and 2) did not affect normal cells but selectively induced cell death in cancer cells depending on their PKR expression status. Pac 1 significantly inhibited the growth of human lung and breast xenograft tumors in mice with no toxicity. Pac 1 binds to PI4K2A and disrupts the PKR/PI4K2A-associated lysosome complex, contributing to destabilization of cancer cell lysosomes and triggering cell death. We observed that PKR and PI4K2A play significant prognostic roles in breast cancer patients. These results demonstrate that targeting of a PI4K2A/PKR lysosome complex may be an effective approach for cancer therapy.

RNA editing in the forefront of epitranscriptomics and human health

Post-transcriptional modifications have been recently expanded with the addition of RNA editing, which is predominantly mediated by adenosine and cytidine deaminases acting on DNA and RNA. Here, we review the full spectrum of physiological processes in which these modifiers are implicated, among different organisms. Adenosine to inosine (A-to-I) editors, members of the ADAR and ADAT protein families are important regulators of alternative splicing and transcriptional control. On the other hand, cytidine to uridine (C-to-U) editors, members of the AID/APOBEC family, are heavily implicated in innate and adaptive immunity with important roles in antibody diversification and antiviral response. Physiologically, these enzymes are present in the nucleus and/or the cytoplasm, where they modify various RNA molecules, including miRNAs, tRNAs apart from mRNAs, whereas DNA editing is also possible by some of them. The expansion of next generation sequencing technologies provided a wealth of data regarding such modifications. RNA editing has been implicated in various disorders including cancer, and neurological diseases of the brain or the central nervous system. It is also related to cancer heterogeneity and the onset of carcinogenesis. Response to treatment can also be affected by the RNA editing status where drug efficacy is significantly compromised. Studying RNA editing events can pave the way to the identification of new disease biomarkers, and provide a more personalised therapy to various diseases.

Circular RNA biogenesis is decreased in postmortem cortical gray matter in schizophrenia and may alter the bioavailability of associated miRNA

Circular RNAs (circRNAs) are a covalently closed subclass of non-coding RNA molecules formed by back splicing of linear precursor RNA. These molecules are relatively stable and particularly abundant in the mammalian brain and therefore may participate in neural development and function. With the emergence of circRNAs activity in gene regulation, these molecules have been implicated in several biological processes, including synaptic plasticity, and we therefore suspect they may have a role in neurobehavioral disorders. Here, we profile cortical circRNAs expression in 35 postmortem cortical gray matter (BA46) schizophrenia and a non-psychiatric comparison group, using circRNA enrichment sequencing. While more than 90,000 circRNAs species were identified in the dorsolateral prefrontal cortex (DLPFC), we observed lower complexity and substantial depletion in subjects with the disorder. Although circRNAs expression was independent of their host gene transcription, alternative splicing rates were lower in samples from cases compared to controls. Gene set analysis of differentially expressed circRNAs host genes revealed significant enrichment of neural functions and neurological disorders. Many of these depleted circRNAs are also predicted to sequester miRNAs that were shown previously to be increased in the disorder, potentially exacerbating the functional impact of their dysregulation through posttranscriptional gene silencing. While this is the first reported exploration of circRNAs in schizophrenia, there is significant potential for dysregulation more broadly in other major mental illnesses and behavioral disorders. Given their capacity for modulating miRNA function, circRNA may play a significant role in the pathophysiology of disease and even be targeted for therapeutic manipulation.

Adenosine Deaminase Acting on RNA 1 Associates with Orf Virus OV20.0 and Enhances Viral Replication

Orf virus (ORFV) infects sheep and goats and is also an important zoonotic pathogen. The viral protein OV20.0 has been shown to suppress innate immunity by targeting the double-stranded RNA (dsRNA)-activated protein kinase (PKR) by multiple mechanisms. These mechanisms include a direct interaction with PKR and binding with two PKR activators, dsRNA and the cellular PKR activator (PACT), which ultimately leads to the inhibition of PKR activation. In the present study, we identified a novel association between OV20.0 and adenosine deaminase acting on RNA 1 (ADAR1). OV20.0 bound directly to the dsRNA binding domains (RBDs) of ADAR1 in the absence of dsRNA. Additionally, OV20.0 preferentially interacted with RBD1 of ADAR1, which was essential for its dsRNA binding ability and for the homodimerization that is critical for intact adenosine-to-inosine (A-to-I)-editing activity. Finally, the association with OV20.0 suppressed the A-to-I-editing ability of ADAR1, while ADAR1 played a proviral role during ORFV infection by inhibiting PKR phosphorylation. These observations revealed a new strategy used by OV20.0 to evade antiviral responses via PKR.IMPORTANCE Viruses evolve specific strategies to counteract host innate immunity. ORFV, an important zoonotic pathogen, encodes OV20.0 to suppress PKR activation via multiple mechanisms, including interactions with PKR and two PKR activators. In this study, we demonstrated that OV20.0 interacts with ADAR1, a cellular enzyme responsible for converting adenosine (A) to inosine (I) in RNA. The RNA binding domains, but not the catalytic domain, of ADAR1 are required for this interaction. The OV20.0-ADAR1 association affects the functions of both proteins; OV20.0 suppressed the A-to-I editing of ADAR1, while ADAR1 elevated OV20.0 expression. The proviral role of ADAR1 is likely due to the inhibition of PKR phosphorylation. As RNA editing by ADAR1 contributes to the stability of the genetic code and the structure of RNA, these observations suggest that in addition to serving as a PKR inhibitor, OV20.0 might modulate ADAR1-dependent gene expression to combat antiviral responses or achieve efficient viral infection.

Chikungunya Virus Fidelity Variants Exhibit Differential Attenuation and Population Diversity in Cell Culture and Adult Mice

Chikungunya virus (CHIKV) is a reemerging global health threat that produces debilitating arthritis in people. Like other RNA viruses with high mutation rates, CHIKV produces populations of genetically diverse genomes within a host. While several known CHIKV mutations influence disease severity in vertebrates and transmission by mosquitoes, the role of intrahost diversity in chikungunya arthritic disease has not been studied. In this study, high- and low-fidelity CHIKV variants, previously characterized by altered in vitro population mutation frequencies, were used to evaluate how intrahost diversity influences clinical disease, CHIKV replication, and antibody neutralization in immunocompetent adult mice inoculated in the rear footpads. Both high- and low-fidelity mutations were hypothesized to attenuate CHIKV arthritic disease, replication, and neutralizing antibody levels compared to wild-type (WT) CHIKV. Unexpectedly, high-fidelity mutants elicited more severe arthritic disease than the WT despite comparable CHIKV replication, whereas a low-fidelity mutant produced attenuated disease and replication. Serum antibody developed against both high- and low-fidelity CHIKV exhibited reduced neutralization of WT CHIKV. Using next-generation sequencing (NGS), the high-fidelity mutations were demonstrated to be genetically stable but produced more genetically diverse populations than WT CHIKV in mice. This enhanced diversification was subsequently reproduced after serial in vitro passage. The NGS results contrast with previously reported population diversities for fidelity variants, which focused mainly on part of the E1 gene, and highlight the need for direct measurements of mutation rates to clarify CHIKV fidelity phenotypes.IMPORTANCE CHIKV is a reemerging global health threat that elicits debilitating arthritis in humans. There are currently no commercially available CHIKV vaccines. Like other RNA viruses, CHIKV has a high mutation rate and is capable of rapid intrahost diversification during an infection. In other RNA viruses, virus population diversity associates with disease progression; however, potential impacts of intrahost viral diversity on CHIKV arthritic disease have not been studied. Using previously characterized CHIKV fidelity variants, we addressed whether CHIKV population diversity influences the severity of arthritis and host antibody response in an arthritic mouse model. Our findings show that CHIKV populations with greater genetic diversity can cause more severe disease and stimulate antibody responses with reduced neutralization of low-diversity virus populations in vitro The discordant high-fidelity phenotypes in this study highlight the complexity of inferring replication fidelity indirectly from population diversity.

Abortigenic but Not Neurotropic Equine Herpes Virus 1 Modulates the Interferon Antiviral Defense

Equine herpesvirus 1 (EHV1) is considered as a major pathogen of Equidae, causing symptoms from mild respiratory disease to late-term abortion and neurological disorders. Different EHV1 strains circulating in the field have been characterized to be of abortigenic or neurovirulent phenotype. Both variants replicate in a plaque-wise manner in the epithelium of the upper respiratory tract (URT), where the abortigenic strains induce more prominent viral plaques, compared to the neurovirulent strains. Considering the differences in replication at the URT, we hypothesized that abortigenic strains may show an increased ability to modulate the type I IFN secretion/signaling pathway, compared to strains that display the neurovirulent phenotype. Here, we analyze IFN levels induced by abortigenic and neurovirulent EHV1 using primary respiratory epithelial cells (EREC) and respiratory mucosa ex vivo explants. Similar levels of IFNα (~70 U/ml) were detected in explants inoculated with both types of EHV1 strains from 48 to 72 hpi. Second, EREC and mucosa explants were treated with recombinant equine IFNα (rEqIFNα) or Ruxolitinib (Rux), an IFN signaling inhibitor, prior to and during inoculation with abortigenic or neurovirulent EHV1. Replication of both EHV1 variants was suppressed by rEqIFNα. Further, addition of Rux increased replication in a concentration-dependent manner, indicating an IFN-susceptibility for both variants. However, in two out of three horses, at a physiological concentration of 100 U/ml of rEqIFNα, an increase in abortigenic EHV1 replication was observed compared to 10 U/ml of rEqIFNα, which was not observed for the neurovirulent strains. Moreover, in the presence of Rux, the plaque size of the abortigenic variants remained unaltered, whereas the typically smaller viral plaques induced by the neurovirulent variants became larger. Overall, our results demonstrate the importance of IFNα in the control of EHV1 replication in the URT for both abortigenic and neurovirulent variants. In addition, our findings support the speculation that abortigenic variants of EHV1 may have developed anti-IFN mechanisms that appear to be absent or less pronounced in neurovirulent EHV1 strains.

Metagenomic sequencing suggests a diversity of RNA interference-like responses to viruses across multicellular eukaryotes

RNA interference (RNAi)-related pathways target viruses and transposable element (TE) transcripts in plants, fungi, and ecdysozoans (nematodes and arthropods), giving protection against infection and transmission. In each case, this produces abundant TE and virus-derived 20-30nt small RNAs, which provide a characteristic signature of RNAi-mediated defence. The broad phylogenetic distribution of the Argonaute and Dicer-family genes that mediate these pathways suggests that defensive RNAi is ancient, and probably shared by most animal (metazoan) phyla. Indeed, while vertebrates had been thought an exception, it has recently been argued that mammals also possess an antiviral RNAi pathway, although its immunological relevance is currently uncertain and the viral small RNAs (viRNAs) are not easily detectable. Here we use a metagenomic approach to test for the presence of viRNAs in five species from divergent animal phyla (Porifera, Cnidaria, Echinodermata, Mollusca, and Annelida), and in a brown alga-which represents an independent origin of multicellularity from plants, fungi, and animals. We use metagenomic RNA sequencing to identify around 80 virus-like contigs in these lineages, and small RNA sequencing to identify viRNAs derived from those viruses. We identified 21U small RNAs derived from an RNA virus in the brown alga, reminiscent of plant and fungal viRNAs, despite the deep divergence between these lineages. However, contrary to our expectations, we were unable to identify canonical (i.e. Drosophila- or nematode-like) viRNAs in any of the animals, despite the widespread presence of abundant micro-RNAs, and somatic transposon-derived piwi-interacting RNAs. We did identify a distinctive group of small RNAs derived from RNA viruses in the mollusc. However, unlike ecdysozoan viRNAs, these had a piRNA-like length distribution but lacked key signatures of piRNA biogenesis. We also identified primary piRNAs derived from putatively endogenous copies of DNA viruses in the cnidarian and the echinoderm, and an endogenous RNA virus in the mollusc. The absence of canonical virus-derived small RNAs from our samples may suggest that the majority of animal phyla lack an antiviral RNAi response. Alternatively, these phyla could possess an antiviral RNAi response resembling that reported for vertebrates, with cryptic viRNAs not detectable through simple metagenomic sequencing of wild-type individuals. In either case, our findings show that the antiviral RNAi responses of arthropods and nematodes, which are highly divergent from each other and from that of plants and fungi, are also highly diverged from the most likely ancestral metazoan state.

Natural Variation in Resistance to Virus Infection in Dipteran Insects

The power and ease of Drosophila genetics and the medical relevance of mosquito-transmitted viruses have made dipterans important model organisms in antiviral immunology. Studies of virus-host interactions at the molecular and population levels have illuminated determinants of resistance to virus infection. Here, we review the sources and nature of variation in antiviral immunity and virus susceptibility in model dipteran insects, specifically the fruit fly Drosophila melanogaster and vector mosquitoes of the genera Aedes and Culex. We first discuss antiviral immune mechanisms and describe the virus-specificity of these responses. In the following sections, we review genetic and microbiota-dependent variation in antiviral immunity. In the final sections, we explore less well-studied sources of variation, including abiotic factors, sexual dimorphism, infection history, and endogenous viral elements. We borrow from work on other pathogen types and non-dipteran species when it parallels or complements studies in dipterans. Understanding natural variation in virus-host interactions may lead to the identification of novel restriction factors and immune mechanisms and shed light on the molecular determinants of vector competence.

Codon usage of HIV regulatory genes is not determined by nucleotide composition

Codon usage bias can be a result of either mutational bias or selection for translational efficiency and/or accuracy. Previous data has suggested that nucleotide composition constraint was the main determinant of HIV codon usage, and that nucleotide composition and codon usage were different between the regulatory genes, tat and rev, and other viral genes. It is not clear whether translational selection contributed to the codon usage difference and how nucleotide composition and translational selection interact to determine HIV codon usage. In this study, a model of codon bias due to GC composition with modification for the A-rich third codon position was used to calculate predicted HIV codon frequencies based on its nucleotide composition. The predicted codon usage of each gene was compared with the actual codon frequency. The predicted codon usage based on GC composition matched well with the actual codon frequencies for the structural genes (gag, pol and env). However, the codon usage of the regulatory genes (tat and rev) could not be predicted. Codon usage of the regulatory genes was also relatively unbiased showing the highest effective number of codons (ENC). Moreover, the codon adaptation index (CAI) of the regulatory genes showed better adaptation to human codons when compared to other HIV genes. Therefore, the early expressed genes responsible for regulation of the replication cycle, tat and rev, were more similar to humans in terms of codon usage and GC content than other HIV genes. This may help these genes to be expressed efficiently during the early stages of infection.

A-to-I RNA editing - thinking beyond the single nucleotide

Adenosine-to-inosine RNA editing is a conserved process, which is performed by ADAR enzymes. By changing nucleotides in coding regions of genes and altering codons, ADARs expand the cell's protein repertoire. This function of the ADAR enzymes is essential for human brain development. However, most of the known editing sites are in non-coding repetitive regions in the transcriptome and the purpose of editing in these regions is unclear. Recent studies, which have shown that editing levels of transcripts vary between tissues and developmental stages in many organisms, suggest that the targeted RNA and ADAR editing are both regulated. We discuss the implications of these findings, and the possible role of RNA editing in innate immunity.

Methods and Applications of CRISPR-Mediated Base Editing in Eukaryotic Genomes

The past several years have seen an explosion in development of applications for the CRISPR-Cas9 system, from efficient genome editing, to high-throughput screening, to recruitment of a range of DNA and chromatin-modifying enzymes. While homology-directed repair (HDR) coupled with Cas9 nuclease cleavage has been used with great success to repair and re-write genomes, recently developed base-editing systems present a useful orthogonal strategy to engineer nucleotide substitutions. Base editing relies on recruitment of cytidine deaminases to introduce changes (rather than double-stranded breaks and donor templates) and offers potential improvements in efficiency while limiting damage and simplifying the delivery of editing machinery. At the same time, these systems enable novel mutagenesis strategies to introduce sequence diversity for engineering and discovery. Here, we review the different base-editing platforms, including their deaminase recruitment strategies and editing outcomes, and compare them to other CRISPR genome-editing technologies. Additionally, we discuss how these systems have been applied in therapeutic, engineering, and research settings. Lastly, we explore future directions of this emerging technology.

Methods and Applications of CRISPR-Mediated Base Editing in Eukaryotic Genomes

The past several years have seen an explosion in development of applications for the CRISPR-Cas9 system, from efficient genome editing, to high-throughput screening, to recruitment of a range of DNA and chromatin-modifying enzymes. While homology-directed repair (HDR) coupled with Cas9 nuclease cleavage has been used with great success to repair and re-write genomes, recently developed base-editing systems present a useful orthogonal strategy to engineer nucleotide substitutions. Base editing relies on recruitment of cytidine deaminases to introduce changes (rather than double-stranded breaks and donor templates) and offers potential improvements in efficiency while limiting damage and simplifying the delivery of editing machinery. At the same time, these systems enable novel mutagenesis strategies to introduce sequence diversity for engineering and discovery. Here, we review the different base-editing platforms, including their deaminase recruitment strategies and editing outcomes, and compare them to other CRISPR genome-editing technologies. Additionally, we discuss how these systems have been applied in therapeutic, engineering, and research settings. Lastly, we explore future directions of this emerging technology.

Flavivirus and Filovirus EvoPrinters: New alignment tools for the comparative analysis of viral evolution

Background

Flavivirus and Filovirus infections are serious epidemic threats to human populations. Multi-genome comparative analysis of these evolving pathogens affords a view of their essential, conserved sequence elements as well as progressive evolutionary changes. While phylogenetic analysis has yielded important insights, the growing number of available genomic sequences makes comparisons between hundreds of viral strains challenging. We report here a new approach for the comparative analysis of these hemorrhagic fever viruses that can superimpose an unlimited number of one-on-one alignments to identify important features within genomes of interest.

Methodology/Principal finding

We have adapted EvoPrinter alignment algorithms for the rapid comparative analysis of Flavivirus or Filovirus sequences including Zika and Ebola strains. The user can input a full genome or partial viral sequence and then view either individual comparisons or generate color-coded readouts that superimpose hundreds of one-on-one alignments to identify unique or shared identity SNPs that reveal ancestral relationships between strains. The user can also opt to select a database genome in order to access a library of pre-aligned genomes of either 1,094 Flaviviruses or 460 Filoviruses for rapid comparative analysis with all database entries or a select subset. Using EvoPrinter search and alignment programs, we show the following: 1) superimposing alignment data from many related strains identifies lineage identity SNPs, which enable the assessment of sublineage complexity within viral outbreaks; 2) whole-genome SNP profile screens uncover novel Dengue2 and Zika recombinant strains and their parental lineages; 3) differential SNP profiling identifies host cell A-to-I hyper-editing within Ebola and Marburg viruses, and 4) hundreds of superimposed one-on-one Ebola genome alignments highlight ultra-conserved regulatory sequences, invariant amino acid codons and evolutionarily variable protein-encoding domains within a single genome.

Conclusions/Significance

EvoPrinter allows for the assessment of lineage complexity within Flavivirus or Filovirus outbreaks, identification of recombinant strains, highlights sequences that have undergone host cell A-to-I editing, and identifies unique input and database SNPs within highly conserved sequences. EvoPrinter’s ability to superimpose alignment data from hundreds of strains onto a single genome has allowed us to identify unique Zika virus sublineages that are currently spreading in South, Central and North America, the Caribbean, and in China. This new set of integrated alignment programs should serve as a useful addition to existing tools for the comparative analysis of these viruses.

Flaviviruses, including Zika and Dengue viruses, and Filoviruses, including Ebola and Marburg viruses, are significant global public health threats. Genetic surveillance of viral isolates provides important insights into the origin of outbreaks, reveals lineage heterogeneity and diversification, and facilitates identification of novel recombinant strains and host cell modified viral genomes. We report the development of EvoPrinter, a web-accessed alignment tool for the rapid comparative analysis of viral genomes. EvoPrinter superimposes alignment data from multiple pairwise comparisons onto a single reference sequence of interest, to reveal both similarities and differences detected in hundreds of selected viral isolates. Evoprinter databases provide easy access to hundreds of non-redundant Flavivirus and Filovirus genomes. allowing the user to distinguish between sublineage identity SNPs and unique strain-specific SNPs, thus facilitating analysis of the history of viral diversification during an epidemic. EvoPrinter also proves useful in identifying recombinant strains and their parental lineages and detecting host-cell genomic editing. EvoPrinter should serve as a useful addition to existing tools for the comparative analysis of these viruses.

Methods for the Detection of Adenosine-to-Inosine Editing Events in Cellular RNA

Modification of RNA is essential for properly expressing the repertoire of RNA transcripts necessary for both cell type and developmental specific functions. RNA modifications serve to dynamically re-wire and fine-tune the genetic information carried by an invariable genome. One important type of RNA modification is RNA editing and the most common and well-studied type of RNA editing is the hydrolytic deamination of adenosine to inosine. Inosine is a biological mimic of guanosine; therefore, when RNA is reverse transcribed, inosine is recognized as guanosine by the reverse transcriptase and a cytidine is incorporated into the complementary DNA (cDNA) strand. During PCR amplification, guanosines pair with the newly incorporated cytidines. As a result, the adenosine-to-inosine (A-to-I) editing events are recognized as adenosine to guanosine changes when comparing the sequences of the genomic DNA to the cDNA. This chapter describes the methods for extracting endogenous RNA for subsequent analyses of A-to-I RNA editing using reverse transcriptase-based approaches. We discuss techniques for the detection of A-to-I RNA editing events in messenger RNA (mRNA), including analyzing editing levels at specific adenosines within the total pool of mRNA versus analyzing editing patterns that occur in individual transcripts and a method for detecting editing events across the entire transcriptome. The detection of RNA editing events and editing levels can be used to better understand normal biological processes and disease states.

Genetic Instability of RNA Viruses

Despite having very limited coding capacity, RNA viruses are able to withstand challenge of antiviral drugs, cause epidemics in previously exposed human populations, and, in some cases, infect multiple host species. They are able to achieve this by virtue of their ability to multiply very rapidly, coupled with their extraordinary degree of genetic heterogeneity. RNA viruses exist not as single genotypes, but as a swarm of related variants, and this genomic diversity is an essential feature of their biology. RNA viruses have a variety of mechanisms that act in combination to determine their genetic heterogeneity. These include polymerase fidelity, error-mitigation mechanisms, genomic recombination, and different modes of genome replication. RNA viruses can vary in their ability to tolerate mutations, or “genetic robustness,” and several factors contribute to this. Finally, there is evidence that some RNA viruses exist close to a threshold where polymerase error rate has evolved to maximize the possible sequence space available, while avoiding the accumulation of a lethal load of deleterious mutations. We speculate that different viruses have evolved different error rates to complement the different “life-styles” they possess.

Diverse selective regimes shape genetic diversity at ADAR genes and at their coding targets

A-to-I RNA editing operated by ADAR enzymes is extremely common in mammals. Several editing events in coding regions have pivotal physiological roles and affect protein sequence (recoding events) or function. We analyzed the evolutionary history of the 3 ADAR family genes and of their coding targets. Evolutionary analysis indicated that ADAR evolved adaptively in primates, with the strongest selection in the unique N-terminal domain of the interferon-inducible isoform. Positively selected residues in the human lineage were also detected in the ADAR deaminase domain and in the RNA binding domains of ADARB1 and ADARB2. During the recent history of human populations distinct variants in the 3 genes increased in frequency as a result of local selective pressures. Most selected variants are located within regulatory regions and some are in linkage disequilibrium with eQTLs in monocytes. Finally, analysis of conservation scores of coding editing sites indicated that editing events are counter-selected within regions that are poorly tolerant to change. Nevertheless, a minority of recoding events occurs at highly conserved positions and possibly represents the functional fraction. These events are enriched in pathways related to HIV-1 infection and to epidermis/hair development. Thus, both ADAR genes and their targets evolved under variable selective regimes, including purifying and positive selection. Pressures related to immune response likely represented major drivers of evolution for ADAR genes. As for their coding targets, we suggest that most editing events are slightly deleterious, although a minority may be beneficial and contribute to antiviral response and skin homeostasis.

A-to-I editing of coding and non-coding RNAs by ADARs

Adenosine deaminases acting on RNA (ADARs) convert adenosine to inosine in double-stranded RNA. This A-to-I editing occurs not only in protein-coding regions of mRNAs, but also frequently in non-coding regions that contain inverted Alu repeats. Editing of coding sequences can result in the expression of functionally altered proteins that are not encoded in the genome, whereas the significance of Alu editing remains largely unknown. Certain microRNA (miRNA) precursors are also edited, leading to reduced expression or altered function of mature miRNAs. Conversely, recent studies indicate that ADAR1 forms a complex with Dicer to promote miRNA processing, revealing a new function of ADAR1 in the regulation of RNA interference.

Activation of Tomato Bushy Stunt Virus RNA-Dependent RNA Polymerase by Cellular Heat Shock Protein 70 Is Enhanced by Phospholipids In Vitro

Similar to other positive-strand RNA viruses, tombusviruses are replicated by the membrane-bound viral replicase complex (VRC). The VRC consists of the p92 virus-coded RNA-dependent RNA polymerase (RdRp), the viral p33 RNA chaperone, and several co-opted host proteins. In order to become a functional RdRp after its translation, the p92 replication protein should be incorporated into the VRC, followed by its activation. We have previously shown in a cell-free yeast extract-based assay that the activation of the Tomato bushy stunt virus (TBSV) RdRp requires a soluble host factor(s). In this article, we identify the cellular heat shock protein 70 (Hsp70) as the co-opted host factor required for the activation of an N-terminally truncated recombinant TBSV RdRp. In addition, small-molecule-based blocking of Hsp70 function inhibits RNA synthesis by the tombusvirus RdRp in vitro. Furthermore, we show that neutral phospholipids, namely, phosphatidylethanolamine (PE) and phosphatidylcholine (PC), enhance RdRp activation in vitro. In contrast, phosphatidylglycerol (PG) shows a strong and dominant inhibitory effect on in vitro RdRp activation. We also demonstrate that PE and PC stimulate RdRp-viral plus-strand RNA [(+)RNA] interaction, while PG inhibits the binding of the viral RNA to the RdRp. Based on the stimulatory versus inhibitory roles of various phospholipids in tombusvirus RdRp activation, we propose that the lipid composition of targeted subcellular membranes might be utilized by tombusviruses to regulate new VRC assembly during the course of infection.

IMPORTANCE

The virus-coded RNA-dependent RNA polymerase (RdRp), which is responsible for synthesizing the viral RNA progeny in infected cells of several positive-strand RNA viruses, is initially inactive. This strategy is likely to avoid viral RNA synthesis in the cytosol that would rapidly lead to induction of RNA-triggered cellular antiviral responses. During the assembly of the membrane-bound replicase complex, the viral RdRp becomes activated through an incompletely understood process that makes the RdRp capable of RNA synthesis. By using TBSV RdRp, we show that the co-opted cellular Hsp70 chaperone and neutral phospholipids facilitate RdRp activation in vitro. In contrast, phosphatidylglycerol (PG) has a dominant inhibitory effect on in vitro RdRp activation and RdRp-viral RNA interaction, suggesting that the membranous microdomain surrounding the RdRp greatly affects its ability for RNA synthesis. Thus, the activation of the viral RdRp likely depends on multiple host components in infected cells.

ADAR1 enhances HTLV-1 and HTLV-2 replication through inhibition of PKR activity

Background

The role of innate immunity in general and of type I interferon (IFN-I) in particular in HTLV-1 pathogenesis is still a matter of debate. ADAR1-p150 is an Interferon Stimulated Gene (ISG) induced by IFN-I that can edit viral RNAs. We therefore investigated whether it could play the role of an anti-HTLV factor.

Results

We demonstrate here that ADAR1 is also expressed in the absence of IFN stimulation in activated primary T-lymphocytes that are the natural target of this virus and in HTLV-1 or HTLV-2 chronically infected T-cells. ADAR1 expression is also increased in primary lymphocytes obtained from HTLV-1 infected individuals. We show that ADAR1 enhances HTLV-1 and HTLV-2 infection in T-lymphocytes and that this proviral effect is independent from its editing activity. ADAR1 expression suppresses IFN-α inhibitory effect on HTLV-1 and HTLV-2 and acts through the repression of PKR phosphorylation.

Discussion

This study demonstrates that two interferon stimulated genes, i.e. PKR and ADAR1 have opposite effects on HTLV replication in vivo. The balanced expression of those proteins could determine the fate of the viral cycle in the course of infection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-014-0093-9) contains supplementary material, which is available to authorized users.

Template role of double-stranded RNA in tombusvirus replication

Replication of plus-strand RNA [(+)RNA] viruses of plants is a relatively simple process that involves complementary minus-strand RNA [(-)RNA] synthesis and subsequent (+)RNA synthesis. However, the actual replicative form of the (-)RNA template in the case of plant (+)RNA viruses is not yet established unambiguously. In this paper, using a cell-free replication assay supporting a full cycle of viral replication, we show that replication of Tomato bushy stunt virus (TBSV) leads to the formation of double-stranded RNA (dsRNA). Using RNase digestion, DNAzyme, and RNA mobility shift assays, we demonstrate the absence of naked (-)RNA templates during replication. Time course experiments showed the rapid appearance of dsRNA earlier than the bulk production of new (+)RNAs, suggesting an active role for dsRNA in replication. Radioactive nucleotide chase experiments showed that the mechanism of TBSV replication involves the use of dsRNA templates in strand displacement reactions, where the newly synthesized plus strand replaces the original (+)RNA in the dsRNA. We propose that the use of dsRNA as a template for (+)RNA synthesis by the viral replicase is facilitated by recruited host DEAD box helicases and the viral p33 RNA chaperone protein. Altogether, this replication strategy allows TBSV to separate minus- and plus-strand syntheses in time and regulate asymmetrical RNA replication that leads to abundant (+)RNA progeny.

IMPORTANCE

Positive-stranded RNA viruses of plants use their RNAs as the templates for replication. First, the minus strand is synthesized by the viral replicase complex (VRC), which then serves as a template for new plus-strand synthesis. To characterize the nature of the (-)RNA in the membrane-bound viral replicase, we performed complete RNA replication of Tomato bushy stunt virus (TBSV) in yeast cell-free extracts and in plant extracts. The experiments demonstrated that the TBSV (-)RNA is present as a double-stranded RNA that serves as the template for TBSV replication. During the production of new plus strands, the viral replicase displaces the old plus strand in the dsRNA template, leading to asymmetrical RNA synthesis. The presented data are in agreement with the model that the dsRNA is present in nuclease-resistant membranous VRCs. This strategy likely allows TBSV to protect the replicating viral RNA from degradation as well as to evade the early detection of viral dsRNAs by the host surveillance system.

MicroRNA editing facilitates immune elimination of HCMV infected cells

The human cytomegalovirus (HCMV) is extremely prevalent in the human population. Infection by HCMV is life threatening in immune compromised individuals and in immune competent individuals it can cause severe birth defects, developmental retardation and is even associated with tumor development. While numerous mechanisms were developed by HCMV to interfere with immune cell activity, much less is known about cellular mechanisms that operate in response to HCMV infection. Here we demonstrate that in response to HCMV infection, the expression of the short form of the RNA editing enzyme ADAR1 (ADAR1-p110) is induced. We identified the specific promoter region responsible for this induction and we show that ADAR1-p110 can edit miR-376a. Accordingly, we demonstrate that the levels of the edited-miR-376a (miR-376a(e)) increase during HCMV infection. Importantly, we show that miR-376a(e) downregulates the immune modulating molecule HLA-E and that this consequently renders HCMV infected cells susceptible to elimination by NK cells.

KSHV 2.0: a comprehensive annotation of the Kaposi's sarcoma-associated herpesvirus genome using next-generation sequencing reveals novel genomic and functional features

Productive herpesvirus infection requires a profound, time-controlled remodeling of the viral transcriptome and proteome. To gain insights into the genomic architecture and gene expression control in Kaposi's sarcoma-associated herpesvirus (KSHV), we performed a systematic genome-wide survey of viral transcriptional and translational activity throughout the lytic cycle. Using mRNA-sequencing and ribosome profiling, we found that transcripts encoding lytic genes are promptly bound by ribosomes upon lytic reactivation, suggesting their regulation is mainly transcriptional. Our approach also uncovered new genomic features such as ribosome occupancy of viral non-coding RNAs, numerous upstream and small open reading frames (ORFs), and unusual strategies to expand the virus coding repertoire that include alternative splicing, dynamic viral mRNA editing, and the use of alternative translation initiation codons. Furthermore, we provide a refined and expanded annotation of transcription start sites, polyadenylation sites, splice junctions, and initiation/termination codons of known and new viral features in the KSHV genomic space which we have termed KSHV 2.0. Our results represent a comprehensive genome-scale image of gene regulation during lytic KSHV infection that substantially expands our understanding of the genomic architecture and coding capacity of the virus.

Kaposi's sarcoma-associated herpesvirus (KSHV) is a cancer-causing agent in immunocompromised patients that establishes long-lasting infections in its hosts. Initially described in 1994 and extensively studied ever since, KSHV molecular biology is understood in broad outline, but many detailed questions are still to be resolved. After almost two decades, specific aspects pertaining to the organization of the KSHV genome as well as the fate of the viral transcripts during the productive stages of infection remain unexplored. Here we use a systematic genome-wide approach to investigate changes in gene and protein expression during the productive stage of infection known as the lytic cycle. We found that the viral genome has a large coding capacity, capable of generating at least 45% more products than initially anticipated by bioinformatic analyses alone, and that it uses multiple strategies to expand its coding capacity well beyond what is determined solely by the DNA sequence of its genome. We also provide an expanded and highly detailed annotation of known and new genomic features in KSHV. We have termed this new architectural and functional annotation KSHV 2.0. Our results indicate that viral genomes are more complex than anticipated, and that they are subject to tight mechanisms of regulation to ensure correct gene expression.

Adenosine-to-inosine RNA editing and human disease

A-to-I RNA editing is a post-transcriptional modification that converts adenosines to inosines in both coding and noncoding RNA transcripts. It is catalyzed by ADAR (adenosine deaminase acting on RNA) enzymes, which exist throughout the body but are most prevalent in the central nervous system. Inosines exhibit properties that are most similar to those of guanosines. As a result, ADAR-mediated editing can post-transcriptionally alter codons, introduce or remove splice sites, or affect the base pairing of the RNA molecule with itself or with other RNAs. A-to-I editing is a mechanism that regulates and diversifies the transcriptome, but the full biological significance of ADARs is not understood. ADARs are highly conserved across vertebrates and are essential for normal development in mammals. Aberrant ADAR activity has been associated with a wide range of human diseases, including cancer, neurological disorders, metabolic diseases, viral infections and autoimmune disorders. ADARs have been shown to contribute to disease pathologies by editing of glutamate receptors, editing of serotonin receptors, mutations in ADAR genes, and by other mechanisms, including recently identified regulatory roles in microRNA processing. Advances in research into many of these diseases may depend on an improved understanding of the biological functions of ADARs. Here, we review recent studies investigating connections between ADAR-mediated RNA editing and human diseases.