HIV-1 RNAs whose transcription initiates from the third deoxyguanosine of GGG tract in the 5' long terminal repeat serve as a dominant genome for efficient provirus DNA formation.

Unspliced HIV-1 RNAs function as messenger RNAs for Gag or Gag-Pol polyproteins and progeny genomes packaged into virus particles. Recently, it has been reported that fate of the RNAs might be primarily determined, depending on transcriptional initiation sites among three consecutive deoxyguanosine residues (GGG tract) downstream of TATA-box in the 5' long terminal repeat (LTR). Although HIV-1 RNA transcription starts mostly from the first deoxyguanosine of the GGG tract and often from the second or third deoxyguanosine, RNAs beginning with one guanosine (G1-form RNAs), whose transcription initiates from the third deoxyguanosine, were predominant in HIV-1 particles. Despite selective packaging of G1-form RNAs into virus particles, its biological impact during viral replication remains to be determined. In this study, we revealed that G1-form RNAs are primarily selected as a template for provirus DNA rather than other RNAs. In competitions between HIV-1 and lentiviral vector transcripts in virus-producing cells, approximately 80% of infectious particles were found to generate provirus using HIV-1 transcripts, while lentiviral vector transcripts were conversely selected when we used HIV-1 mutants in which the third deoxyguanosine in the GGG tract was replaced with deoxythymidine or deoxycytidine (GGT or GGC mutants, respectively). In the other analyses of proviral sequences after infection with an HIV-1 mutant in which the GGG tract in 3' LTR was replaced with TTT, most proviral sequences of the GGG-tract region in 5' LTR were found to be TTG, which is reasonably generated using the G1-form transcripts. Our results indicate that the G1-form RNAs serve as a dominant genome to establish provirus DNA.IMPORTANCESince the promoter for transcribing HIV-1 RNA is unique, all viral elements including genomic RNA and viral proteins have to be generated by the unique transcripts through ingenious mechanisms including RNA splicing and frameshifting during protein translation. Previous studies suggested a new mechanism for diversification of HIV-1 RNA functions by heterogeneous transcriptional initiation site usage; HIV-1 RNAs whose transcription initiates from a certain nucleotide were predominant in virus particles. In this study, we established two methods to analyze heterogenous transcriptional initiation site usage by HIV-1 during viral infection and showed that RNAs beginning with one guanosine (G1-form RNAs), whose transcription initiates from the third deoxyguanosine of the GGG tract in 5' LTR, were primarily selected as viral genome in infectious particles and thus are used as a template to generate provirus for continuous replication. This study provides insights into the mechanism for diversification of unspliced RNA functions and requisites of lentivirus infectivity.

m7GHub V2.0: an updated database for decoding the N7-methylguanosine (m7G) epitranscriptome.

With recent progress in mapping N7-methylguanosine (m7G) RNA methylation sites, tens of thousands of experimentally validated m7G sites have been discovered in various species, shedding light on the significant role of m7G modification in regulating numerous biological processes including disease pathogenesis. An integrated resource that enables the sharing, annotation and customized analysis of m7G data will greatly facilitate m7G studies under various physiological contexts. We previously developed the m7GHub database to host mRNA m7G sites identified in the human transcriptome. Here, we present m7GHub v.2.0, an updated resource for a comprehensive collection of m7G modifications in various types of RNA across multiple species: an m7GDB database containing 430 898 putative m7G sites identified in 23 species, collected from both widely applied next-generation sequencing (NGS) and the emerging Oxford Nanopore direct RNA sequencing (ONT) techniques; an m7GDiseaseDB hosting 156 206 m7G-associated variants (involving addition or removal of an m7G site), including 3238 disease-relevant m7G-SNPs that may function through epitranscriptome disturbance; and two enhanced analysis modules to perform interactive analyses on the collections of m7G sites (m7GFinder) and functional variants (m7GSNPer). We expect that m7Ghub v.2.0 should serve as a valuable centralized resource for studying m7G modification. It is freely accessible at: www.rnamd.org/m7GHub2.

The potential role of N7-methylguanosine (m7G) in cancer.

N7-methylguanosine (m7G), one of the most prevalent RNA modifications, has recently attracted significant attention. The m7G modification actively participates in biological and pathological functions by affecting the metabolism of various RNA molecules, including messenger RNA, ribosomal RNA, microRNA, and transfer RNA. Increasing evidence indicates a critical role for m7G in human disease development, especially cancer, and aberrant m7G levels are closely associated with tumorigenesis and progression via regulation of the expression of multiple oncogenes and tumor suppressor genes. Currently, the underlying molecular mechanisms of m7G modification in cancer are not comprehensively understood. Here, we review the current knowledge regarding the potential function of m7G modifications in cancer and discuss future m7G-related diagnostic and therapeutic strategies.

Oxidative damage and DNA repair in desiccated recalcitrant embryonic axes of Acer pseudoplatanus L.

BACKGROUND: Most plants encounter water stress at one or more different stages of their life cycle. The maintenance of genetic stability is the integral component of desiccation tolerance that defines the storage ability and long-term survival of seeds. Embryonic axes of desiccation-sensitive recalcitrant seeds of Acer pseudoplatnus L. were used to investigate the genotoxic effect of desiccation. Alkaline single-cell gel electrophoresis (comet assay) methodology was optimized and used to provide unique insights into the onset and repair of DNA strand breaks and 8-oxo-7,8-dihydroguanine (8-oxoG) formation during progressive steps of desiccation and rehydration. RESULTS: The loss of DNA integrity and impairment of damage repair were significant predictors of the viability of embryonic axes. In contrast to the comet assay, automated electrophoresis failed to detect changes in DNA integrity resulting from desiccation. Notably, no significant correlation was observed between hydroxyl radical (Ù OH) production and 8-oxoG formation, although the former is regarded to play a major role in guanine oxidation. CONCLUSIONS: The high-throughput comet assay represents a sensitive tool for monitoring discrete changes in DNA integrity and assessing the viability status in plant germplasm processed for long-term storage.

Identifying N7-methylguanosine sites by integrating multiple features.

Recent studies reported that N7-methylguanosine (m7G) plays a vital role in gene expression regulation. As a consequence, determining the distribution of m7G is a crucial step towards further understanding its biological functions. Although biological experimental approaches are capable of accurately locating m7G sites, they are labor-intensive, costly, and time-consuming. Therefore, it is necessary to develop more effective and robust computational methods to replace, or at least complement current experimental methods. In this study, we developed a novel sequence-based computational tool to identify RNA m7G sites. In this model, 22 kinds of dinucleotide physicochemical (PC) properties were employed to encode the RNA sequence. Three types of descriptors, including auto-covariance, cross-covariance, and discrete wavelet transform were adopted to extract effective features from the PC matrix. The least absolute shrinkage and selection operator (LASSO) algorithm was utilized to reduce the influence of irrelevant or redundant features. Finally, these selected features were fed into a support vector machine (SVM) for distinguishing m7G from non-m7G sites. The proposed method significantly outperforms existing predictors across all evaluation metrics. It indicates that the approach is effective in identifying RNA m7G sites.

RFhy-m2G: Identification of RNA N2-methylguanosine modification sites based on random forest and hybrid features.

N2-methylguanosine is a post-transcriptional modification of RNA that is found in eukaryotes and archaea. The biological function of m2G modification discovered so far is to control and stabilize the three-dimensional structure of tRNA and the dynamic barrier of reverse transcription. To discover additional biological functions of m2G, it is necessary to develop time-saving and labor-saving calculation tools to identify m2G. In this paper, based on hybrid features and a random forest, a novel predictor, RFhy-m2G, was developed to identify the m2G modification sites for three species. The hybrid feature used by the predictor is used to fuse the three features of ENAC, PseDNC, and NPPS. These three features include primary sequence derivation properties, physicochemical properties, and position-specific properties. Since there are redundant features in hybrid features, MRMD2.0 is used for optimal feature selection. Through feature analysis, it is found that the optimal hybrid features obtained still contain three kinds of properties, and the hybrid features can more accurately identify m2G modification sites and improve prediction performance. Based on five-fold cross-validation and independent testing to evaluate the prediction model, the accuracies obtained were 0.9982 and 0.9417, respectively. The robustness of the predictor is demonstrated by comparisons with other predictors.

m7G-DPP: Identifying N7-methylguanosine sites based on dinucleotide physicochemical properties of RNA.

N7-methylguanosine (m7G) modification is one of the most common post-transcriptional RNA modifications, which play vital role in the regulation of gene expression. Dysfunction of m7G may result to developmental defects and the appearance of some serious diseases. Thus, it is an urgent task to fast and accurate identifying m7G sites. In view of experimental approaches are costly and time-consuming, researchers focused their attention on computational models. Hence, in current study, we proposed a novel predictor called m7G-DPP to identify m7G sites. In the predictor, the RNA sequences were firstly encoded by physicochemical (PC) properties of dinucleotide. Then, sliding window approach was adopted to divide PC matrix into multiple matrixes, and Pearson's correlation coefficient (PCC), dynamic time warping (DTW), and distance correlation (DC) were employed to extract classification features at each window. Next, the least absolute shrinkage and selection operator (LASSO) algorithm was applied to select discriminative features. Finally, these selected features were fed into support vector machine to identify m7G sites. Experimental results showed that the proposed method is effective, which may play a complementary role in current m7G sites prediction studies. The MATLAB codes and dataset can be obtained from website at https://figshare.com/articles/online_resource/m7G-DPP/15000348.

BERT-m7G: A Transformer Architecture Based on BERT and Stacking Ensemble to Identify RNA N7-Methylguanosine Sites from Sequence Information.

As one of the most prevalent posttranscriptional modifications of RNA, N7-methylguanosine (m7G) plays an essential role in the regulation of gene expression. Accurate identification of m7G sites in the transcriptome is invaluable for better revealing their potential functional mechanisms. Although high-throughput experimental methods can locate m7G sites precisely, they are overpriced and time-consuming. Hence, it is imperative to design an efficient computational method that can accurately identify the m7G sites. In this study, we propose a novel method via incorporating BERT-based multilingual model in bioinformatics to represent the information of RNA sequences. Firstly, we treat RNA sequences as natural sentences and then employ bidirectional encoder representations from transformers (BERT) model to transform them into fixed-length numerical matrices. Secondly, a feature selection scheme based on the elastic net method is constructed to eliminate redundant features and retain important features. Finally, the selected feature subset is input into a stacking ensemble classifier to predict m7G sites, and the hyperparameters of the classifier are tuned with tree-structured Parzen estimator (TPE) approach. By 10-fold cross-validation, the performance of BERT-m7G is measured with an ACC of 95.48% and an MCC of 0.9100. The experimental results indicate that the proposed method significantly outperforms state-of-the-art prediction methods in the identification of m7G modifications.

Defects in 8-oxo-guanine repair pathway cause high frequency of C > A substitutions in neuroblastoma.

Neuroblastomas are childhood tumors with frequent fatal relapses after induction treatment, which is related to tumor evolution with additional genomic events. Our whole-genome sequencing data analysis revealed a high frequency of somatic cytosine > adenine (C > A) substitutions in primary neuroblastoma tumors, which was associated with poor survival. We showed that increased levels of C > A substitutions correlate with copy number loss (CNL) of OGG1 or MUTYH Both genes encode DNA glycosylases that recognize 8-oxo-guanine (8-oxoG) lesions as a first step of 8-oxoG repair. Tumor organoid models with CNL of OGG1 or MUTYH show increased 8-oxoG levels compared to wild-type cells. We used CRISPR-Cas9 genome editing to create knockout clones of MUTYH and OGG1 in neuroblastoma cells. Whole-genome sequencing of single-cell OGG1 and MUTYH knockout clones identified an increased accumulation of C > A substitutions. Mutational signature analysis of these OGG1 and MUTYH knockout clones revealed enrichment for C > A signatures 18 and 36, respectively. Clustering analysis showed that the knockout clones group together with tumors containing OGG1 or MUTYH CNL. In conclusion, we demonstrate that defects in 8-oxoG repair cause accumulation of C > A substitutions in neuroblastoma, which contributes to mutagenesis and tumor evolution.

N7-Methylguanosine tRNA modification enhances oncogenic mRNA translation and promotes intrahepatic cholangiocarcinoma progression.

Cancer cells selectively promote translation of specific oncogenic transcripts to facilitate cancer survival and progression, but the underlying mechanisms are poorly understood. Here, we find that N7-methylguanosine (m7G) tRNA modification and its methyltransferase complex components, METTL1 and WDR4, are significantly upregulated in intrahepatic cholangiocarcinoma (ICC) and associated with poor prognosis. We further reveal the critical role of METTL1/WDR4 in promoting ICC cell survival and progression using loss- and gain-of-function assays in vitro and in vivo. Mechanistically, m7G tRNA modification selectively regulates the translation of oncogenic transcripts, including cell-cycle and epidermal growth factor receptor (EGFR) pathway genes, in m7G-tRNA-decoded codon-frequency-dependent mechanisms. Moreover, using overexpression and knockout mouse models, we demonstrate the crucial oncogenic function of Mettl1-mediated m7G tRNA modification in promoting ICC tumorigenesis and progression in vivo. Our study uncovers the important physiological function and mechanism of METTL1-mediated m7G tRNA modification in the regulation of oncogenic mRNA translation and cancer progression.

METTL1-mediated m7G modification of Arg-TCT tRNA drives oncogenic transformation.

The emerging "epitranscriptomics" field is providing insights into the biological and pathological roles of different RNA modifications. The RNA methyltransferase METTL1 catalyzes N7-methylguanosine (m7G) modification of tRNAs. Here we find METTL1 is frequently amplified and overexpressed in cancers and is associated with poor patient survival. METTL1 depletion causes decreased abundance of m7G-modified tRNAs and altered cell cycle and inhibits oncogenicity. Conversely, METTL1 overexpression induces oncogenic cell transformation and cancer. Mechanistically, we find increased abundance of m7G-modified tRNAs, in particular Arg-TCT-4-1, and increased translation of mRNAs, including cell cycle regulators that are enriched in the corresponding AGA codon. Accordingly, Arg-TCT expression is elevated in many tumor types and is associated with patient survival, and strikingly, overexpression of this individual tRNA induces oncogenic transformation. Thus, METTL1-mediated tRNA modification drives oncogenic transformation through a remodeling of the mRNA "translatome" to increase expression of growth-promoting proteins and represents a promising anti-cancer target.

Preferential import of queuosine-modified tRNAs into Trypanosoma brucei mitochondrion is critical for organellar protein synthesis.

Transfer RNAs (tRNAs) are key players in protein synthesis. To be fully active, tRNAs undergo extensive post-transcriptional modifications, including queuosine (Q), a hypermodified 7-deaza-guanosine present in the anticodon of several tRNAs in bacteria and eukarya. Here, molecular and biochemical approaches revealed that in the protozoan parasite Trypanosoma brucei, Q-containing tRNAs have a preference for the U-ending codons for asparagine, aspartate, tyrosine and histidine, analogous to what has been described in other systems. However, since a lack of tRNA genes in T. brucei mitochondria makes it essential to import a complete set from the cytoplasm, we surprisingly found that Q-modified tRNAs are preferentially imported over their unmodified counterparts. In turn, their absence from mitochondria has a pronounced effect on organellar translation and affects function. Although Q modification in T. brucei is globally important for codon selection, it is more so for mitochondrial protein synthesis. These results provide a unique example of the combined regulatory effect of codon usage and wobble modifications on protein synthesis; all driven by tRNA intracellular transport dynamics.

Transcriptome-Wide Detection of Internal N7-Methylguanosine.

m7G-seq detects internal 7-methylguanosine (m7G) sites within mRNAs and noncoding RNAs by misincorporation signatures. A chemical-assisted sequencing approach selectively converts internal m7G sites into abasic sites, triggering misincorporation at these sites in the presence of a specific reverse transcriptase. The further enrichment of m7G-induced abasic sites by biotin pull-down reveals hundreds of internal m7G sites in human mRNA. The misincorporation ratio before pull-down enrichment can be used for estimating the methylation fraction of some highly methylated m7G sites.

AlkAniline-Seq: A Highly Sensitive and Specific Method for Simultaneous Mapping of 7-Methyl-guanosine (m7G) and 3-Methyl-cytosine (m3C) in RNAs by High-Throughput Sequencing.

Epitranscriptomics is an emerging field where the development of high-throughput analytical technologies is essential to profile the dynamics of RNA modifications under different conditions. Despite important advances during the last 10 years, the number of RNA modifications detectable by next-generation sequencing is restricted to a very limited subset. Here, we describe a highly efficient and fast method called AlkAniline-Seq to map simultaneously two different RNA modifications: 7-methyl-guanosine (m7G) and 3-methyl-cytosine (m3C) in RNA. Our protocol is based on three subsequent chemical/enzymatic steps allowing the enrichment of RNA fragments ending at position n + 1 to the modified nucleotide, without any prior RNA selection. Therefore, AlkAniline-Seq demonstrates an outstanding sensitivity and specificity for these two RNA modifications. We have validated AlkAniline-Seq using bacterial, yeast, and human total RNA, and here we present, as an example, a synthetic view of the complete profiling of these RNA modifications in S. cerevisiae tRNAs.

Detecting Internal N7-Methylguanosine mRNA Modifications by Differential Enzymatic Digestion Coupled with Mass Spectrometry Analysis.

The recent discovery of reversible chemical modifications on mRNA has opened a new era of post-transcriptional gene regulation in eukaryotes. Among these modifications identified in eukaryotic mRNA, N7-methylguanosine (m7G) is unique owing to its presence in the 5' cap structure. Recently, it has been reported that m7G also exists internally in mRNA. Here, we describe a protocol of combining differential enzymatic digestion with liquid chromatography-electrospray ionization-tandem mass spectrometry (LC-ESI-MS/MS) analysis to detect internal m7G modification in mRNA. This protocol can also be used to quantify the level of m7G at both the 5' cap and internal positions of mRNA.

Iterative feature representation algorithm to improve the predictive performance of N7-methylguanosine sites.

MOTIVATION: N7-methylguanosine (m7G) is an important epigenetic modification, playing an essential role in gene expression regulation. Therefore, accurate identification of m7G modifications will facilitate revealing and in-depth understanding their potential functional mechanisms. Although high-throughput experimental methods are capable of precisely locating m7G sites, they are still cost ineffective. Therefore, it's necessary to develop new methods to identify m7G sites. RESULTS: In this work, by using the iterative feature representation algorithm, we developed a machine learning based method, namely m7G-IFL, to identify m7G sites. To demonstrate its superiority, m7G-IFL was evaluated and compared with existing predictors. The results demonstrate that our predictor outperforms existing predictors in terms of accuracy for identifying m7G sites. By analyzing and comparing the features used in the predictors, we found that the positive and negative samples in our feature space were more separated than in existing feature space. This result demonstrates that our features extracted more discriminative information via the iterative feature learning process, and thus contributed to the predictive performance improvement.

Bacterial tRNA 2'-O-methylation is dynamically regulated under stress conditions and modulates innate immune response.

RNA modifications are a well-recognized way of gene expression regulation at the post-transcriptional level. Despite the importance of this level of regulation, current knowledge on modulation of tRNA modification status in response to stress conditions is far from being complete. While it is widely accepted that tRNA modifications are rather dynamic, such variations are mostly assessed in terms of total tRNA, with only a few instances where changes could be traced to single isoacceptor species. Using Escherichia coli as a model system, we explored stress-induced modulation of 2'-O-methylations in tRNAs by RiboMethSeq. This analysis and orthogonal analytical measurements by LC-MS show substantial, but not uniform, increase of the Gm18 level in selected tRNAs under mild bacteriostatic antibiotic stress, while other Nm modifications remain relatively constant. The absence of Gm18 modification in tRNAs leads to moderate alterations in E. coli mRNA transcriptome, but does not affect polysomal association of mRNAs. Interestingly, the subset of motility/chemiotaxis genes is significantly overexpressed in ΔTrmH mutant, this corroborates with increased swarming motility of the mutant strain. The stress-induced increase of tRNA Gm18 level, in turn, reduced immunostimulation properties of bacterial tRNAs, which is concordant with the previous observation that Gm18 is a suppressor of Toll-like receptor 7 (TLR7)-mediated interferon release. This documents an effect of stress induced modulation of tRNA modification that acts outside protein translation.

Experimental and theoretical rationalization for the base pairing abilities of inosine, guanosine, adenosine, and their corresponding 8-oxo-7,8-dihydropurine, and 8-bromopurine analogues within A-form duplexes of RNA.

Inosine is an important RNA modification, furthermore RNA oxidation has gained interest due, in part, to its potential role in the development/progression of disease as well as on its impact on RNA structure and function. In this report we established the base pairing abilities of purine nucleobases G, I, A, as well as their corresponding, 8-oxo-7,8-dihydropurine (common products of oxidation at the C8-position of purines), and 8-bromopurine (as probes to explore conformational changes), derivatives, namely 8-oxoG, 8-oxoI, 8-oxoA, 8-BrG, and 8-BrI. Dodecamers of RNA were obtained using standard phosphoramidite chemistry via solid-phase synthesis, and used as models to establish the impact that each of these nucleobases have on the thermal stability of duplexes, when base pairing to canonical and noncanonical nucleobases. Thermal stabilities were obtained from thermal denaturation transition (Tm ) measurements, via circular dichroism (CD). The results were then rationalized using models of base pairs between two monomers, via density functional theory (DFT), that allowed us to better understand potential contributions from H-bonding patterns arising from distinct conformations. Overall, some of the important results indicate that: (a) an anti-I:syn-A base pair provides thermal stability, due to the absence of the exocyclic amine; (b) 8-oxoG base pairs like U, and does not induce destabilization within the duplex when compared to the pyrimidine ring; (c) a U:G wobble-pair is only stabilized by G; and (d) 8-oxoA displays an inherited base pairing promiscuity in this sequence context. Gaining a better understanding of how this oxidatively generated lesions potentially base pair with other nucleobases will be useful to predict various biological outcomes, as well as in the design of biomaterials and/or nucleotide derivatives with biological potential.

Loss of the RNA trimethylguanosine cap is compatible with nuclear accumulation of spliceosomal snRNAs but not pre-mRNA splicing or snRNA processing during animal development.

The 2,2,7-trimethylguanosine (TMG) cap is one of the first identified modifications on eukaryotic RNAs. TMG, synthesized by the conserved Tgs1 enzyme, is abundantly present on snRNAs essential for pre-mRNA splicing. Results from ex vivo experiments in vertebrate cells suggested that TMG ensures nuclear localization of snRNAs. Functional studies of TMG using tgs1 mutations in unicellular organisms yield results inconsistent with TMG being indispensable for either nuclear import or splicing. Utilizing a hypomorphic tgs1 mutation in Drosophila, we show that TMG reduction impairs germline development by disrupting the processing, particularly of introns with smaller sizes and weaker splice sites. Unexpectedly, loss of TMG does not disrupt snRNAs localization to the nucleus, disputing an essential role of TMG in snRNA transport. Tgs1 loss also leads to defective 3' processing of snRNAs. Remarkably, stronger tgs1 mutations cause lethality without severely disrupting splicing, likely due to the preponderance of TMG-capped snRNPs. Tgs1, a predominantly nucleolar protein in Drosophila, likely carries out splicing-independent functions indispensable for animal development. Taken together, our results suggest that nuclear import is not a conserved function of TMG. As a distinctive structure on RNA, particularly non-coding RNA, we suggest that TMG prevents spurious interactions detrimental to the function of RNAs that it modifies.

CMTR1-Catalyzed 2'-O-Ribose Methylation Controls Neuronal Development by Regulating Camk2α Expression Independent of RIG-I Signaling.

Eukaryotic mRNAs are 5' end capped with a 7-methylguanosine, which is important for processing and translation of mRNAs. Cap methyltransferase 1 (CMTR1) catalyzes 2'-O-ribose methylation of the first transcribed nucleotide (N1 2'-O-Me) to mask mRNAs from innate immune surveillance by retinoic-acid-inducible gene-I (RIG-I). Nevertheless, whether this modification regulates gene expression for neuronal functions remains unexplored. Here, we find that knockdown of CMTR1 impairs dendrite development independent of secretory cytokines and RIG-I signaling. Using transcriptomic analyses, we identify altered gene expression related to dendrite morphogenesis instead of RIG-I-activated interferon signaling, such as decreased calcium/calmodulin-dependent protein kinase 2α (Camk2α). In line with these molecular changes, dendritic complexity in CMTR1-insufficient neurons is rescued by ectopic expression of CaMK2α but not by inactivation of RIG-I signaling. We further generate brain-specific CMTR1-knockout mice to validate these findings in vivo. Our study reveals the indispensable role of CMTR1-catalyzed N1 2'-O-Me in gene regulation for brain development.