Mutational Pressure in Zika Virus: Local ADAR-Editing Areas Associated with Pauses in Translation and Replication

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.

The hidden RNA code: implications of the RNA epitranscriptome in the context of viral infections

Emerging evidence highlights the multifaceted roles of the RNA epitranscriptome during viral infections. By modulating the modification landscape of viral and host RNAs, viruses enhance their propagation and elude host surveillance mechanisms. Here, we discuss how specific RNA modifications, in either host or viral RNA molecules, impact the virus-life cycle and host antiviral responses, highlighting the potential of targeting the RNA epitranscriptome for novel antiviral therapies.

Host-mediated RNA editing in viruses

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

Spatio-temporal dynamics of intra-host variability in SARS-CoV-2 genomes

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying variants of concern (VOC). Viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host single nucleotide variations (iSNVs). Analysing 1347 samples collected till June 2020, we recorded 16 410 iSNV sites throughout the SARS-CoV-2 genome. We found ∼42% of the iSNV sites to be reported as SNVs by 30 September 2020 in consensus sequences submitted to GISAID, which increased to ∼80% by 30th June 2021. Following this, analysis of another set of 1774 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) lineage-defining variations appeared as iSNVs before getting fixed in the population. Besides, mutations in RdRp as well as RNA-editing by APOBEC and ADAR deaminases seem to contribute to the differential prevalence of iSNVs in hosts. We also observe hyper-variability at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions.

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.

Mutational spectrum of SARS-CoV-2 during the global pandemic

Viruses accumulate mutations under the influence of natural selection and host-virus interactions. Through a systematic comparison of 351,525 full viral genome sequences collected during the recent COVID-19 pandemic, we reveal the spectrum of SARS-CoV-2 mutations. Unlike those of other viruses, the mutational spectrum of SARS-CoV-2 exhibits extreme asymmetry, with a much higher rate of C>U than U>C substitutions, as well as a higher rate of G>U than U>G substitutions. This suggests directional genome sequence evolution during transmission. The substantial asymmetry and directionality of the mutational spectrum enable pseudotemporal tracing of SARS-CoV-2 without prior information about the root sequence, collection time, and sampling region. This shows that the viral genome sequences collected in Asia are similar to the original genome sequence. Adjusted estimation of the dN/dS ratio accounting for the asymmetrical mutational spectrum also shows evidence of negative selection on viral genes, consistent with previous reports. Our findings provide deep insights into the mutational processes in SARS-CoV-2 viral infection and advance the understanding of the history and future evolution of the virus.

From RNA World to SARS-CoV-2: The Edited Story of RNA Viral Evolution

The current SARS-CoV-2 pandemic underscores the importance of understanding the evolution of RNA genomes. While RNA is subject to the formation of similar lesions as DNA, the evolutionary and physiological impacts RNA lesions have on viral genomes are yet to be characterized. Lesions that may drive the evolution of RNA genomes can induce breaks that are repaired by recombination or can cause base substitution mutagenesis, also known as base editing. Over the past decade or so, base editing mutagenesis of DNA genomes has been subject to many studies, revealing that exposure of ssDNA is subject to hypermutation that is involved in the etiology of cancer. However, base editing of RNA genomes has not been studied to the same extent. Recently hypermutation of single-stranded RNA viral genomes have also been documented though its role in evolution and population dynamics. Here, we will summarize the current knowledge of key mechanisms and causes of RNA genome instability covering areas from the RNA world theory to the SARS-CoV-2 pandemic of today. We will also highlight the key questions that remain as it pertains to RNA genome instability, mutations accumulation, and experimental strategies for addressing these questions.

Adenosine Deaminases Acting on RNA (ADARs) and Viral Infections

C6 deamination of adenosine (A) to inosine (I) in double-stranded RNA (dsRNA) is catalyzed by a family of enzymes known as ADARs (adenosine deaminases acting on RNA) encoded by three genes in mammals. Alternative promoters and splicing produce two ADAR1 proteins, an interferon-inducible cytoplasmic p150 and a constitutively expressed p110 that like ADAR2 is a nuclear enzyme. ADAR3 lacks deaminase activity. A-to-I editing occurs with both viral and cellular RNAs. Deamination activity is dependent on dsRNA substrate structure and regulatory RNA-binding proteins and ranges from highly site selective with hepatitis D RNA and glutamate receptor precursor messenger RNA (pre-mRNA) to hyperediting of measles virus and polyomavirus transcripts and cellular inverted Alu elements. Because I base-pairs as guanosine instead of A, editing can alter mRNA decoding, pre-mRNA splicing, and microRNA silencing. Editing also alters dsRNA structure, thereby suppressing innate immune responses including interferon production and action. Expected final online publication date for the Annual Review of Virology, Volume 8 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Translation-Associated Mutational U-Pressure in the First ORF of SARS-CoV-2 and Other Coronaviruses

Within 4 months of the ongoing COVID-19 pandemic caused by SARS-CoV-2, more than 250 nucleotide mutations have been detected in ORF1ab of the virus isolated from infected persons from different parts of the globe. These observations open up an obvious question about the rate and direction of mutational pressure for further vaccine and therapeutics designing. In this study, we did a comparative analysis of ORF1a and ORF1b by using the first isolate (Wuhan strain) as the parent sequence. We observed that most of the nucleotide mutations are C to U transitions. The rate of synonymous C to U transitions is significantly higher than the rate of non-synonymous ones, indicating negative selection on amino acid substitutions. Further, trends in nucleotide usage bias have been investigated in 49 coronaviruses species. A strong bias in nucleotide usage in fourfold degenerate sites toward uracil residues is seen in ORF1ab of all the studied coronaviruses: both in the ORF1a and in the ORF1b translated thanks to the programmed ribosomal frameshifting that has an efficiency of 14 – 45% in different species. A more substantial mutational U-pressure is observed in ORF1a than in ORF1b perhaps because ORF1a is translated more frequently than ORF1b. Mutational U-pressure is there even in ORFs that are not translated from genomic RNA plus strands, but the bias is weaker than in ORF1ab. Unlike other nucleotide mutations, mutational U-pressure caused by cytosine deamination, mostly occurring during the RNA plus strand replication and also translation, cannot be corrected by the proof-reading machinery of coronaviruses. The knowledge generated on the mutational U-pressure that becomes stronger during translation of viral RNA plus strands has implications for vaccine and nucleoside analog development for treating COVID-19 and other coronavirus infections.

It's the Little Things (in Viral RNA)

Chemical modifications of viral RNA are an integral part of the viral life cycle and are present in most classes of viruses. To date, more than 170 RNA modifications have been discovered in all types of cellular RNA. Only a few, however, have been found in viral RNA, and the function of most of these has yet to be elucidated. Those few we have discovered and whose functions we understand have a varied effect on each virus. They facilitate RNA export from the nucleus, aid in viral protein synthesis, recruit host enzymes, and even interact with the host immune machinery. The most common methods for their study are mass spectrometry and antibody assays linked to next-generation sequencing. However, given that the actual amount of modified RNA can be very small, it is important to pair meticulous scientific methodology with the appropriate detection methods and to interpret the results with a grain of salt. Once discovered, RNA modifications enhance our understanding of viruses and present a potential target in combating them. This review provides a summary of the currently known chemical modifications of viral RNA, the effects they have on viral machinery, and the methods used to detect them.

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.

Codon usage similarity between viral and some host genes suggests a codon-specific translational regulation

The codon usage pattern is a specific characteristic of each species; however, the codon usage of all of the genes in a genome is not uniform. Intriguingly, most viruses have codon usage patterns that are vastly different from the optimal codon usage of their hosts. How viral genes with different codon usage patterns are efficiently expressed during a viral infection is unclear. An analysis of the similarity between viral codon usage and the codon usage of the individual genes of a host genome has never been performed. In this study, we demonstrated that the codon usage of human RNA viruses is similar to that of some human genes, especially those involved in the cell cycle. This finding was substantiated by its concordance with previous reports of an upregulation at the protein level of some of these biological processes. It therefore suggests that some suboptimal viral codon usage patterns may actually be compatible with cellular translational machineries in infected conditions.

Epitranscriptomic marks: Emerging modulators of RNA virus gene expression

Epitranscriptomics, the study of posttranscriptional chemical moieties placed on RNA, has blossomed in recent years. This is due in part to the emergence of high‐throughput detection methods as well as the burst of discoveries showing biological function of select chemical marks. RNA modifications have been shown to affect RNA structure, localization, and functions such as alternative splicing, stabilizing transcripts, nuclear export, cap‐dependent and cap‐independent translation, microRNA biogenesis and binding, RNA degradation, and immune regulation. As such, the deposition of chemical marks on RNA has the unique capability to spatially and temporally regulate gene expression. The goal of this article is to present the exciting convergence of the epitranscriptomic and virology fields, specifically the deposition and biological impact of N7‐methylguanosine, ribose 2′‐O‐methylation, pseudouridine, inosine, N6‐methyladenosine, and 5‐methylcytosine epitranscriptomic marks on gene expression of RNA viruses.

This article is categorized under:

RNA in Disease and Development > RNA in Disease

RNA Interactions with Proteins and Other Molecules > Protein–RNA Interactions: Functional Implications

The history of mutational pressure changes during the evolution of adeno-associated viruses: A message to gene therapy and DNA-vaccine vectors designers

The use of virus-associated vectors for gene therapy and vaccination have emerged as safe and effective delivery system. Like all other genetic materials, these vehicles are also prone to spontaneous mutations. To understand what types of nucleotide mutations are expected in the vector, one needs to know distinct characteristics of mutational process in the corresponding virus. In this study we analyzed mutational pressure directions along the length of the genomes of all types of primate adeno-associated viruses (AAV) that are frequently used in gene therapy or DNA-vaccines. We observed clear evidences of transcription-associated mutational pressure in AAV: nucleotide usage biases are changing drastically after each of the three promoters: the higher the rate of transcription, the stronger the bias towards GC to AT mutations. Moreover, the usage of G decreased at the lower transcription rate (after P19 promoter) than the usage of C (after P40 promoter). Since nucleotide usage biases are retrospective indices, we created a scenario of changes in transcriptional map during the AAV evolution. Current mutational pressure directions are different for AAV types, while all of them demonstrate high rates of T to C transitions in the second long ORF. Since transcription rate and cell tropism are the main factors determining the preferable direction of nucleotide mutations in AAV, mutational pressure should be checked experimentally in DNA vectors before their final design with the aim to make the transferred gene more stable against those mutations.

Double-stranded RNA deaminase ADAR1 promotes the Zika virus replication by inhibiting the activation of protein kinase PKR

Zika virus (ZIKV) is a mosquito-borne flavivirus that has emerged as a threat to global health. The family of adenosine deaminases acting on dsRNA (ADARs) are human host factors important for the genetic diversity and evolution of ZIKV. Here, we further investigated the role of ADAR1 in ZIKV replication by utilizing CRISPR/Cas9-based gene editing and RNAi-based gene knockdown techniques. Both ADAR1 knockout and knockdown significantly reduced ZIKV RNA synthesis, protein levels, and viral titers in several human cell lines. Trans-complementation with the full-length ADAR1 form p150 or the shorter form p110 lacking the Zα domain restored viral replication levels suppressed by the ADAR1 knockout. Moreover, we observed that the nuclear p110 form was redistributed to the cytoplasm in response to ZIKV infection. ADAR1 was not involved in viral entry but promoted viral protein translation by impairing ZIKV-induced activation of protein kinase regulated by dsRNA (PKR). Of note, the RNA-editing activity of ADAR1 was not required to promote ZIKV replication. We also found that the proviral role of ADAR1 was partially mediated through its ability to suppress IFN production and PKR activation. Our work identifies ADAR1 as a proviral factor involved in ZIKV replication, suggesting that ADAR1 could be a potential antiviral target.

Explaining Pathogenicity of Congenital Zika and Guillain-Barré Syndromes: Does Dysregulation of RNA Editing Play a Role?

Previous studies of Zika virus (ZIKV) pathogenesis have focused primarily on virus-driven pathology and neurotoxicity, as well as host-related changes in cell proliferation, autophagy, immunity, and uterine function. It is now hypothesized that ZIKV pathogenesis arises instead as an (unintended) consequence of host innate immunity, specifically, as the side effect of an otherwise well-functioning machine. The hypothesis presented here suggests a new way of thinking about the role of host immune mechanisms in disease pathogenesis, focusing on dysregulation of post-transcriptional RNA editing as a candidate driver of a broad range of observed neurodevelopmental defects and neurodegenerative clinical symptoms in both infants and adults linked with ZIKV infections. The authors collect and synthesize existing evidence of ZIKV-mediated changes in the expression of adenosine deaminases acting on RNA (ADARs), known links between abnormal RNA editing and pathogenesis, as well as ideas for future research directions, including potential treatment strategies.

RNA regulatory processes in RNA virus biology

Numerous post‐transcriptional RNA processes play a major role in regulating the quantity, quality and diversity of gene expression in the cell. These include RNA processing events such as capping, splicing, polyadenylation and modification, but also aspects such as RNA localization, decay, translation, and non‐coding RNA‐associated regulation. The interface between the transcripts of RNA viruses and the various RNA regulatory processes in the cell, therefore, has high potential to significantly impact virus gene expression, regulation, cytopathology and pathogenesis. Furthermore, understanding RNA biology from the perspective of an RNA virus can shed considerable light on the broad impact of these post‐transcriptional processes in cell biology. Thus the goal of this article is to provide an overview of the richness of cellular RNA biology and how RNA viruses use, usurp and/or avoid the associated machinery to impact the outcome of infection.

This article is categorized under:

RNA in Disease and Development > RNA in Disease

ScanFold: an approach for genome-wide discovery of local RNA structural elements-applications to Zika virus and HIV

In addition to encoding RNA primary structures, genomes also encode RNA secondary and tertiary structures that play roles in gene regulation and, in the case of RNA viruses, genome replication. Methods for the identification of functional RNA structures in genomes typically rely on scanning analysis windows, where multiple partially-overlapping windows are used to predict RNA structures and folding metrics to deduce regions likely to form functional structure. Separate structural models are produced for each window, where the step size can greatly affect the returned model. This makes deducing unique local structures challenging, as the same nucleotides in each window can be alternatively base paired. We are presenting here a new approach where all base pairs from analysis windows are considered and weighted by favorable folding. This results in unique base pairing throughout the genome and the generation of local regions/structures that can be ranked by their propensity to form unusually thermodynamically stable folds. We applied this approach to the Zika virus (ZIKV) and HIV-1 genomes. ZIKV is linked to a variety of neurological ailments including microcephaly and Guillain-Barré syndrome and its (+)-sense RNA genome encodes two, previously described, functionally essential structured RNA regions. HIV, the cause of AIDS, contains multiple functional RNA motifs in its genome, which have been extensively studied. Our approach is able to successfully identify and model the structures of known functional motifs in both viruses, while also finding additional regions likely to form functional structures. All data have been archived at the RNAStructuromeDB (www.structurome.bb.iastate.edu), a repository of RNA folding data for humans and their pathogens.

An Evolutionary Insight into Zika Virus Strains Isolated in the Latin American Region

Zika virus (ZIKV) is an emerging pathogen member of the Flaviviridae family. ZIKV has spread rapidly in the Latin American region, causing hundreds of thousands of cases of ZIKV disease, as well as microcephaly in congenital infections. Detailed studies on the pattern of evolution of ZIKV strains have been extremely important to our understanding of viral survival, fitness, and evasion of the host’s immune system. For these reasons, we performed a comprehensive phylogenetic analysis of ZIKV strains recently isolated in the Americas. The results of these studies revealed evidence of diversification of ZIKV strains circulating in the Latin American region into at least five different genetic clusters. This diversification was also reflected in the different trends in dinucleotide bias and codon usage variation. Amino acid substitutions were found in E and prM proteins of the ZIKV strains isolated in this region, revealing the presence of novel genetic variants circulating in Latin America.

Mechanisms and Concepts in RNA Virus Population Dynamics and Evolution

RNA viruses are unique in their evolutionary capacity, exhibiting high mutation rates and frequent recombination. They rapidly adapt to environmental changes, such as shifts in immune pressure or pharmacological challenge. The evolution of RNA viruses has been brought into new focus with the recent developments of genetic and experimental tools to explore and manipulate the evolutionary dynamics of viral populations. These studies have uncovered new mechanisms that enable viruses to overcome evolutionary challenges in the environment and have emphasized the intimate relationship of viral populations with evolution. Here, we review some of the emerging viral and host mechanisms that underlie the evolution of RNA viruses. We also discuss new studies that demonstrate that the relationship between evolutionary dynamics and virus biology spans many spatial and temporal scales, affecting transmission dynamics within and between hosts as well as pathogenesis.

Amino acid compositions contribute to the proteins' evolution under the influence of their abundances and genomic GC content

Inconsistent results on the association between evolutionary rates and amino acid composition of proteins have been reported in eukaryotes. However, there are few studies of how amino acid composition can influence evolutionary rates in bacteria. Thus, we constructed linear regression models between composition frequencies of amino acids and evolutionary rates for bacteria. Compositions of all amino acids can on average explain 21.5% of the variation in evolutionary rates among 273 investigated bacterial organisms. In five model organisms, amino acid composition contributes more to variation in evolutionary rates than protein abundance, and frequency of optimal codons. The contribution of individual amino acid composition to evolutionary rate varies among organisms. The closer the GC-content of genome to its maximum or minimum, the better the correlation between the amino acid content and the evolutionary rate of proteins would appear in that genome. The types of amino acids that significantly contribute to evolutionary rates can be grouped into GC-rich and AT-rich amino acids. Besides, the amino acid with high composition also contributes more to evolutionary rates than amino acid with low composition in proteome. In summary, amino acid composition significantly contributes to the rate of evolution in bacterial organisms and this in turn is impacted by GC-content.

Deciphering the dark proteome of Chikungunya virus

Chikungunya virus (CHIKV) is a mosquito-borne alphavirus. The outbreak of CHIKV infection has been seen in many tropical and subtropical regions of the biosphere. Current reports evidenced that after outbreaks in 2005-06, the fitness of this virus propagating in Aedes albopictus enhanced due to the epistatic mutational changes in its envelope protein. In our study, we evaluated the prevalence of intrinsically disordered proteins (IDPs) and IDP regions (IDPRs) in CHIKV proteome. IDPs/IDPRs are known as members of a 'Dark Proteome' that defined as a set of polypeptide segments or whole protein without unique three-dimensional structure within the cellular milieu but with significant biological functions, such as cell cycle regulation, control of signaling pathways, and maintenance of viral proteomes. However, the intrinsically disordered aspects of CHIKV proteome and roles of IDPs/IDPRs in the pathogenic mechanism of this important virus have not been evaluated as of yet. There are no existing reports on the analysis of intrinsic disorder status of CHIKV. To fulfil this goal, we have analyzed the abundance and functionality of IDPs/IDPRs in CHIKV proteins, involved in the replication and maturation. It is likely that these IDPs/IDPRs can serve as novel targets for disorder based drug design.

Pregnancy and infection: using disease pathogenesis to inform vaccine strategy

Vaccination is the mainstay of preventative medicine for many infectious diseases. Pregnant women, unborn fetuses, and neonates represent three at-risk populations that can be simultaneously protected by strategic vaccination protocols. Because the pathogenesis of different infectious microbes varies based on tissue tropism, timing of infection, and host susceptibility, the goals of immunization are not uniform across all vaccines. Mechanistic understanding of infectious disease pathogenesis and immune responses is therefore essential to inform vaccine design and the implementation of appropriate immunization protocols that optimize protection of pregnant women, fetuses, and neonates.

A comprehensive study on cellular RNA editing activity in response to infections with different subtypes of influenza a viruses

Background

RNA editing is an important mechanism that expands the diversity and complexity of genetic codes. The conversions of adenosine (A) to inosine (I) and cytosine (C) to uridine (U) are two prominent types of RNA editing in animals. The roles of RNA editing events have been implicated in important biological pathways. Cellular RNA editing activity in response to influenza A virus infection has not been fully characterized in human and avian hosts. This study was designed as a big data analysis to investigate the role and response of RNA editing in epithelial cells during the course of infection with various subtypes of influenza A viruses.

Results

Using a bioinformatics pipeline modified from our previous study, we characterized the profiles of A-to-I and C-to-U RNA editing events in human epithelial cells during the course of influenza A virus infection. Our results revealed a striking diversity of A-to-I RNA editing activities in human epithelial cells in responses to different subtypes of influenza A viruses. The infection of H1N1 and H3N2 significantly up-regulated normalized A-to-I RNA editing levels in human epithelial cells, whereas that of H5N1 did not change it and H7N9 infection significantly down-regulated normalized A-to-I editing level in A549 cells. Next, the expression levels of ADAR and APOBEC enzymes responsible for A-to-I and C-to-U RNA editing during the course of virus infection were examined. The increase of A-to-I RNA editing activities in infections with some influenza A viruses (H1N1 and H3N2) is linked to the up-regulation of ADAR1 but not ADAR2. Further, the pattern recognition receptors of human epithelial cells infected with H1N1, H3N2, H5N1 and H7N9 were examined. Variable responsive changes in gene expression were observed with RIG-I like receptors and Toll like receptors. Finally, the effect of influenza A virus infection on cellular RNA editing activity was also analyzed in avian hosts.

Conclusion

This work represents the first comprehensive study of cellular RNA editing activity in response to different influenza A virus infections in human and avian hosts, highlighting the critical role of RNA editing in innate immune response and the pathogenicity of different subtypes of influenza A viruses.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4330-1) contains supplementary material, which is available to authorized users.

Transcription-associated mutational pressure in the Parvovirus B19 genome: Reactivated genomes contribute to the variability of viral populations

In this study we used non-overlapping parts of the two long open reading frames coding for nonstructural (NS) and capsid (VP) proteins of all available sequences of the Parvovirus B19 subgenotype 1a genome and found out that the rates of A to G, C to T and A to T mutations are higher in the first long reading frame (NS) of the virus than in the second one (VP). This difference in mutational pressure directions for two parts of the same viral genome can be explained by the fact of transcription of just the first long reading frame during the lifelong latency in nonerythroid cells. Adenine deamination (producing A to G and A to T mutations) and cytosine deamination (producing C to T mutations) occur more frequently in transcriptional bubbles formed by DNA "plus" strand of the first open reading frame. These mutations can be inherited only in case of reactivation of the infectious virus due to the help of Adenovirus that allows latent Parvovirus B19 to start transcription of the second reading frame and then to replicate its genome by the rolling circle mechanism using the specific origin. Results of this study provide evidence that the genomes reactivated from latency make significant contributions to the variability of Parvovirus B19.

Functional RNA during Zika virus infection

Zika virus (ZIKV; family Flaviviridae; genus Flavivirus) is a pathogenic mosquito-borne RNA virus that currently threatens human health in the Americas, large parts of Asia and occasionally elsewhere in the world. ZIKV infection is often asymptomatic but can cause severe symptoms including congenital microcephaly and Guillain-Barré syndrome. The positive single-stranded RNA genome of the mosquito-borne ZIKV requires effective replication in two evolutionary distinct hosts - mosquitoes and primates. In addition to some of the viral proteins, the ZIKV genomic RNA and functional RNAs produced thereof aid in the establishment of productive infection and the evasion of host cell antiviral responses. ZIKV has evolved to contain a nucleotide composition and RNA modifications, such as methylation and the formation of G-quadruplexes that allow effective replication in both hosts. Furthermore, a number of host factors interact with the viral genome to modulate RNA replication. Importantly, the ZIKV genome produces non-coding subgenomic flavivirus RNA (sfRNA) due to stalling of host 5'- 3' ribonucleases on viral RNA structures in the 3' untranslated region (UTR). This sfRNA (sfRNA) exerts important proviral functions such as antagonizing the innate interferon response and RNA interference. Here, we discuss the ZIKV genomic RNA and functional RNAs thereof to assess their significance during ZIKV infection. Understanding the details of the ZIKV infection cycle will aid in the development of effective antiviral strategies and safe vaccines.

Epigallocatechin gallate, an active green tea compound inhibits the Zika virus entry into host cells via binding the envelope protein

Emerging infections of Zika virus (ZIKV) are associated with serious consequences like microcephaly and Guillain-Barré syndrome. It leads to a situation of global health emergency and demand an intensive research investigation to develop safe and effective therapeutics. Various efforts have been made to reduce the pathological pressure of ZIKV, but no effective drug has been introduced against ZIKV infections. A recent study has reported the inhibition of ZIKV entry into the host cells by an active green tea ingredient, Epigallocatechin Gallate (EGCG) in Vero E6cells. The effect of EGCG seems remarkable but lacking the information of the mechanism of action. In this study, we have investigated the binding site (Site1) of EGCG on envelope protein and provided the insights into various interactions of molecule with the binding site using molecular docking studies. Further, using molecular dynamics approaches we proposed the possible associated mechanism of inhibition of ZIKV entry by EGCG molecule. EGCG has found to interact with several residues and providing stability to the protein conformations up to 50ns simulations.

RNA editing by the host ADAR system affects the molecular evolution of the Zika virus

Zika virus (ZIKV) is a mosquito‐transmitted flavivirus, linked to microcephaly and fetal death in humans. Here, we investigate whether host‐mediated RNA editing of adenosines (ADAR) plays a role in the molecular evolution of ZIKV. Using complete coding sequences for the ZIKV polyprotein, we show that potential ADAR substitutions are underrepresented at the ADAR‐resistant GA dinucleotides of both the positive and negative strands, that these changes are spatially and temporally clustered (as expected of ADAR editing) for three branches of the viral phylogeny, and that ADAR mutagenesis can be linked to its codon usage. Furthermore, resistant GA dinucleotides are enriched on the positive (but not negative) strand, indicating that the former is under stronger purifying selection than the latter. ADAR editing also affects the evolution of the rhabdovirus sigma. Our study now documents that host ADAR editing is a mutation and evolutionary force of positive‐ as well as negative‐strand RNA viruses.