Bioinformatics for Inosine: Tools and Approaches to Trace This Elusive RNA Modification.

Inosine is a nucleotide resulting from the deamination of adenosine in RNA. This chemical modification process, known as RNA editing, is typically mediated by a family of double-stranded RNA binding proteins named Adenosine Deaminase Acting on dsRNA (ADAR). While the presence of ADAR orthologs has been traced throughout the evolution of metazoans, the existence and extension of RNA editing have been characterized in a more limited number of animals so far. Undoubtedly, ADAR-mediated RNA editing plays a vital role in physiology, organismal development and disease, making the understanding of the evolutionary conservation of this phenomenon pivotal to a deep characterization of relevant biological processes. However, the lack of direct high-throughput methods to reveal RNA modifications at single nucleotide resolution limited an extended investigation of RNA editing. Nowadays, these methods have been developed, and appropriate bioinformatic pipelines are required to fully exploit this data, which can complement existing approaches to detect ADAR editing. Here, we review the current literature on the "bioinformatics for inosine" subject and we discuss future research avenues in the field.

Insights into ADAR gene complement, expression patterns, and RNA editing landscape in Chlamys farreri.

Adenosine Deaminases Acting on RNA (ADARs) are evolutionarily conserved enzymes known to convert adenosine to inosine in double-stranded RNAs and participate in host-virus interactions. Conducting a meta-analysis of available transcriptome data, we identified and characterised eight ADAR transcripts in Chlamys farreri, a farmed marine scallop susceptible to Acute viral necrosis virus (AVNV) infections and mortality outbreaks. Accordingly, we identified six ADAR genes in the Zhikong scallop genome, revised previous gene annotations, and traced alternative splicing variants. In detail, each ADAR gene encodes a unique combination of functional domains, always including the Adenosine deaminase domain, RNA binding domains and, in one case, two copies of a Z-DNA binding domain. After phylogenetic analysis, five C. farreri ADARs clustered in the ADAR1 clade along with sequences from diverse animal phyla. Gene expression analysis indicated CF051320 as the most expressed ADAR, especially in the eye and male gonad. The other four ADAR1 genes and one ADAR2 gene exhibited variable expression levels, with CF105370 and CF051320 significantly increasing during early scallop development. ADAR-mediated single-base editing, evaluated across adult C. farreri tissues and developmental stages, was mainly detectable in intergenic regions (83 % and 85 %, respectively). Overall, the expression patterns of the six ADAR genes together with the editing and hyper-editing values computed on scallops RNA-seq samples support the adaptive value of ADAR1-mediated editing, particularly in the pre-settling larval stages.

Elucidation of the Epitranscriptomic RNA Modification Landscape of Chikungunya Virus.

The genomes of positive-sense (+) single-stranded RNA (ssRNA) viruses are believed to be subjected to a wide range of RNA modifications. In this study, we focused on the chikungunya virus (CHIKV) as a model (+) ssRNA virus to study the landscape of viral RNA modification in infected human cells. Among the 32 distinct RNA modifications analysed by mass spectrometry, inosine was found enriched in the genomic CHIKV RNA. However, orthogonal validation by Illumina RNA-seq analyses did not identify any inosine modification along the CHIKV RNA genome. Moreover, CHIKV infection did not alter the expression of ADAR1 isoforms, the enzymes that catalyse the adenosine to inosine conversion. Together, this study highlights the importance of a multidisciplinary approach to assess the presence of RNA modifications in viral RNA genomes.

Meta-Analysis of Keratoconus Transcriptomic Data Revealed Altered RNA Editing Levels Impacting Keratin Genomic Clusters.

Introduction: Keratoconus (KC) is an ocular disorder with a multifactorial origin. Transcriptomic analyses (RNA-seq) revealed deregulations of coding (mRNA) and non-coding RNAs (ncRNAs) in KC, suggesting that mRNA-ncRNA co-regulations can promote the onset of KC. The present study investigates the modulation of RNA editing mediated by the adenosine deaminase acting on dsRNA (ADAR) enzyme in KC. Materials: The level of ADAR-mediated RNA editing in KC and healthy corneas were determined by two indexes in two different sequencing datasets. REDIportal was used to localize known editing sites, whereas new putative sites were de novo identified in the most extended dataset only and their possible impact was evaluated. Western Blot analysis was used to measure the level of ADAR1 in the cornea from independent samples. Results: KC was characterized by a statistically significant lower RNA-editing level compared to controls, resulting in a lower editing frequency, and less edited bases. The distribution of the editing sites along the human genome showed considerable differences between groups, particularly relevant in the chromosome 12 regions encoding for Keratin type II cluster. A total of 32 recoding sites were characterized, 17 representing novel sites. JUP, KRT17, KRT76, and KRT79 were edited with higher frequencies in KC than in controls, whereas BLCAP, COG3, KRT1, KRT75, and RRNAD1 were less edited. Both gene expression and protein levels of ADAR1 appeared not regulated between diseased and controls. Conclusions: Our findings demonstrated an altered RNA-editing in KC possibly linked to the peculiar cellular conditions. The functional implications should be further investigated.

ADAR-Editing during Ostreid Herpesvirus 1 Infection in Crassostrea gigas: Facts and Limitations.

Ostreid herpesvirus-1 (OsHV-1) RNAs are enzymatically modified by A-to-I conversions during the infection of Crassostrea gigas. The increase of ADAR1 expression and hyper-editing activity parallel to OsHV-1 RNAs suggests a functional connection between dsRNA editing and antiviral responses. We analyzed 87 RNA-seq data sets from immuno-primed, resistant, and susceptible oysters exposed to OsHV-1 to compare the ADAR hyper-editing levels on host and viral transcripts and trace hyper-editing on the oyster genes. Host RNAs were more hyper-edited than viral RNAs, despite the increased editing of viral RNAs in late infection phases. A set of genes, representing ∼0.5% of the oyster transcriptome and including several tripartite motif-containing sequences, were constantly hyper-edited. Conversely, we identified genes involved in antiviral response, miRNA maturation, and epigenetic regulation that were hyper-edited in specific conditions only. Despite technical and biological bottlenecks that hamper the understanding of the bivalve "RNA editome," available tools and technologies can be adapted to bivalve mollusks. IMPORTANCE Ostreid herpesvirus-1 (OsHV-1) is a harmful pathogen of bivalve species, such as oysters. However, knowledge is lacking about host-virus interactions at the molecular level, hampering the possibility of a correct management of viral outbreaks and related massive mortalities. Notably, OsHV-1 transcripts are massively modified by host RNA editing enzyme during infection, resulting in multiple A-to-I variations along RNAs assuming double-strand conformations. The impact of these modifications on host transcripts is, however, not completely clear. Analyzing RNA-seq data of oysters infected with OsHV-1, we revealed that ∼0.5% of the oyster transcriptome is always enzymatically modified by ADAR, whereas genes involved in antiviral response, miRNA maturation, and epigenetic regulation were hyper-edited in specific conditions only. Despite our results, relevant technical bottlenecks impair an accurate quantification of RNA editing events, making necessary an approach specifically dedicated to the progressive understanding of oyster "RNA editome."

A glimpse on metazoan ZNFX1 helicases, ancient players of antiviral innate immunity.

The human zinc finger NFX1-type containing 1 (ZNFX1) is an interferon-stimulated protein associated to the outer mitochondrial membrane, able to bind dsRNAs and interact with MAVS proteins, promoting type I IFN response in the early stage of viral infection. An N-terminal Armadillo (ARM)-type fold and a large helicase core (P-loop) and zinc fingers confer RNA-binding and ATPase activities to ZNFX1. We studied the phylogenetic distribution of metazoan ZNFX1s, ZNFX1 gene expression trends and genomic and protein signatures during viral infection of invertebrates. Based on 221 ZNFX1 sequences, we obtained a polyphyletic tree with a taxonomy-consistent branching at the phylum-level only. In metazoan genomes, ZNFX1 genes were found either in single copy, with up to some tens of exons in vertebrates, or in multiple copies, with one or a few exons and one of them sometimes encompassing most of the coding sequence, in invertebrates like sponges, sea urchins and mollusks. Structural analyses of selected ZNFX1 proteins showed high conservation of the helicase region (P-loop), an overall conserved region and domain architecture, an ARM-fold mostly traceable, and the presence of intrinsically disordered regions of varying length and position. The remarkable over-expression of ZNFX1 in bivalve and gastropod mollusks infected with dsDNA viruses underscores the antiviral role of ZNFX1, whereas nothing similar was found in virus-infected nematodes and corals. Whether the functional diversification reported in the C. elegans ZNFX1 occurs in other metazoan proteins remains to be established.

Viral Decoys: The Only Two Herpesviruses Infecting Invertebrates Evolved Different Transcriptional Strategies to Deflect Post-Transcriptional Editing.

The highly versatile group of Herpesviruses cause disease in a wide range of hosts. In invertebrates, only two herpesviruses are known: the malacoherpesviruses HaHV-1 and OsHV-1 infecting gastropods and bivalves, respectively. To understand viral transcript architecture and diversity we first reconstructed full-length viral genomes of HaHV-1 infecting Haliotis diversicolor supertexta and OsHV-1 infecting Scapharca broughtonii by DNA-seq. We then used RNA-seq over the time-course of experimental infections to establish viral transcriptional dynamics, followed by PacBio long-read sequencing of full-length transcripts to untangle viral transcript architectures at two selected time points. Despite similarities in genome structure, in the number of genes and in the diverse transcriptomic architectures, we measured a ten-fold higher transcript variability in HaHV-1, with more extended antisense gene transcription. Transcriptional dynamics also appeared different, both in timing and expression trends. Both viruses were heavily affected by post-transcriptional modifications performed by ADAR1 affecting sense-antisense gene pairs forming dsRNAs. However, OsHV-1 concentrated these modifications in a few genomic hotspots, whereas HaHV-1 diluted ADAR1 impact by elongated and polycistronic transcripts distributed over its whole genome. These transcriptional strategies might thus provide alternative potential roles for sense-antisense transcription in viral transcriptomes to evade the host's immune response in different virus-host combinations.

NMR Metabolite Profiles of the Bivalve Mollusc Mytilus galloprovincialis Before and After Immune Stimulation With Vibrio splendidus.

The hemolymph metabolome of Mytilus galloprovincialis injected with live Vibrio splendidus bacteria was analyzed by 1H-NMR spectrometry. Changes in spectral hemolymph profiles were already detected after mussel acclimation (3 days at 18 or 25 °C). A significant decrease of succinic acid was accompanied by an increase of most free amino acids, mytilitol, and, to a smaller degree, osmolytes. These metabolic changes are consistent with effective osmoregulation, and the restart of aerobic respiration after the functional anaerobiosis occurred during transport. The injection of Vibrio splendidus in mussels acclimated at 18°C caused a significant decrease of several amino acids, sugars, and unassigned chemical species, more pronounced at 24 than at 12 h postinjection. Correlation heatmaps indicated dynamic metabolic adjustments and the relevance of protein turnover in maintaining the homeostasis during the response to stressful stimuli. This study confirms NMR-based metabolomics as a feasible analytical approach complementary to other omics techniques in the investigation of the functional mussel responses to environmental challenges.

Digging into bivalve miRNAomes: between conservation and innovation.

Bivalves are a diverse mollusc group of economic and ecological importance. An evident resilience to pollution, parasites and extreme environments makes some bivalve species important models for studying adaptation and immunity. Despite substantial progress in sequencing projects of bivalves, information on non-coding genes and gene-regulatory aspects is still lacking. Here, we review the current repertoire of bivalve microRNAs (miRNAs), important regulators of gene expression in Metazoa. We exploited available short non-coding RNA (sncRNA) data for Pinctada martensii, Crassostrea gigas, Corbicula fluminea, Tegillarca granosa and Ruditapes philippinarum, and we produced new sncRNA data for two additional bivalves, the Mediterranean mussel Mytilus galloprovincialis and the blood clam Scapharca broughtonii. We found substantial heterogeneity and incorrect annotations of miRNAs; hence, we reannotated conserved miRNA families using recently established criteria for bona fide microRNA annotation. We found 106 miRNA families missing in the previously published bivalve datasets and 89 and 87 miRNA complements were identified in the two additional species. The overall results provide a homogeneous and evolutionarily consistent picture of miRNAs in bivalves and enable future comparative studies. The identification of two bivalve-specific miRNA families sheds further light on the complexity of transcription and its regulation in bivalve molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.