APOBEC3A is implicated in a novel class of G-to-A mRNA editing in WT1 transcripts

The role of distinct APOBEC/ADAR mRNA levels in mutational signatures linked to aging and ultraviolet radiation

The APOBEC/AID family is known for its mutator activity, and recent evidence also supports the potential impact of ADARs. Furthermore, the mutator impacts of APOBEC/ADAR mutations have not yet been investigated. Assessment of pancancer TCGA exomes identified enriched somatic variants among exomes with nonsynonymous APOBEC1, APOBEC3B, APOBEC3C, ADAR, and ADARB1 mutations, compared to exomes with synonymous ones. Principal component (PC) analysis reduced the number of potential players to eight in cancer exomes/genomes, and to five in cancer types. Multivariate regression analysis was used to assess the impact of the PCs on each COSMIC mutational signature among pancancer exomes/genomes and particular cancers, identifying several novel links, including SBS17b, SBS18, and ID7 mainly determined by APOBEC1 mRNA levels; SBS40, ID1, and ID2 by age; SBS3 and SBS16 by APOBEC3A/APOBEC3B mRNA levels; ID5 and DBS9 by DNA repair/replication (DRR) defects; and SBS7a-d, SBS38, ID4, ID8, ID13, and DBS1 by ultraviolet (UV) radiation/ADARB1 mRNA levels. APOBEC/ADAR mutations appeared to potentiate the impact of DRR defects on several mutational signatures, and some factors seemed to inversely affect certain signatures. These findings potentially implicate certain APOBEC/ADAR mutations/mRNA levels in distinct mutational signatures, particularly APOBEC1 mRNA levels in aging-related signatures and ADARB1 mRNA levels in UV radiation-related signatures.

The cytidine deaminase APOBEC3A regulates nucleolar function to promote cell growth and ribosome biogenesis

Cancer initiates as a consequence of genomic mutations and its subsequent progression relies in part on increased production of ribosomes to maintain high levels of protein synthesis for unchecked cell growth. Recently, cytidine deaminases have been uncovered as sources of mutagenesis in cancer. In an attempt to form a connection between these 2 cancer driving processes, we interrogated the cytidine deaminase family of proteins for potential roles in human ribosome biogenesis. We identified and validated APOBEC3A and APOBEC4 as novel ribosome biogenesis factors through our laboratory's established screening platform for the discovery of regulators of nucleolar function in MCF10A cells. Through siRNA depletion experiments, we highlight APOBEC3A's requirement in making ribosomes and specific role within the processing and maturation steps that form the large subunit 5.8S and 28S ribosomal (r)RNAs. We demonstrate that a subset of APOBEC3A resides within the nucleolus and associates with critical ribosome biogenesis factors. Mechanistic insight was revealed by transient overexpression of both wild-type and a catalytically dead mutated APOBEC3A, which both increase cell growth and protein synthesis. Through an innovative nuclear RNA sequencing methodology, we identify only modest predicted APOBEC3A C-to-U target sites on the pre-rRNA and pre-mRNAs. Our work reveals a potential direct role for APOBEC3A in ribosome biogenesis likely independent of its editing function. More broadly, we found an additional function of APOBEC3A in cancer pathology through its function in ribosome biogenesis, expanding its relevance as a target for cancer therapeutics.

RNA base editors: The emerging approach of RNA therapeutics

RNA-based therapeutics offer a flexible and reversible approach for treating genetic disorders, such as antisense oligonucleotides, RNA interference, aptamers, mRNA vaccines, and RNA editing. In recent years, significant advancements have been made in RNA base editing to correct disease-relevant point mutations. These achievements have significantly influenced the fields of biotechnology, biomedical research and therapeutics development. In this article, we provide a comprehensive overview of the design and performance of contemporary RNA base editors, including A-to-I, C-to-U, A-to-m6A, and U-to-Ψ. We compare recent innovative developments and highlight their applications in disease-relevant contexts. Lastly, we discuss the limitations and future prospects of utilizing RNA base editing for therapeutic purposes. This article is categorized under: RNA Processing > RNA Editing and Modification RNA in Disease and Development > RNA in Development.

Programmable RNA base editing via targeted modifications

Clustered regularly interspaced short palindromic repeats (CRISPR)-based genome editors are powerful tools in biology and hold great promise for the treatment of human diseases. Advanced DNA base editing tools, such as cytosine base editor and adenine base editor, have been developed to correct permanent mistakes in genetic material. However, undesired off-target edits would also be permanent, which poses a considerable risk for therapeutics. Alternatively, base editing at the RNA level is capable of correcting disease-causing mutations but does not lead to lasting genotoxic effects. RNA base editors offer temporary and reversible therapies and have been catching on in recent years. Here, we summarize some emerging RNA editors based on A-to-inosine, C-to-U and U-to-pseudouridine changes. We review the programmable RNA-targeting systems as well as modification enzyme-based effector proteins and highlight recent technological breakthroughs. Finally, we compare editing tools, discuss limitations and opportunities, and provide insights for the future directions of RNA base editing.

Detection of Developmental Asexual Stage-Specific RNA Editing Events in Plasmodium falciparum 3D7 Malaria Parasite

Transcriptional variation has been studied but post-transcriptional modification due to RNA editing has not been investigated in Plasmodium. We investigated developmental stage-specific RNA editing in selected genes in Plasmodium falciparum 3D7. We detected extensive amination- and deamination-type RNA editing at 8, 16, 24, 32, 40, and 46 h in tightly synchronized Plasmodium. Most of the editing events were observed in 8 and 16 h ring-stage parasites. Extensive A-to-G deamination-type editing was detected more during the 16 h ring stage (25%) than the 8 h ring stage (20%). Extensive U-to-C amination-type editing was detected more during the 16 h ring stage (31%) than the 8 h ring stage (22%). In 28S, rRNA editing converted the loop structure to the stem structure. The hemoglobin binding activity of PF3D7_0216900 was also altered due to RNA editing. Among the expressed 28S rRNA genes, PF3D7_0532000 and PF3D7_0726000 expression was higher. Increased amounts of the transcripts of these two genes were found, particularly PF3D7_0726000 in the ring stage and PF3D7_0532000 in the trophozoite and schizont stages. Adenosine deaminase (ADA) expression did not correlate with the editing level. This first experimental report of RNA editing will help to identify the editing machinery that might be useful for antimalarial drug discovery and malaria control.

Vaccinia-related kinase 2 variants differentially affect breast cancer growth by regulating kinase activity

Genetic information is transcribed from genomic DNA to mRNA, which is then translated into three-dimensional proteins. mRNAs can undergo various post-transcriptional modifications, including RNA editing that alters mRNA sequences, ultimately affecting protein function. In this study, RNA editing was identified at the 499th base (c.499) of human vaccinia-related kinase 2 (VRK2). This RNA editing changes the amino acid in the catalytic domain of VRK2 from isoleucine (with adenine base) to valine (with guanine base). Isoleucine-containing VRK2 has higher kinase activity than the valine-containing VRK2, which leads to an increase in tumor cell proliferation. Earlier we reported that VRK2 directly interacts with dystrobrevin-binding protein (dysbindin) and results in reducing its stability. Herein, we demonstrate that isoleucine-containing VRK2 decreases the level of dysbindin than valine-containing VRK2. Dysbindin interacts with cyclin D and thereby regulates its expression and function. The reduction in the level of dysbindin by isoleucine-containing VRK2 further enhances the cyclin D expression, resulting in increased tumor growth and reduction in survival rates. It has also been observed that in patient samples, VRK2 level was elevated in breast cancer tissue compared to normal breast tissue. Additionally, the isoleucine form of VRK2 exhibited a greater increase in breast cancer tissue. Therefore, it is concluded that VRK2, especially dependent on the 167th variant amino acid, can be one of the indexes of tumor progression and proliferation.

RNA Editing in Cancer Progression

Coding and noncoding RNA molecules play their roles in ensuring cell function and tissue homeostasis in an ordered and systematic fashion. RNA chemical modifications can occur both at bases and ribose sugar, and, similarly to DNA and histone modifications, can be written, erased, and recognized by the corresponding enzymes, thus modulating RNA activities and fine-tuning gene expression programs. RNA editing is one of the most prevalent and abundant forms of post-transcriptional RNA modification in normal physiological processes. By altering the sequences of mRNAs, it makes them different from the corresponding genomic template. Hence, edited mRNAs can produce protein isoforms that are functionally different from the corresponding genome-encoded variants. Abnormalities in regulatory enzymes and changes in RNA-modification patterns are closely associated with the occurrence and development of various human diseases, including cancer. To date, the roles played by RNA modifications in cancer are gathering increasing interest. In this review, we focus on the role of RNA editing in cancer transformation and provide a new perspective on its impact on tumorigenesis, by regulating cell proliferation, differentiation, invasion, migration, stemness, metabolism, and drug resistance.

Unraveling the Enzyme-Substrate Properties for APOBEC3A-Mediated RNA Editing

The APOBEC3 family of human cytidine deaminases is involved in various cellular processes, including the innate and acquired immune system, mostly through inducing C-to-U in single-stranded DNA and/or RNA mutations. Although recent studies have examined RNA editing by APOBEC3A (A3A), its intracellular target specificity are not fully characterized. To address this gap, we performed in-depth analysis of cellular RNA editing using our recently developed sensitive cell-based fluorescence assay. Our findings demonstrate that A3A and an A3A-loop1-containing APOBEC3B (A3B) chimera are capable of RNA editing. We observed that A3A prefers to edit specific RNA substrates which are not efficiently deaminated by other APOBEC members. The editing efficiency of A3A is influenced by the RNA sequence contexts and distinct stem-loop secondary structures. Based on the identified RNA specificity features, we predicted potential A3A-editing targets in the encoding region of cellular mRNAs and discovered novel RNA transcripts that are extensively edited by A3A. Furthermore, we found a trend of increased synonymous mutations at the sites for more efficient A3A-editing, indicating evolutionary adaptation to the higher editing rate by A3A. Our results shed light on the intracellular RNA editing properties of A3A and provide insights into new RNA targets and potential impact of A3A-mediated RNA editing.

APOBEC3G protects the genome of human cultured cells and mice from radiation-induced damage

Cytosine deaminases AID/APOBEC proteins act as potent nucleic acid editors, playing important roles in innate and adaptive immunity. However, the mutagenic effects of some of these proteins compromise genomic integrity and may promote tumorigenesis. Here, we demonstrate that human APOBEC3G (A3G), in addition to its role in innate immunity, promotes repair of double-strand breaks (DSBs) in vitro and in vivo. Transgenic mice expressing A3G successfully survived lethal irradiation, whereas wild-type controls quickly succumbed to radiation syndrome. Mass spectrometric analyses identified the differential upregulation of a plethora of proteins involved in DSB repair pathways in A3G-expressing cells early following irradiation to facilitate repair. Importantly, we find that A3G not only accelerates DSB repair but also promotes deamination-dependent error-free rejoining. These findings have two implications: (a) strategies aimed at inhibiting A3G may improve the efficacy of genotoxic therapies used to cure malignant tumours; and (b) enhancing A3G activity may reduce acute radiation syndrome in individuals exposed to ionizing radiation.

Molecular mechanism of Wilms' tumor (Wt1) (+/-KTS) variants promoting proliferation and migration of ovarian epithelial cells by bioinformatics analysis

Epithelial ovarian cancer (EOC) is a gynecological disease with the highest mortality. With the lack of understanding of its pathogenesis, no accurate early diagnosis and screening method has been established for EOC. Studies revealed the multi-faceted function of Wilms' tumor (Wt1) genes in cancer, which may be related to the existence of multiple alternative splices. Our results show that Wt1 (+KTS) or Wt1 (-KTS) overexpression can significantly promote the proliferation and migration of human ovarian epithelial cells HOSEpiC, and Wt1 (+KTS) effects were more evident. To explore the Wt1 (+/-KTS) variant mechanism in HOSEpiC proliferation and migration and ovarian cancer (OC) occurrence and development, this study explored the differential regulation of Wt1 (+/-KTS) in HOSEpiC proliferation and migration by transcriptome sequencing. OC-related hub genes were screened by bioinformatics analysis to further explore the differential molecular mechanism of Wt1 (+/-KTS) in the occurrence of OC. Finally, we found that the regulation of Wt1 (+/-KTS) variants on the proliferation and migration of HOSEpiC may act through different genes and signaling pathways and screened out key genes and differentially regulated genes that regulate the malignant transformation of ovarian epithelial cells. The implementation of this study will provide new clues for the early diagnosis and precise treatment of OC.

Epitranscriptome Analysis of Oxidative Stressed Retinal Epithelial Cells Depicted a Possible RNA Editing Landscape of Retinal Degeneration

Oxidative stress represents one of the principal causes of inherited retinal dystrophies, with many related molecular mechanisms still unknown. We investigated the posttranscriptional RNA editing landscape of human retinal pigment epithelium cells (RPE) exposed to the oxidant agent N-retinylidene-N-retinyl ethanolamine (A2E) for 1 h, 2 h, 3 h and 6 h. Using a transcriptomic approach, refined with a specific multialgorithm pipeline, 62,880 already annotated and de novo RNA editing sites within about 3000 genes were identified among all samples. Approximately 19% of these RNA editing sites were found within 3' UTR, including sites common to all time points that were predicted to change the binding capacity of 359 miRNAs towards 9654 target genes. A2E exposure also determined significant gene expression differences in deaminase family ADAR, APOBEC and ADAT members, involved in canonical and tRNA editing events. On GO and KEGG enrichment analyses, genes that showed different RNA editing levels are mainly involved in pathways strongly linked to a possible neovascularization of retinal tissue, with induced apoptosis mediated by the ECM and surface protein altered signaling. Collectively, this work demonstrated dynamic RNA editome profiles in RPE cells for the first time and shed more light on new mechanisms at the basis of retinal degeneration.

A Comprehensive Evolutionary Study of Chloroplast RNA Editing in Gymnosperms: A Novel Type of G-to-A RNA Editing Is Common in Gymnosperms

Although more than 9100 plant plastomes have been sequenced, RNA editing sites of the whole plastome have been experimentally verified in only approximately 21 species, which seriously hampers the comprehensive evolutionary study of chloroplast RNA editing. We investigated the evolutionary pattern of chloroplast RNA editing sites in 19 species from all 13 families of gymnosperms based on a combination of genomic and transcriptomic data. We found that the chloroplast C-to-U RNA editing sites of gymnosperms shared many common characteristics with those of other land plants, but also exhibited many unique characteristics. In contrast to that noted in angiosperms, the density of RNA editing sites in ndh genes was not the highest in the sampled gymnosperms, and both loss and gain events at editing sites occurred frequently during the evolution of gymnosperms. In addition, GC content and plastomic size were positively correlated with the number of chloroplast RNA editing sites in gymnosperms, suggesting that the increase in GC content could provide more materials for RNA editing and facilitate the evolution of RNA editing in land plants or vice versa. Interestingly, novel G-to-A RNA editing events were commonly found in all sampled gymnosperm species, and G-to-A RNA editing exhibits many different characteristics from C-to-U RNA editing in gymnosperms. This study revealed a comprehensive evolutionary scenario for chloroplast RNA editing sites in gymnosperms, and reported that a novel type of G-to-A RNA editing is prevalent in gymnosperms.

Intra-host evolution during SARS-CoV-2 prolonged infection

Long-term infection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) represents a challenge to virus dispersion and the control of coronavirus disease 2019 (COVID-19) pandemic. The reason why some people have prolonged infection and how the virus persists for so long are still not fully understood. Recent studies suggested that the accumulation of intra-host single nucleotide variants (iSNVs) over the course of the infection might play an important role in persistence as well as emergence of mutations of concern. For this reason, we aimed to investigate the intra-host evolution of SARS-CoV-2 during prolonged infection. Thirty-three patients who remained reverse transcription polymerase chain reaction (RT-PCR) positive in the nasopharynx for on average 18 days from the symptoms onset were included in this study. Whole-genome sequences were obtained for each patient at two different time points. Phylogenetic, populational, and computational analyses of viral sequences were consistent with prolonged infection without evidence of coinfection in our cohort. We observed an elevated within-host genomic diversity at the second time point samples positively correlated with cycle threshold (Ct) values (lower viral load). Direct transmission was also confirmed in a small cluster of healthcare professionals that shared the same workplace by the presence of common iSNVs. A differential accumulation of missense variants between the time points was detected targeting crucial structural and non-structural proteins such as Spike and helicase. Interestingly, longitudinal acquisition of iSNVs in Spike protein coincided in many cases with SARS-CoV-2 reactive and predicted T cell epitopes. We observed a distinguishing pattern of mutations over the course of the infection mainly driven by increasing A→U and decreasing G→A signatures. G→A mutations may be associated with RNA-editing enzyme activities; therefore, the mutational profiles observed in our analysis were suggestive of innate immune mechanisms of the host cell defense. Therefore, we unveiled a dynamic and complex landscape of host and pathogen interaction during prolonged infection of SARS-CoV-2, suggesting that the host’s innate immunity shapes the increase of intra-host diversity. Our findings may also shed light on possible mechanisms underlying the emergence and spread of new variants resistant to the host immune response as recently observed in COVID-19 pandemic.

An integrated DNA and RNA variant detector identifies a highly conserved three base exon in the MAP4K5 kinase locus

RNA variants that emerge from editing and alternative splicing form important regulatory stages in protein signalling. In this report, we apply an integrated DNA and RNA variant detection workbench to define the range of RNA variants that deviate from the reference genome in a human melanoma cell model. The RNA variants can be grouped into (i) classic ADAR-like or APOBEC-like RNA editing events and (ii) multiple-nucleotide variants (MNVs) including three and six base pair in-frame non-canonical unmapped exons. We focus on validating representative genes of these classes. First, clustered non-synonymous RNA edits (A-I) in the CDK13 gene were validated by Sanger sequencing to confirm the integrity of the RNA variant detection workbench. Second, a highly conserved RNA variant in the MAP4K5 gene was detected that results most likely from the splicing of a non-canonical three-base exon. The two RNA variants produced from the MAP4K5 locus deviate from the genomic reference sequence and produce V569E or V569del isoform variants. Low doses of splicing inhibitors demonstrated that the MAP4K5-V569E variant emerges from an SF3B1-dependent splicing event. Mass spectrometry of the recombinant SBP-tagged MAP4K5^V569E and MAP4K5^V569del proteins pull-downs in transfected cell systems was used to identify the protein-protein interactions of these two MAP4K5 isoforms and propose possible functions. Together these data highlight the utility of this integrated DNA and RNA variant detection platform to detect RNA variants in cancer cells and support future analysis of RNA variant detection in cancer tissue.

Transcriptome-Wide Identification of G-to-A RNA Editing in Chronic Social Defeat Stress Mouse Models

Emerging evidence suggests that RNA editing is associated with stress, neurological diseases, and psychiatric disorders. However, the role of G-to-A RNA editing in chronic social defeat stress (CSDS) remains unclear. We herein identified G-to-A RNA editing and its changes in the ventral tegmental area (VTA), a key region of the brain reward system, in CSDS mouse models under emotional stress (ES) and physiological stress (PS) conditions. Our results revealed 3812 high-confidence G-to-A editing events. Among them, 56 events were significantly downregulated while 23 significantly upregulated in CSDS compared to controls. Moreover, divergent editing patterns were observed between CSDS mice under ES and PS conditions, with 42 and 21 events significantly upregulated in PS and ES, respectively. Interestingly, differential RNA editing was enriched in genes with multiple editing events. Genes differentially edited in CSDS included those genetically associated with mental or neurodevelopmental disorders, especially mood disorders, such as FAT atypical cadherin 1 and solute carrier family 6 member 1. Notably, changes of G-to-A RNA editing were also implicated in ionotropic glutamate receptors, a group of well-known targets of adenosine-to-inosine RNA editing. Such results demonstrate dynamic G-to-A RNA editing changes in the brain of CSDS mouse models, underlining its role as a potential molecular mechanism of CSDS and stress-related diseases.

APOBEC3G rescues cells from the deleterious effects of DNA damage

Human apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like 3G (hA3G), a member of the APOBEC family, was described as an anti-HIV-1 restriction factor, deaminating reverse transcripts of the HIV-1 genome. Several types of cancer cells that express high levels of A3G, such as diffuse large B-cell lymphoma cells and glioblastomas, show enhanced cell survival after ionizing radiation and chemotherapy treatments. Previously, we showed that hA3G promotes (DNA) double-strand breaks repair in cultured cells and rescues transgenic mice from a lethal dose of ionizing radiation. Here, we show that A3G rescues cells from the detrimental effects of DNA damage induced by ultraviolet irradiation and by combined bromodeoxyuridine and ultraviolet treatments. The combined treatments stimulate the synthesis of cellular proteins, which are exclusively associated with A3G expression. These proteins participate mainly in nucleotide excision repair and homologous recombination DNA repair pathways. Our results implicate A3G inhibition as a potential strategy for increasing tumor cell sensitivity to genotoxic treatments.

Predicting clinical outcomes using cancer progression associated signatures

Somatic mutation signatures are an informative facet of cancer aetiology, however they are rarely useful for predicting patient outcome. The aim of this study is to evaluate the utility of a panel of 142 mutation-signature–associated metrics (P142) for predicting cancer progression in patients from a ‘TCGA PanCancer Atlas’ cohort. The P142 metrics are comprised of AID/APOBEC and ADAR deaminase associated SNVs analyzed for codon context, strand bias, and transitions/transversions. TCGA tumor-normal mutation data was obtained for 10,437 patients, representing 31 of the most prevalent forms of cancer. Stratified random sampling was used to split patients into training, tuning and validation cohorts for each cancer type. Cancer specific machine learning (XGBoost) models were built using the output from the P142 panel to predict patient Progression Free Survival (PFS) status as either “High PFS” or “Low PFS”. Predictive performance of each model was evaluated using the validation cohort. Models accurately predicted PFS status for several cancer types, including adrenocortical carcinoma, glioma, mesothelioma, and sarcoma. In conclusion, the P142 panel of metrics successfully predicted cancer progression status in patients with some, but not all cancer types analyzed. These results pave the way for future studies on cancer progression associated signatures.

The interesting relationship between APOBEC3 deoxycytidine deaminases and cancer: a long road ahead

Cancer is considered a group of diseases characterized by uncontrolled growth and spread of abnormal cells and is propelled by somatic mutations. Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) family of enzymes are endogenous sources of somatic mutations found in multiple human cancers. While these enzymes normally act as an intrinsic immune defence against viruses, they can also catalyse 'off-target' cytidine deamination in genomic single-stranded DNA intermediates. The deamination of cytosine forms uracil, which is promutagenic in DNA. Key factors to trigger the APOBEC 'off-target' activity are overexpression in a non-normal cell type, nuclear localization and replication stress. The resulting uracil-induced mutations contribute to genomic variation, which may result in neutral, beneficial or harmful consequences for the cancer. This review summarizes the functional and biochemical basis of the APOBEC3 enzyme activity and highlights their relationship with the most well-studied cancers in this particular context such as breast, lung, bladder, and human papillomavirus-associated cancers. We focus on APOBEC3A, APOBEC3B and APOBEC3H haplotype I because they are the leading candidates as sources of somatic mutations in these and other cancers. Also, we discuss the prognostic value of the APOBEC3 expression in drug resistance and response to therapies.

Similarity between mutation spectra in hypermutated genomes of rubella virus and in SARS-CoV-2 genomes accumulated during the COVID-19 pandemic

Genomes of tens of thousands of SARS-CoV2 isolates have been sequenced across the world and the total number of changes (predominantly single base substitutions) in these isolates exceeds ten thousand. We compared the mutational spectrum in the new SARS-CoV-2 mutation dataset with the previously published mutation spectrum in hypermutated genomes of rubella-another positive single stranded (ss) RNA virus. Each of the rubella virus isolates arose by accumulation of hundreds of mutations during propagation in a single subject, while SARS-CoV-2 mutation spectrum represents a collection events in multiple virus isolates from individuals across the world. We found a clear similarity between the spectra of single base substitutions in rubella and in SARS-CoV-2, with C to U as well as A to G and U to C being the most prominent in plus strand genomic RNA of each virus. Of those, U to C changes universally showed preference for loops versus stems in predicted RNA secondary structure. Similarly, to what was previously reported for rubella virus, C to U changes showed enrichment in the uCn motif, which suggested a subclass of APOBEC cytidine deaminase being a source of these substitutions. We also found enrichment of several other trinucleotide-centered mutation motifs only in SARS-CoV-2-likely indicative of a mutation process characteristic to this virus. Altogether, the results of this analysis suggest that the mutation mechanisms that lead to hypermutation of the rubella vaccine virus in a rare pathological condition may also operate in the background of the SARS-CoV-2 viruses currently propagating in the human population.

Similarity between mutation spectra in hypermutated genomes of rubella virus and in SARS-CoV-2 genomes accumulated during the COVID-19 pandemic

Genomes of tens of thousands of SARS-CoV2 isolates have been sequenced across the world and the total number of changes (predominantly single base substitutions) in these isolates exceeds ten thousand. We compared the mutational spectrum in the new SARS-CoV-2 mutation dataset with the previously published mutation spectrum in hypermutated genomes of rubella - another positive single stranded (ss) RNA virus. Each of the rubella isolates arose by accumulation of hundreds of mutations during propagation in a single subject, while SARS-CoV-2 mutation spectrum represents a collection events in multiple virus isolates from individuals across the world. We found a clear similarity between the spectra of single base substitutions in rubella and in SARS-CoV-2, with C to U as well as A to G and U to C being the most prominent in plus strand genomic RNA of each virus. Of those, U to C changes universally showed preference for loops versus stems in predicted RNA secondary structure. Similarly, to what was previously reported for rubella, C to U changes showed enrichment in the uCn motif, which suggested a subclass of APOBEC cytidine deaminase being a source of these substitutions. We also found enrichment of several other trinucleotide-centered mutation motifs only in SARS-CoV-2 - likely indicative of a mutation process characteristic to this virus. Altogether, the results of this analysis suggest that the mutation mechanisms that lead to hypermutation of the rubella vaccine virus in a rare pathological condition may also operate in the background of the SARS-CoV-2 viruses currently propagating in the human population.

A comprehensive overview of bull sperm-borne small non-coding RNAs and their diversity across breeds

Background

Mature sperm carry thousands of RNAs, including mRNAs, lncRNAs, tRNAs, rRNAs and sncRNAs, though their functional significance is still a matter of debate. Growing evidence suggests that sperm RNAs, especially sncRNAs, are selectively retained during spermiogenesis or specifically transferred during epididymis maturation, and are thus delivered to the oocyte at fertilization, providing resources for embryo development. However , a deep characterization of the sncRNA content of bull sperm and its expression profile across breeds is currently lacking. To fill this gap, we optimized a guanidinium–Trizol total RNA extraction protocol to prepare high-quality RNA from frozen bull sperm collected from 40 representative bulls from six breeds. Deep sequencing was performed (40 M single 50-bp reads per sample) to establish a comprehensive repertoire of cattle sperm sncRNA.

Results

Our study showed that it comprises mostly piRNAs (26%), rRNA fragments (25%), miRNAs (20%) and tRNA fragments (tsRNA, 14%). We identified 5p-halves as the predominant tsRNA subgroup in bull sperm, originating mostly from Gly and Glu isoacceptors. Our study also increased by ~ 50% the sperm repertoire of known miRNAs and identified 2022 predicted miRNAs. About 20% of sperm miRNAs were located within genomic clusters, expanding the list of known polycistronic pri-miRNA clusters and defining several networks of co-expressed miRNAs. Strikingly, our study highlighted the great diversity of isomiRs, resulting mainly from deletions and non-templated additions (A and U) at the 3p end. Substitutions within miRNA sequence accounted for 40% of isomiRs, with G>A, U>C and C>U substitutions being the most frequent variations. In addition, many sncRNAs were found to be differentially expressed across breeds.

Conclusions

Our study provides a comprehensive overview of cattle sperm sncRNA, and these findings will pave the way for future work on the role of sncRNAs in embryo development and their relevance as biomarkers of semen fertility.

Single-nucleotide variants in human RNA: RNA editing and beyond

Through analysis of paired high-throughput DNA-Seq and RNA-Seq data, researchers quickly recognized that RNA-Seq can be used for more than just gene expression quantification. The alternative applications of RNA-Seq data are abundant, and we are particularly interested in its usefulness for detecting single-nucleotide variants, which arise from RNA editing, genomic variants and other RNA modifications. A stunning discovery made from RNA-Seq analyses is the unexpectedly high prevalence of RNA-editing events, many of which cannot be explained by known RNA-editing mechanisms. Over the past 6-7 years, substantial efforts have been made to maximize the potential of RNA-Seq data. In this review we describe the controversial history of mining RNA-editing events from RNA-Seq data and the corresponding development of methodologies to identify, predict, assess the quality of and catalog RNA-editing events as well as genomic variants.

Infectious vaccine-derived rubella viruses emerge, persist, and evolve in cutaneous granulomas of children with primary immunodeficiencies

Rubella viruses (RV) have been found in an association with granulomas in children with primary immune deficiencies (PID). Here, we report the recovery and characterization of infectious immunodeficiency-related vaccine-derived rubella viruses (iVDRV) from diagnostic skin biopsies of four patients. Sequence evolution within PID hosts was studied by comparison of the complete genomic sequences of the iVDRVs with the genome of the vaccine virus RA27/3. The degree of divergence of each iVDRV correlated with the duration of persistence indicating continuous intrahost evolution. The evolution rates for synonymous and nonsynonymous substitutions were estimated to be 5.7 x 10-3 subs/site/year and 8.9 x 10-4 subs/site/year, respectively. Mutational spectra and signatures indicated a major role for APOBEC cytidine deaminases and a secondary role for ADAR adenosine deaminases in generating diversity of iVDRVs. The distributions of mutations across the genes and 3D hotspots for amino acid substitutions in the E1 glycoprotein identified regions that may be under positive selective pressure. Quasispecies diversity was higher in granulomas than in recovered infectious iVDRVs. Growth properties of iVDRVs were assessed in WI-38 fibroblast cultures. None of the iVDRV isolates showed complete reversion to wild type phenotype but the replicative and persistence characteristics of iVDRVs were different from those of the RA27/3 vaccine strain, making predictions of iVDRV transmissibility and teratogenicity difficult. However, detection of iVDRV RNA in nasopharyngeal specimen and poor neutralization of some iVDRV strains by sera from vaccinated persons suggests possible public health risks associated with iVDRV carriers. Detection of IgM antibody to RV in sera of two out of three patients may be a marker of virus persistence, potentially useful for identifying patients with iVDRV before development of lesions. Studies of the evolutionary dynamics of iVDRV during persistence will contribute to development of infection control strategies and antiviral therapies.

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.

Widespread cis-regulation of RNA editing in a large mammal

Post-transcriptional RNA editing may regulate transcript expression and diversity in cells, with potential impacts on various aspects of physiology and environmental adaptation. A small number of recent genome-wide studies in Drosophila, mouse, and human have shown that RNA editing can be genetically modulated, highlighting loci that quantitatively impact editing of transcripts. The potential gene expression and physiological consequences of these RNA-editing quantitative trait loci (edQTL), however, are almost entirely unknown. Here, we present analyses of RNA editing in a large domestic mammal (Bos taurus), where we use whole-genome and high-depth RNA sequencing to discover, characterize, and conduct genetic mapping studies of novel transcript edits. Using a discovery population of nine deeply sequenced cows, we identify 2413 edit sites in the mammary transcriptome, the majority of which are adenosine to inosine edits (98.6%). Most sites are predicted to reside in double-stranded secondary structures (85.1%), and quantification of the rates of editing in an additional 355 cows reveals editing is negatively correlated with gene expression in the majority of cases. Genetic analyses of RNA editing and gene expression highlight 152 cis-regulated edQTL, of which 15 appear to cosegregate with expression QTL effects. Trait association analyses in a separate population of 9989 lactating cows also shows 12 of the cis-edQTL coincide with at least one cosegregating lactation QTL. Together, these results enhance our understanding of RNA-editing dynamics in mammals, and suggest mechanistic links by which loci may impact phenotype through RNA editing mediated processes.

The Role of RNA Editing in Cancer Development and Metabolic Disorders

Numerous human diseases arise from alterations of genetic information, most notably DNA mutations. Thought to be merely the intermediate between DNA and protein, changes in RNA sequence were an afterthought until the discovery of RNA editing 30 years ago. RNA editing alters RNA sequence without altering the sequence or integrity of genomic DNA. The most common RNA editing events are A-to-I changes mediated by adenosine deaminase acting on RNA (ADAR), and C-to-U editing mediated by apolipoprotein B mRNA editing enzyme, catalytic polypeptide 1 (APOBEC1). Both A-to-I and C-to-U editing were first identified in the context of embryonic development and physiological homeostasis. The role of RNA editing in human disease has only recently started to be understood. In this review, the impact of RNA editing on the development of cancer and metabolic disorders will be examined. Distinctive functions of each RNA editase that regulate either A-to-I or C-to-U editing will be highlighted in addition to pointing out important regulatory mechanisms governing these processes. The potential of developing novel therapeutic approaches through intervention of RNA editing will be explored. As the role of RNA editing in human disease is elucidated, the clinical utility of RNA editing targeted therapies will be needed. This review aims to serve as a bridge of information between past findings and future directions of RNA editing in the context of cancer and metabolic disease.

[Donor chimerism and minimal residual disease monitoring in leukemia patients after allo-HSCT]

In present research the comparative analysis of donor chimerism (DC) using different tests was performed to improve the diagnostic tool in patients with malignant hematological disorders after allo-HSCT. The RBC antigen typing, identification of ABO blood type and quantitative analysis of InDel-, STR-, Y-polymorphisms were carried out for detection of DC. In addition, the expression of well-known oncogenes and CD-markers for monitoring MRD was evaluated to predict relapse and clinical outcome. According to our research, the analysis of InDel polymorphism using AlleleSEQR-PCR is more sensitive test for estimation of DC as compared with other assays. Moreover, the sensitivity of AlleleSEQR-PCR may be increased after isolation of the CD34 cell population in bone marrow. Nevertheless, observation of high levels in DC (³95%) in some leukemia patients (ALL, Ph+, bcr-abl/p190+) during first 6 months after HSCT cannot exclude the possibility of relapse. Thus, the combined monitoring of both DC (InDel) and MRD (oncogenes, WT1 and CD-markers) is a more advisable and useful test in managing hematologic malignancies and predicting relapse risk after allo-HSCT.

Roles of APOBEC3A and APOBEC3B in Human Papillomavirus Infection and Disease Progression

The apolipoprotein B messenger RNA-editing, enzyme-catalytic, polypeptide-like 3 (APOBEC3) family of cytidine deaminases plays an important role in the innate immune response to viral infections by editing viral genomes. However, the cytidine deaminase activity of APOBEC3 enzymes also induces somatic mutations in host genomes, which may drive cancer progression. Recent studies of human papillomavirus (HPV) infection and disease outcome highlight this duality. HPV infection is potently inhibited by one family member, APOBEC3A. Expression of APOBEC3A and APOBEC3B is highly elevated by the HPV oncoproteins E6 and E7 during persistent virus infection and disease progression. Furthermore, there is a high prevalence of APOBEC3A and APOBEC3B mutation signatures in HPV-associated cancers. These findings suggest that induction of an APOBEC3-mediated antiviral response during HPV infection may inadvertently contribute to cancer mutagenesis and virus evolution. Here, we discuss current understanding of APOBEC3A and APOBEC3B biology in HPV restriction, evolution, and associated cancer mutagenesis.

APOBEC3B, a molecular driver of mutagenesis in human cancers

Human cancers results in large part from the accumulation of multiple mutations. The progression of premalignant cells is an evolutionary process in which mutations provide the fundamental driving force for genetic diversity. The increased mutation rate in premalignant cells allows selection for increased proliferation and survival and ultimately leads to invasion, metastasis, recurrence, and therapeutic resistance. Therefore, it is important to understand the molecular determinants of the mutational processes. Recent genome-wide sequencing data showed that apolipoprotein B mRNA editing catalytic polypeptide-like 3B (APOBEC3B) is a key molecular driver inducing mutations in multiple human cancers. APOBEC3B, a DNA cytosine deaminase, is overexpressed in a wide spectrum of human cancers. Its overexpression and aberrant activation lead to unexpected clusters of mutations in the majority of cancers. This phenomenon of clustered mutations, termed kataegis (from the Greek word for showers), forms unique mutation signatures. In this review, we will discuss the biological function of APOBEC3B, its tumorigenic role in promoting mutational processes in cancer development and the clinical potential to develop novel therapeutics by targeting APOBEC3B.

In vitro recapitulation of the site-specific editing (to wild-type) of mutant IDS mRNA transcripts, and the characterization of IDS protein translated from the edited mRNAs

The transfer of genomic information into the primary RNA sequence can be altered by RNA editing. We have previously shown that genomic variants can be RNA-edited to wild-type. The presence of distinct "edited" iduronate 2-sulfatase (IDS) mRNA transcripts ex vivo evidenced the correction of a nonsense and frameshift variant, respectively, in three unrelated Hunter syndrome patients. This phenomenon was confirmed in various patient samples by a variety of techniques, and was quantified by single-nucleotide primer extension. Western blotting also confirmed the presence of IDS protein similar in size to the wild-type. Since preliminary experimental evidence suggested that the "corrected" IDS proteins produced by the patients were similar in molecular weight and net charge to their wild-type counterparts, an in vitro system employing different cell types was established to recapitulate the site-specific editing of IDS RNA (uridine to cytidine conversion and uridine deletion), and to confirm the findings previously observed ex vivo in the three patients. In addition, confocal microscopy and flow cytometry analyses demonstrated the expression and lysosomal localization in HEK293 cells of GFP-labeled proteins translated from edited IDS mRNAs. Confocal high-content analysis of the two patients' cells expressing wild-type or mutated IDS confirmed lysosomal localization and showed no accumulation in the Golgi or early endosomes.

RNA Editing-Systemic Relevance and Clue to Disease Mechanisms?

Recent advances in sequencing technologies led to the identification of a plethora of different genes and several hundreds of amino acid recoding edited positions. Changes in editing rates of some of these positions were associated with diseases such as atherosclerosis, myopathy, epilepsy, major depression disorder, schizophrenia and other mental disorders as well as cancer and brain tumors. This review article summarizes our current knowledge on that front and presents glycine receptor C-to-U RNA editing as a first example of disease-associated increased RNA editing that includes assessment of disease mechanisms of the corresponding gene product in an animal model.

Functions and Malfunctions of Mammalian DNA-Cytosine Deaminases

The AID/APOBEC family enzymes convert cytosines in single-stranded DNA to uracils, causing base substitutions and strand breaks. They are induced by cytokines produced during the body's inflammatory response to infections, and they help combat infections through diverse mechanisms. AID is essential for the maturation of antibodies and causes mutations and deletions in antibody genes through somatic hypermutation (SHM) and class-switch recombination (CSR) processes. One member of the APOBEC family, APOBEC1, edits mRNA for a protein involved in lipid transport. Members of the APOBEC3 subfamily in humans (APOBEC3A, APOBEC3B, APOBEC3C, APOBEC3D, APOBEC3F, APOBEC3G, and APOBEC3H) inhibit infections of viruses such as HIV-1, HBV, and HCV, and retrotransposition of endogenous retroelements through mutagenic and nonmutagenic mechanisms. There is emerging consensus that these enzymes can cause mutations in the cellular genome at replication forks or within transcription bubbles depending on the physiological state of the cell and the phase of the cell cycle during which they are expressed. We describe here the state of knowledge about the structures of these enzymes, regulation of their expression, and both the advantageous and deleterious consequences of their expression, including carcinogenesis. We highlight similarities among them and present a holistic view of their regulation and function.

RNA-DNA sequence differences in Saccharomyces cerevisiae

Alterations of RNA sequences and structures, such as those from editing and alternative splicing, result in two or more RNA transcripts from a DNA template. It was thought that in yeast, RNA editing only occurs in tRNAs. Here, we found that Saccharomyces cerevisiae have all 12 types of RNA–DNA sequence differences (RDDs) in the mRNA. We showed these sequence differences are propagated to proteins, as we identified peptides encoded by the RNA sequences in addition to those by the DNA sequences at RDD sites. RDDs are significantly enriched at regions with R-loops. A screen of yeast mutants showed that RDD formation is affected by mutations in genes regulating R-loops. Loss-of-function mutations in ribonuclease H, senataxin, and topoisomerase I that resolve RNA–DNA hybrids lead to increases in RDD frequency. Our results demonstrate that RDD is a conserved process that diversifies transcriptomes and proteomes and provide a mechanistic link between R-loops and RDDs.

A chimeric human APOBEC3A protein with a three amino acid insertion confers differential HIV-1 and adeno-associated virus restriction

Old World monkey (OWM) and hominid APOBEC3Aproteins exhibit differential restriction activities against lentiviruses and DNA viruses. Human APOBEC3A(hA3A)has weak restriction activity against HIV-1Δvifbut is efficiently restricted by an artificially generated chimeric from mandrills (mndA3A/G). We show that a chimeric hA3Acontaining the "WVS" insertion (hA3A[(27)WVS(29)]) conferred potent HIV-1restriction activity. Analysis of each amino acid of the "WVS" motif show that the length and not necessarily the charge or hydrophobicity of the amino acids accounted for restriction activity. Our results suggest that hA3A[(27)WVS(29)]restricts HIV-1at the level of reverse transcription in target cells. Finally, our results suggest that insertion of "WVS" into hA3Amodestly reduces restriction of adeno-associated virus 2(AAV-2)while insertion of the AC Loop1region of the mndA3A/G into hA3A abolished AAV-2 restriction, strengthening the role of this molecular interface in the functional evolution of primate A3A.

Translating the epitranscriptome

RNA modifications are indispensable for the translation machinery to provide accurate and efficient protein synthesis. Whereas the importance of transfer RNA (tRNA) and ribosomal RNA (rRNA) modifications has been well described and is unquestioned for decades, the significance of internal messenger RNA (mRNA) modifications has only recently been revealed. Novel experimental methods have enabled the identification of thousands of modified sites within the untranslated and translated regions of mRNAs. Thus far, N⁶‐methyladenosine (m⁶A), pseudouridine (Ψ), 5‐methylcytosine (m⁵C) and N¹‐methyladenosine (m¹A) were identified in eukaryal, and to some extent in prokaryal mRNAs. Several of the functions of these mRNA modifications have previously been reported, but many aspects remain elusive. Modifications can be important factors for the direct regulation of protein synthesis. The potential diversification of genomic information and regulation of RNA expression through editing and modifying mRNAs is versatile and many questions need to be addressed to completely elucidate the role of mRNA modifications. Herein, we summarize and highlight some recent findings on various co‐ and post‐transcriptional modifications, describing the impact of these processes on gene expression, with emphasis on protein synthesis. WIREs RNA 2017, 8:e1375. doi: 10.1002/wrna.1375

For further resources related to this article, please visit the WIREs website.

Functional requirements of AID's higher order structures and their interaction with RNA-binding proteins

Activation-induced cytidine deaminase (AID) is essential for the somatic hypermutation (SHM) and class-switch recombination (CSR) of Ig genes. Although both the N and C termini of AID have unique functions in DNA cleavage and recombination, respectively, during SHM and CSR, their molecular mechanisms are poorly understood. Using a bimolecular fluorescence complementation (BiFC) assay combined with glycerol gradient fractionation, we revealed that the AID C terminus is required for a stable dimer formation. Furthermore, AID monomers and dimers form complexes with distinct heterogeneous nuclear ribonucleoproteins (hnRNPs). AID monomers associate with DNA cleavage cofactor hnRNP K whereas AID dimers associate with recombination cofactors hnRNP L, hnRNP U, and Serpine mRNA-binding protein 1. All of these AID/ribonucleoprotein associations are RNA-dependent. We propose that AID's structure-specific cofactor complex formations differentially contribute to its DNA-cleavage and recombination functions.

Role of the single deaminase domain APOBEC3A in virus restriction, retrotransposition, DNA damage and cancer

The apolipoprotein mRNA editing enzyme catalytic polypeptide-like 3 (APOBEC3; A3) proteins are a family of seven cytidine deaminases (A3A, A3B, A3C, A3D, A3F, A3G and A3H) that restrict certain viral infections. These innate defence factors are best known for their ability to restrict the replication of human immunodeficiency virus type 1 (HIV-1) lacking a functional Vif protein (HIV-1Δvif) through the deamination of cytidine residues to uridines during reverse transcription, ultimately leading to lethal G → A changes in the viral genome. The best studied of the A3 proteins has been APOBEC3G because of its potent activity against HIV-1Δvif. However, one member of this family, A3A, has biological properties that make it unique among the A3 proteins. In this review, we will focus on the structural and phylogenetic features of the human and non-human primate A3A proteins, their role in the restriction of retroviruses and other viruses, and current findings on other biological properties affected by this protein.