Systematic identification of A-to-I RNA editing in zebrafish development and adult organs

miR-430 regulates zygotic mRNA during zebrafish embryogenesis

BACKGROUND

Early embryonic developmental programs are guided by the coordinated interplay between maternally inherited and zygotically manufactured RNAs and proteins. Although these processes happen concomitantly and affecting gene function during this period is bound to affect both pools of mRNAs, it has been challenging to study their expression dynamics separately.

RESULTS

By employing SLAM-seq, a nascent mRNA labeling transcriptomic approach, in a developmental time series we observe that over half of the early zebrafish embryo transcriptome consists of maternal-zygotic genes, emphasizing their pivotal role in early embryogenesis. We provide an hourly resolution of de novo transcriptional activation events and follow nascent mRNA trajectories, finding that most de novo transcriptional events are stable throughout this period. Additionally, by blocking microRNA-430 function, a key post transcriptional regulator during zebrafish embryogenesis, we directly show that it destabilizes hundreds of de novo transcribed mRNAs from pure zygotic as well as maternal-zygotic genes. This unveils a novel miR-430 function during embryogenesis, fine-tuning zygotic gene expression.

CONCLUSION

These insights into zebrafish early embryo transcriptome dynamics emphasize the significance of post-transcriptional regulators in zygotic genome activation. The findings pave the way for future investigations into the coordinated interplay between transcriptional and post-transcriptional landscapes required for the establishment of animal cell identities and functions.

Tissue-specific temperature dependence of RNA editing levels in zebrafish

BACKGROUND

RNA editing by adenosine deaminase acting on RNA (ADAR) occurs in all metazoans and fulfils several functions. Here, we examined effects of acclimation temperature (27 °C, 18 °C,13 °C) on editing patterns in six tissues of zebrafish (Danio rerio).

RESULTS

Sites and total amounts of editing differed among tissues. Brain showed the highest levels, followed by gill and skin. In these highly edited tissues, decreases in temperatures led to large increases in total amounts of editing and changes in specific edited sites. Gene ontology analysis showed both similarities (e.g., endoplasmic reticulum stress response) and differences in editing among tissues. The majority of edited sites were in transcripts of transposable elements and the 3'UTR regions of protein coding genes. By experimental validation, translation efficiency was directly related to extent of editing of the 3'UTR region of an mRNA.

CONCLUSIONS

RNA editing increases 3'UTR polymorphism and affects efficiency of translation. Such editing may lead to temperature-adaptive changes in the proteome through altering relative amounts of synthesis of different proteins.

High-level RNA editing diversifies the coleoid cephalopod brain proteome

Coleoid cephalopods (octopus, squid and cuttlefish) have unusually complex nervous systems. The coleoid nervous system is also the only one currently known to recode the majority of expressed proteins through A-to-I RNA editing. The deamination of adenosine by adenosine deaminase acting on RNA (ADAR) enzymes produces inosine, which is interpreted as guanosine during translation. If this occurs in an open reading frame, which is the case for tens of thousands of editing sites in coleoids, it can recode the encoded protein. Here, we describe recent findings aimed at deciphering the mechanisms underlying high-level recoding and its adaptive potential. We describe the complement of ADAR enzymes in cephalopods, including a recently discovered novel domain in sqADAR1. We further summarize current evidence supporting an adaptive role of high-level RNA recoding in coleoids, and review recent studies showing that a large proportion of recoding sites is temperature-sensitive. Despite these new findings, the mechanisms governing the high level of RNA recoding in coleoid cephalopods remain poorly understood. Recent advances using genome editing in squid may provide useful tools to further study A-to-I RNA editing in these animals.

Deep transcriptome profiling reveals limited conservation of A-to-I RNA editing in Xenopus

BACKGROUND

Xenopus has served as a valuable model system for biomedical research over the past decades. Notably, ADAR was first detected in frog oocytes and embryos as an activity that unwinds RNA duplexes. However, the scope of A-to-I RNA editing by the ADAR enzymes in Xenopus remains underexplored.

RESULTS

Here, we identify millions of editing events in Xenopus with high accuracy and systematically map the editome across developmental stages, adult organs, and species. We report diverse spatiotemporal patterns of editing with deamination activity highest in early embryogenesis before zygotic genome activation and in the ovary. Strikingly, editing events are poorly conserved across different Xenopus species. Even sites that are detected in both X. laevis and X. tropicalis show largely divergent editing levels or developmental profiles. In protein-coding regions, only a small subset of sites that are found mostly in the brain are well conserved between frogs and mammals.

CONCLUSIONS

Collectively, our work provides fresh insights into ADAR activity in vertebrates and suggest that species-specific editing may play a role in each animal's unique physiology or environmental adaptation.

Identification of Bona Fide RNA Editing Sites: History, Challenges, and Opportunities

Adenosine-to-inosine (A-to-I) RNA editing, catalyzed by the adenosine deaminase acting on the RNA (ADAR) family of enzymes of which there are three members (ADAR1, ADAR2, and ADAR3), is a major gene regulatory mechanism that diversifies the transcriptome. It is widespread in many metazoans, including humans. As inosine is interpreted by cellular machineries mainly as guanosine, A-to-I editing effectively gives A-to-G nucleotide changes. Depending on its location, an editing event can generate new protein isoforms or influence other RNA processing pathways. Researchers have found that ADAR-mediated editing performs diverse functions. For example, it enables living organisms such as cephalopods to adapt rapidly to fluctuating environmental conditions such as water temperature. In development, the loss of ADAR1 is embryonically lethal partly because endogenous double-stranded RNAs (dsRNAs) are no longer marked by inosines, which signal "self", and thus cause the melanoma differentiation-associated protein 5 (MDA5) sensor to trigger a deleterious interferon response. Hence, ADAR1 plays a key role in preventing aberrant activation of the innate immune system. Furthermore, ADAR enzymes have been implicated in myriad human diseases. Intriguingly, some cancer cells are known to exploit ADAR1 activity to dodge immune responses. However, the exact identities of immunogenic RNAs in different biological contexts have remained elusive. Consequently, there is tremendous interest in identifying inosine-containing RNAs in the cell.The identification of A-to-I RNA editing sites is dependent on the sequencing of nucleic acids. Technological and algorithmic advancements over the past decades have revolutionized the way editing events are detected. At the beginning, the discovery of editing sites relies on Sanger sequencing, a first-generation technology. Both RNA, which is reverse transcribed into complementary DNA (cDNA), and genomic DNA (gDNA) from the same source are analyzed. After sequence alignment, one would require an adenosine to be present in the genome but a guanosine to be detected in the RNA sample for a position to be declared as an editing site. However, an issue with Sanger sequencing is its low throughput. Subsequently, Illumina sequencing, a second-generation technology, was invented. By permitting the simultaneous interrogation of millions of molecules, it enables many editing sites to be identified rapidly. However, a key challenge is that the Illumina platform produces short sequencing reads that can be difficult to map accurately. To tackle the challenge, we and others developed computational workflows with a series of filters to discard sites that are likely to be false positives. When Illumina sequencing data sets are properly analyzed, A-to-G variants should emerge as the most dominant mismatch type. Moreover, the quantitative nature of the data allows us to build a comprehensive atlas of editing-level measurements across different biological contexts, providing deep insights into the spatiotemporal dynamics of RNA editing. However, difficulties remain in identifying true A-to-I editing sites in short protein-coding exons or in organisms and diseases where DNA mutations and genomic polymorphisms are prevalent and mostly unknown. Nanopore sequencing, a third-generation technology, promises to address the difficulties, as it allows native RNAs to be sequenced without conversion to cDNA, preserving base modifications that can be directly detected through machine learning. We recently demonstrated that nanopore sequencing could be used to identify A-to-I editing sites in native RNA directly. Although further work is needed to enhance the detection accuracy in single molecules from fewer cells, the nanopore technology holds the potential to revolutionize epitranscriptomic studies.

Genetic perturbation of AMP biosynthesis extends lifespan and restores metabolic health in a naturally short-lived vertebrate

During aging, the loss of metabolic homeostasis drives a myriad of pathologies. A central regulator of cellular energy, the AMP-activated protein kinase (AMPK), orchestrates organismal metabolism. However, direct genetic manipulations of the AMPK complex in mice have, so far, produced detrimental phenotypes. Here, as an alternative approach, we alter energy homeostasis by manipulating the upstream nucleotide pool. Using the turquoise killifish, we mutate APRT, a key enzyme in AMP biosynthesis, and extend the lifespan of heterozygous males. Next, we apply an integrated omics approach to show that metabolic functions are rejuvenated in old mutants, which also display a fasting-like metabolic profile and resistance to high-fat diet. At the cellular level, heterozygous cells exhibit enhanced nutrient sensitivity, reduced ATP levels, and AMPK activation. Finally, lifelong intermittent fasting abolishes the longevity benefits. Our findings suggest that perturbing AMP biosynthesis may modulate vertebrate lifespan and propose APRT as a promising target for promoting metabolic health.

Squid express conserved ADAR orthologs that possess novel features

The coleoid cephalopods display unusually extensive mRNA recoding by adenosine deamination, yet the underlying mechanisms are not well understood. Because the adenosine deaminases that act on RNA (ADAR) enzymes catalyze this form of RNA editing, the structure and function of the cephalopod orthologs may provide clues. Recent genome sequencing projects have provided blueprints for the full complement of coleoid cephalopod ADARs. Previous results from our laboratory have shown that squid express an ADAR2 homolog, with two splice variants named sqADAR2a and sqADAR2b and that these messages are extensively edited. Based on octopus and squid genomes, transcriptomes, and cDNA cloning, we discovered that two additional ADAR homologs are expressed in coleoids. The first is orthologous to vertebrate ADAR1. Unlike other ADAR1s, however, it contains a novel N-terminal domain of 641 aa that is predicted to be disordered, contains 67 phosphorylation motifs, and has an amino acid composition that is unusually high in serines and basic amino acids. mRNAs encoding sqADAR1 are themselves extensively edited. A third ADAR-like enzyme, sqADAR/D-like, which is not orthologous to any of the vertebrate isoforms, is also present. Messages encoding sqADAR/D-like are not edited. Studies using recombinant sqADARs suggest that only sqADAR1 and sqADAR2 are active adenosine deaminases, both on perfect duplex dsRNA and on a squid potassium channel mRNA substrate known to be edited in vivo. sqADAR/D-like shows no activity on these substrates. Overall, these results reveal some unique features in sqADARs that may contribute to the high-level RNA recoding observed in cephalopods.

On the origin and evolution of RNA editing in metazoans

Extensive adenosine-to-inosine (A-to-I) editing of nuclear-transcribed mRNAs is the hallmark of metazoan transcriptional regulation. Here, by profiling the RNA editomes of 22 species that cover major groups of Holozoa, we provide substantial evidence supporting A-to-I mRNA editing as a regulatory innovation originating in the last common ancestor of extant metazoans. This ancient biochemistry process is preserved in most extant metazoan phyla and primarily targets endogenous double-stranded RNA (dsRNA) formed by evolutionarily young repeats. We also find intermolecular pairing of sense-antisense transcripts as an important mechanism for forming dsRNA substrates for A-to-I editing in some but not all lineages. Likewise, recoding editing is rarely shared across lineages but preferentially targets genes involved in neural and cytoskeleton systems in bilaterians. We conclude that metazoan A-to-I editing might first emerge as a safeguard mechanism against repeat-derived dsRNA and was later co-opted into diverse biological processes due to its mutagenic nature.

Extensive Recoding of the Neural Proteome in Cephalopods by RNA Editing

The coleoid cephalopods have the largest brains, and display the most complex behaviors, of all invertebrates. The molecular and cellular mechanisms that underlie these remarkable advancements remain largely unexplored. Early molecular cloning studies of squid ion channel transcripts uncovered an unusually large number of A→I RNA editing sites that recoded codons. Further cloning of other neural transcripts showed a similar pattern. The advent of deep-sequencing technologies and the associated bioinformatics allowed the mapping of RNA editing events across the entire neural transcriptomes of various cephalopods. The results were remarkable: They contained orders of magnitude more recoding editing sites than any other taxon. Although RNA editing sites are abundant in most multicellular metazoans, they rarely recode. In cephalopods, the majority of neural transcripts are recoded. Recent studies have focused on whether these events are adaptive, as well as other noncanonical aspects of cephalopod RNA editing.

Recurrent RNA edits in human preimplantation potentially enhance maternal mRNA clearance

Posttranscriptional modification plays an important role in key embryonic processes. Adenosine-to-inosine RNA editing, a common example of such modifications, is widespread in human adult tissues and has various functional impacts and clinical consequences. However, whether it persists in a consistent pattern in most human embryos, and whether it supports embryonic development, are poorly understood. To address this problem, we compiled the largest human embryonic editome from 2,071 transcriptomes and identified thousands of recurrent embryonic edits (>=50% chances of occurring in a given stage) for each early developmental stage. We found that these recurrent edits prefer exons consistently across stages, tend to target genes related to DNA replication, and undergo organized loss in abnormal embryos and embryos from elder mothers. In particular, these recurrent edits are likely to enhance maternal mRNA clearance, a possible mechanism of which could be introducing more microRNA binding sites to the 3'-untranslated regions of clearance targets. This study suggests a potentially important, if not indispensable, role of RNA editing in key human embryonic processes such as maternal mRNA clearance; the identified editome can aid further investigations.

The role of post-transcriptional modifications during development

While the existence of post-transcriptional modifications of RNA nucleotides has been known for decades, in most RNA species the exact positions of these modifications and their physiological function have been elusive until recently. Technological advances, such as high-throughput next-generation sequencing (NGS) methods and nanopore-based mapping technologies, have made it possible to map the position of these modifications with single nucleotide accuracy, and genetic screens have uncovered the “writer”, “reader” and “eraser” proteins that help to install, interpret and remove such modifications, respectively. These discoveries led to intensive research programmes with the aim of uncovering the roles of these modifications during diverse biological processes. In this review, we assess novel discoveries related to the role of post-transcriptional modifications during animal development, highlighting how these discoveries can affect multiple aspects of development from fertilization to differentiation in many species.

Granulosa Cell Specific Loss of Adar in Mice Delays Ovulation, Oocyte Maturation and Leads to Infertility

Adenosine deaminases acting on RNA-(ADAR) comprise one family of RNA editing enzymes that specifically catalyze adenosine to inosine (A-to-I) editing. A granulosa cell (GC) specific Adar depleted mouse model [Adar flox/flox:Cyp19a1-Cre/+ (gcAdarKO)] was used to evaluate the role of ADAR1 during the periovulatory period. Loss of Adar in GCs led to failure to ovulate at 16 h post-hCG, delayed oocyte germinal vesicle breakdown and severe infertility. RNAseq analysis of GC collected from gcAdarKO and littermate control mice at 0 and 4 h post-hCG following a super-ovulatory dose of eCG (48 h), revealed minimal differences after eCG treatment alone (0 h), consistent with normal folliculogenesis observed histologically and uterine estrogenic responses. In contrast, 300 differential expressed genes (DEGs; >1.5-fold change and FDRP < 0.1) were altered at 4 h post-hCG. Ingenuity pathway analysis identified many downstream targets of estrogen and progesterone pathways, while multiple genes involved in inflammatory responses were upregulated in the gcAdarKO GCs. Temporal expression analysis of GCs at 0, 4, 8, and 12 h post-hCG of Ifi44, Ifit1, Ifit3b, and Oas1g and Ovgp1 confirmed upregulation of these inflammatory and interferon genes and downregulation of Ovgp1 a glycoprotein involved in oocyte zona pellucida stability. Thus, loss of ADAR1 in GCs leads to increased expression of inflammatory and interferon response genes which are temporally linked to ovulation failure, alterations in oocyte developmental progression and infertility.

Proteomic Analysis of Zebrafish Protein Recoding via mRNA Editing by ADAR Enzymes

RNA editing by adenosine deaminases of the ADAR family can lead to protein recoding, since inosine formed from adenosine in mRNA is complementary to cytosine; the resulting codon editing might introduce amino acid substitutions into translated proteins. Proteome recoding can have functional consequences which have been described in many animals including humans. Using protein recoding database derived from publicly available transcriptome data, we identified for the first time the recoding sites in the zebrafish shotgun proteomes. Out of more than a hundred predicted recoding events, ten substitutions were found in six used datasets. Seven of them were in the AMPA glutamate receptor subunits, whose recoding has been well described, and are conserved among vertebrates. Three sites were specific for zebrafish proteins and were found in the transmembrane receptors astrotactin 1 and neuregulin 3b (proteins involved in the neuronal adhesion and signaling) and in the rims2b gene product (presynaptic membrane protein participating in the neurotransmitter release), respectively. Further studies are needed to elucidate the role of recoding of the said three proteins in the zebrafish.

Genome-wide detection of RNA editing events during the hair follicles cycle of Tianzhu white yak

BACKGROUND

The hair coat is available for the yak to live in the harsh environment of the plateau. Besides, improving the hair production of yak is necessary for its textile industry development. Hair grows from hair follicles (HFs). The HFs undergo periodic growth after birth and are regulated by the complex gene regulatory network. However, the molecular mechanism of HFs regeneration in the Tianzhu white yak remains unclear. RNA editing is a post-transcriptional mechanism that regulates gene expression and produces new transcripts. Hence, we investigated the influence of the A-to-I RNA editing events on the HFs cycle of the Tianzhu white yak.

RESULTS

We finally identified 54,707 adenosine-to-inosine (A-to-I) RNA editing sites (RESs) from RNA sequencing data of the HFs cycle in the Tianzhu white yak. Annotation results showed RESs caused missense amino acid changes in 7 known genes. And 202 A-to-I editing sites altered 23 target genes of 140 microRNAs. A total of 1,722 differential RESs were identified during the HFs cycle of Tianzhu white yak. GO and KEGG enrichment analysis revealed several signaling pathways and GO terms involved skin development, hair growth, and HFs cycle. Such as genes with differential RNA editing levels were significantly enriched in the peroxisome, metabolic pathways, Notch signaling pathway, and PPAR signaling pathway. Besides, the editing sites in HFs development-related genes FAS, APCDD1, WWOX, MPZL3, RUNX1, KANK2, DCN, DSC2, LEPR, HEPHL1, and PTK2B were suggested as the potential RESs involving HFs development.

CONCLUSION

This study investigated the global A-to-I RNA editing events during the HFs cycle of yak skin tissue and expanded the knowledge of A-to-I RNA editing on the HFs cycle. Furthermore, this study revealed that RNA editing-influenced genes may regulate the HFs cycle by participating in the HFs development-related pathways. The findings might provide new insight into the regulation of RNA editing in hair growth.

Adar-mediated A-to-I editing is required for embryonic patterning and innate immune response regulation in zebrafish

Adenosine deaminases (ADARs) catalyze the deamination of adenosine to inosine, also known as A-to-I editing, in RNA. Although A-to-I editing occurs widely across animals and is well studied, new biological roles are still being discovered. Here, we study the role of A-to-I editing in early zebrafish development. We demonstrate that Adar, the zebrafish orthologue of mammalian ADAR1, is essential for establishing the antero-posterior and dorso-ventral axes and patterning. Genome-wide editing discovery reveals pervasive editing in maternal and the earliest zygotic transcripts, the majority of which occurred in the 3'-UTR. Interestingly, transcripts implicated in gastrulation as well as dorso-ventral and antero-posterior patterning are found to contain multiple editing sites. Adar knockdown or overexpression affect gene expression by 12 hpf. Analysis of adar-/- zygotic mutants further reveals that the previously described role of Adar in mammals in regulating the innate immune response is conserved in zebrafish. Our study therefore establishes distinct maternal and zygotic functions of RNA editing by Adar in embryonic patterning along the zebrafish antero-posterior and dorso-ventral axes, and in the regulation of the innate immune response, respectively.

Direct identification of A-to-I editing sites with nanopore native RNA sequencing

Inosine is a prevalent RNA modification in animals and is formed when an adenosine is deaminated by the ADAR family of enzymes. Traditionally, inosines are identified indirectly as variants from Illumina RNA-sequencing data because they are interpreted as guanosines by cellular machineries. However, this indirect method performs poorly in protein-coding regions where exons are typically short, in non-model organisms with sparsely annotated single-nucleotide polymorphisms, or in disease contexts where unknown DNA mutations are pervasive. Here, we show that Oxford Nanopore direct RNA sequencing can be used to identify inosine-containing sites in native transcriptomes with high accuracy. We trained convolutional neural network models to distinguish inosine from adenosine and guanosine, and to estimate the modification rate at each editing site. Furthermore, we demonstrated their utility on the transcriptomes of human, mouse and Xenopus. Our approach expands the toolkit for studying adenosine-to-inosine editing and can be further extended to investigate other RNA modifications.

Genome-wide characterization of RNA editing highlights roles of high editing events of glutamatergic synapse during mouse retinal development

Adenosine-to-inosine (A-to-I) RNA editing leads to functional change of neurotransmitter receptor which is essential for neurotransmission and normal neuronal development. As a highly accessible part of central nervous system, retina has been extensively studied, however, it remains largely unknown how RNA editing regulates its development. Here, a genome-wide screening of high-confidence RNA editing events were performed to decipher the dynamic transcriptome regulation by RNA editing during mouse retinal development. 2000 high-confidence editing sites across eight developmental stages of retina were called. Three unique patterns (RNA-editing^high pattern, RNA-editing^medium pattern and RNA-editing^low pattern) were identified by clustering these editing sites based on their editing level during retinal development. Editing events from RNA-editing^high pattern were significantly associated with glutamate receptors and regulated synaptic transmission. Interestingly, most non-synonymous high-editing sites were mapped to ion channel genes of glutamatergic synapse which were associated with neurotransmission by controlling ion channel permeability and affecting exocytosis. Meanwhile, these non-synonymous editing sites were evolutionarily conserved and exhibited a consistently increasing editing levels between mouse and human retinal development. Single-cell RNA-seq data analysis revealed that RNA editing events prefer to occur in two main cell types including bipolar and amacrine cells. Genes with non-synonymous high-editing sites were enriched in both bipolar cells and retina ganglion cells, which may mediate retina ganglion cell differentiation by altering channel ion permeability. Together, our results provide novel insights into mechanism of post-transcriptional regulation during retinal development and help to develop novel RNA editing-guided therapeutic strategies for retinal disorders.