Global Epitranscriptomics Profiling of RNA Post-Transcriptional Modifications as an Effective Tool for Investigating the Epitranscriptomics of Stress Response

Differential Protein Expression in Set5p-Mediated Acetic Acid Stress Response and Novel Targets for Engineering Yeast Stress Tolerance

Acetic acid is a prevalent inhibitor in lignocellulosic hydrolysate, which represses microbial growth and bioproduction. Histone modification and chromatin remodeling have been revealed to be critical for regulating eukaryotic metabolism. However, related studies in chronic acetic acid stress responses remain unclear. Our previous studies revealed that overexpression of the histone H4 methyltransferase Set5p enhanced acetic acid stress tolerance of the budding yeast Saccharomyces cerevisiae. In this study, we examined the role of Set5p in acetic acid stress by analyzing global protein expression. Significant activation of intracellular protein expression under the stress was discovered, and the functions of the differential proteins were mainly involved in chromatin modification, signal transduction, and carbohydrate metabolism. Notably, a substantial increase of Set5p expression was observed in response to acetic acid stress. Functional studies demonstrated that the restriction of the telomere capping protein Rtc3p, as well as Ies3p and Taf14p, which are related to chromatin regulation, was critical for yeast stress response. This study enriches the understanding of the epigenetic regulatory mechanisms underlying yeast stress response mediated by histone-modifying enzymes. The results also benefit the development of robust yeast strains for lignocellulosic bioconversion.

Analyzing RNA posttranscriptional modifications to decipher the epitranscriptomic code

The discovery of RNA silencing has revealed that non-protein-coding sequences (ncRNAs) can cover essential roles in regulatory networks and their malfunction may result in severe consequences on human health. These findings have prompted a general reassessment of the significance of RNA as a key player in cellular processes. This reassessment, however, will not be complete without a greater understanding of the distribution and function of the over 170 variants of the canonical ribonucleotides, which contribute to the breathtaking structural diversity of natural RNA. This review surveys the analytical approaches employed for the identification, characterization, and detection of RNA posttranscriptional modifications (rPTMs). The merits of analyzing individual units after exhaustive hydrolysis of the initial biopolymer are outlined together with those of identifying their position in the sequence of parent strands. Approaches based on next generation sequencing and mass spectrometry technologies are covered in depth to provide a comprehensive view of their respective merits. Deciphering the epitranscriptomic code will require not only mapping the location of rPTMs in the various classes of RNAs, but also assessing the variations of expression levels under different experimental conditions. The fact that no individual platform is currently capable of meeting all such demands implies that it will be essential to capitalize on complementary approaches to obtain the desired information. For this reason, the review strived to cover the broadest possible range of techniques to provide readers with the fundamental elements necessary to make informed choices and design the most effective possible strategy to accomplish the task at hand.

The emerging role of epitranscriptome in shaping stress responses in plants

KEY MESSAGE

RNA modifications and editing changes constitute 'epitranscriptome' and are crucial in regulating the development and stress response in plants. Exploration of the epitranscriptome and associated machinery would facilitate the engineering of stress tolerance in crops. RNA editing and modifications post-transcriptionally decorate almost all classes of cellular RNAs, including tRNAs, rRNAs, snRNAs, lncRNAs and mRNAs, with more than 170 known modifications, among which m6A, Ψ, m5C, 8-OHG and C-to-U editing are the most abundant. Together, these modifications constitute the "epitranscriptome", and contribute to changes in several RNA attributes, thus providing an additional structural and functional diversification to the "cellular messages" and adding another layer of gene regulation in organisms, including plants. Numerous evidences suggest that RNA modifications have a widespread impact on plant development as well as in regulating the response of plants to abiotic and biotic stresses. High-throughput sequencing studies demonstrate that the landscapes of m6A, m5C, Am, Cm, C-to-U, U-to-G, and A-to-I editing are remarkably dynamic during stress conditions in plants. GO analysis of transcripts enriched in Ψ, m6A and m5C modifications have identified bonafide components of stress regulatory pathways. Furthermore, significant alterations in the expression pattern of genes encoding writers, readers, and erasers of certain modifications have been documented when plants are grown in challenging environments. Notably, manipulating the expression levels of a few components of RNA editing machinery markedly influenced the stress tolerance in plants. We provide updated information on the current understanding on the contribution of RNA modifications in shaping the stress responses in plants. Unraveling of the epitranscriptome has opened new avenues for designing crops with enhanced productivity and stress resilience in view of global climate change.

Changes in m6A in Steatotic Liver Disease

Fatty liver disease is one of the major causes of morbidity and mortality worldwide. Fatty liver includes non-alcoholic fatty liver disease (NAFLD) and non-alcoholic steatohepatitis (NASH), now replaced by a consensus group as metabolic dysfunction-associated steatotic liver disease (MASLD). While excess nutrition and obesity are major contributors to fatty liver, the underlying mechanisms remain largely unknown and therapeutic interventions are limited. Reversible chemical modifications in RNA are newly recognized critical regulators controlling post-transcriptional gene expression. Among these modifications, N6-methyladenosine (m6A) is the most abundant and regulates transcript abundance in fatty liver disease. Modulation of m6A by readers, writers, and erasers (RWE) impacts mRNA processing, translation, nuclear export, localization, and degradation. While many studies focus on m6A RWE expression in human liver pathologies, limitations of technology and bioinformatic methods to detect m6A present challenges in understanding the epitranscriptomic mechanisms driving fatty liver disease progression. In this review, we summarize the RWE of m6A and current methods of detecting m6A in specific genes associated with fatty liver disease.

H/ACA snRNP-dependent ribosome biogenesis regulates translation of polyglutamine proteins

Stem cells in many systems, including Drosophila germline stem cells (GSCs), increase ribosome biogenesis and translation during terminal differentiation. Here, we show that the H/ACA small nuclear ribonucleoprotein (snRNP) complex that promotes pseudouridylation of ribosomal RNA (rRNA) and ribosome biogenesis is required for oocyte specification. Reducing ribosome levels during differentiation decreased the translation of a subset of messenger RNAs that are enriched for CAG trinucleotide repeats and encode polyglutamine-containing proteins, including differentiation factors such as RNA-binding Fox protein 1. Moreover, ribosomes were enriched at CAG repeats within transcripts during oogenesis. Increasing target of rapamycin (TOR) activity to elevate ribosome levels in H/ACA snRNP complex-depleted germlines suppressed the GSC differentiation defects, whereas germlines treated with the TOR inhibitor rapamycin had reduced levels of polyglutamine-containing proteins. Thus, ribosome biogenesis and ribosome levels can control stem cell differentiation via selective translation of CAG repeat-containing transcripts.

Disruption of the mouse liver epitranscriptome by long-term aroclor 1260 exposure

Chronic environmental exposure to polychlorinated biphenyls (PCBs) is associated with non-alcoholic fatty liver disease (NAFLD) and exacerbated by a high fat diet (HFD). Here, chronic (34 wks.) exposure of low fat diet (LFD)-fed male mice to Aroclor 1260 (Ar1260), a non-dioxin-like (NDL) mixture of PCBs, resulted in steatohepatitis and NAFLD. Twelve hepatic RNA modifications were altered with Ar1260 exposure including reduced abundance of 2'-O-methyladenosine (Am) and N(6)-methyladenosine (m6A), in contrast to increased Am in the livers of HFD-fed, Ar1260-exposed mice reported previously. Differences in 13 RNA modifications between LFD- and HFD- fed mice, suggest that diet regulates the liver epitranscriptome. Integrated network analysis of epitranscriptomic modifications identified a NRF2 (Nfe2l2) pathway in the chronic, LFD, Ar1260-exposed livers and an NFATC4 (Nfatc4) pathway for LFD- vs. HFD-fed mice. Changes in protein abundance were validated. The results demonstrate that diet and Ar1260 exposure alter the liver epitranscriptome in pathways associated with NAFLD.

Data-Independent Acquisition for the Detection of Mononucleoside RNA Modifications by Mass Spectrometry

RNA is dynamically modified in cells by a plethora of chemical moieties to modulate molecular functions and processes. Over 140 modifications have been identified across species and RNA types, with the highest density and diversity of modifications found in tRNA (tRNA). The methods used to identify and quantify these modifications have developed over recent years and continue to advance, primarily in the fields of next-generation sequencing (NGS) and mass spectrometry (MS). Most current NGS methods are limited to antibody-recognized or chemically derivatized modifications and have limitations in identifying multiple modifications simultaneously. Mass spectrometry can overcome both of these issues, accurately identifying a large number of modifications in a single run. Here, we present advances in MS data acquisition for the purpose of RNA modification identification and quantitation. Using this approach, we identified multiple tRNA wobble position modifications in Arabidopsis thaliana that are upregulated in salt-stressed growth conditions and may stabilize translation of salt stress induced proteins. This work presents improvements in methods for studying RNA modifications and introduces a possible regulatory role of wobble position modifications in A. thaliana translation.

Epitranscriptomic Reprogramming Is Required to Prevent Stress and Damage from Acetaminophen

Epitranscriptomic marks, in the form of enzyme catalyzed RNA modifications, play important gene regulatory roles in response to environmental and physiological conditions. However, little is known with respect to how acute toxic doses of pharmaceuticals influence the epitranscriptome. Here we define how acetaminophen (APAP) induces epitranscriptomic reprogramming and how the writer Alkylation Repair Homolog 8 (Alkbh8) plays a key gene regulatory role in the response. Alkbh8 modifies tRNA selenocysteine (tRNA^Sec) to translationally regulate the production of glutathione peroxidases (Gpx’s) and other selenoproteins, with Gpx enzymes known to play protective roles during APAP toxicity. We demonstrate that APAP increases toxicity and markers of damage, and decreases selenoprotein levels in Alkbh8 deficient mouse livers, when compared to wildtype. APAP also promotes large scale reprogramming of many RNA marks comprising the liver tRNA epitranscriptome including: 5-methoxycarbonylmethyluridine (mcm⁵U), isopentenyladenosine (i⁶A), pseudouridine (Ψ), and 1-methyladenosine (m¹A) modifications linked to tRNA^Sec and many other tRNA’s. Alkbh8 deficiency also leads to wide-spread epitranscriptomic dysregulation in response to APAP, demonstrating that a single writer defect can promote downstream changes to a large spectrum of RNA modifications. Our study highlights the importance of RNA modifications and translational responses to APAP, identifies writers as key modulators of stress responses in vivo and supports the idea that the epitranscriptome may play important roles in responses to pharmaceuticals.

Variations in transfer and ribosomal RNA epitranscriptomic status can adapt eukaryote translation to changing physiological and environmental conditions

The timely reprogramming of gene expression in response to internal and external cues is essential to eukaryote development and acclimation to changing environments. Chemically modifying molecular receptors and transducers of these signals is one way to efficiently induce proper physiological responses. Post-translation modifications, regulating protein biological activities, are central to many well-known signal-responding pathways. Recently, messenger RNA (mRNA) chemical (i.e. epitranscriptomic) modifications were also shown to play a key role in these processes. In contrast, transfer RNA (tRNA) and ribosomal RNA (rRNA) chemical modifications, although critical for optimal function of the translation apparatus, and much more diverse and quantitatively important compared to mRNA modifications, were until recently considered as mainly static chemical decorations. We present here recent observations that are challenging this view and supporting the hypothesis that tRNA and rRNA modifications dynamically respond to various cell and environmental conditions and contribute to adapt translation to these conditions.

Agnostic Framework for the Classification/Identification of Organisms Based on RNA Post-Transcriptional Modifications

We propose a novel approach for building a classification/identification framework based on the full complement of RNA post-transcriptional modifications (rPTMs) expressed by an organism at basal conditions. The approach relies on advanced mass spectrometry techniques to characterize the products of exonuclease digestion of total RNA extracts. Sample profiles comprising identities and relative abundances of all detected rPTM were used to train and test the capabilities of different machine learning (ML) algorithms. Each algorithm proved capable of identifying rigorous decision rules for differentiating closely related classes and correctly assigning unlabeled samples. The ML classifiers resolved different members of the Enterobacteriaceae family, alternative Escherichia coli serotypes, a series of Saccharomyces cerevisiae knockout mutants, and primary cells of the Homo sapiens central nervous system, which shared very similar genetic backgrounds. The excellent levels of accuracy and resolving power achieved by training on a limited number of classes were successfully replicated when the number of classes was significantly increased to escalate complexity. A dendrogram generated from ML-curated data exhibited a hierarchical organization that closely resembled those afforded by established taxonomic systems. Finer clustering patterns revealed the extensive effects induced by the deletion of a single pivotal gene. This information provided a putative roadmap for exploring the roles of rPTMs in their respective regulatory networks, which will be essential to decipher the epitranscriptomics code. The ubiquitous presence of RNA in virtually all living organisms promises to enable the broadest possible range of applications, with significant implications in the diagnosis of RNA-related diseases.

Detecting the epitranscriptome

RNA modifications and their corresponding epitranscriptomic writer and eraser enzymes regulate gene expression. Altered RNA modification levels, dysregulated writers, and sequence changes that disrupt epitranscriptomic marks have been linked to mitochondrial and neurological diseases, cancer, and multifactorial disorders. The detection of epitranscriptomics marks is challenging, but different next generation sequencing (NGS)-based and mass spectrometry-based approaches have been used to identify and quantitate the levels of individual and groups of RNA modifications. NGS and mass spectrometry-based approaches have been coupled with chemical, antibody or enzymatic methodologies to identify modifications in most RNA species, mapped sequence contexts and demonstrated the dynamics of specific RNA modifications, as well as the collective epitranscriptome. While epitranscriptomic analysis is currently limited to basic research applications, specific approaches for the detection of individual RNA modifications and the epitranscriptome have potential biomarker applications in detecting human conditions and diseases. This article is categorized under: RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems RNA Processing > tRNA Processing RNA in Disease and Development > RNA in Disease.

Methylation of rRNA as a host defense against rampant group II intron retrotransposition

BACKGROUND

Group II introns are mobile retroelements, capable of invading new sites in DNA. They are self-splicing ribozymes that complex with an intron-encoded protein to form a ribonucleoprotein that targets DNA after splicing. These molecules can invade DNA site-specifically, through a process known as retrohoming, or can invade ectopic sites through retrotransposition. Retrotransposition, in particular, can be strongly influenced by both environmental and cellular factors.

RESULTS

To investigate host factors that influence retrotransposition, we performed random insertional mutagenesis using the ISS1 transposon to generate a library of over 1000 mutants in Lactococcus lactis, the native host of the Ll.LtrB group II intron. By screening this library, we identified 92 mutants with increased retrotransposition frequencies (RTP-ups). We found that mutations in amino acid transport and metabolism tended to have increased retrotransposition frequencies. We further explored a subset of these RTP-up mutants, the most striking of which is a mutant in the ribosomal RNA methyltransferase rlmH, which exhibited a reproducible 20-fold increase in retrotransposition frequency. In vitro and in vivo experiments revealed that ribosomes in the rlmH mutant were defective in the m3Ψ modification and exhibited reduced binding to the intron RNA.

CONCLUSIONS

Taken together, our results reinforce the importance of the native host organism in regulating group II intron retrotransposition. In particular, the evidence from the rlmH mutant suggests a role for ribosome modification in limiting rampant retrotransposition.

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.

2'-O-ribose methylation of transfer RNA promotes recovery from oxidative stress in Saccharomyces cerevisiae

Chemical modifications that regulate protein expression at the translational level are emerging as vital components of the cellular stress response. Transfer RNAs (tRNAs) are significant targets for methyl-based modifications, which are catalyzed by tRNA methyltransferases (Trms). Here, Saccharomyces cerevisiae served as a model eukaryote system to investigate the role of 2'-O-ribose tRNA methylation in the cell's response to oxidative stress. Using 2'-O-ribose deletion mutants for trms 3, 7, 13, and 44, in acute and chronic exposure settings, we demonstrate a broad cell sensitivity to oxidative stress-inducing toxicants (i.e., hydrogen peroxide, rotenone, and acetic acid). A global analysis of hydrogen peroxide-induced tRNA modifications shows a complex profile of decreased, or undetectable, 2'-O-ribose modification events in 2’-O-ribose trm mutant strains, providing a critical link between this type of modification event and Trm status post-exposure. Based on the pronounced oxidative stress sensitivity observed for trm7 mutants, we used a bioinformatic tool to identify transcripts as candidates for regulation by Trm7-catalyzed modifications (i.e., enriched in UUC codons decoded by tRNA^PheGmAA). This screen identified transcripts linked to diverse biological processes that promote cellular recovery after oxidative stress exposure, including DNA repair, chromatin remodeling, and nutrient acquisition (i.e., CRT10, HIR3, HXT2, and GNP1); moreover, these mutants were also oxidative stress-sensitive. Together, these results solidify a role for TRM3, 7, 13, and 44, in the cellular response to oxidative stress, and implicate 2'-O-ribose tRNA modification as an epitranscriptomic strategy for oxidative stress recovery.

High-resolution ion mobility spectrometry-mass spectrometry of isomeric/isobaric ribonucleotide variants

In this report, we explored the benefits of cyclic ion mobility (cIM) mass spectrometry in the analysis of isomeric post-transcriptional modifications of RNA. Standard methyl-cytidine samples were initially utilized to test the ability to correctly distinguish different structures sharing the same elemental composition and thus molecular mass. Analyzed individually, the analytes displayed characteristic arrival times (tD ) determined by the different positions of the modifying methyl groups onto the common cytidine scaffold. Analyzed in mixture, the widths of the respective signals resulted in significant overlap that initially prevented their resolution on the tD scale. The separation of the four isomers was achieved by increasing the number of passes through the cIM device, which enabled to fully differentiate the characteristic ion mobility behaviors associated with very subtle structural variations. The placement of the cIM device between the mass-selective quadrupole and the time-of-flight analyzer allowed us to perform gas-phase activation of each of these ion populations, which had been first isolated according to a common mass-to-charge ratio and then separated on the basis of different ion mobility behaviors. The observed fragmentation patterns confirmed the structures of the various isomers thus substantiating the benefits of complementing unique tD information with specific fragmentation data to reach more stringent analyte identification. These capabilities were further tested by analyzing natural mono-nucleotide mixtures obtained by exonuclease digestion of total RNA extracts. In particular, the combination of cIM separation and post-mobility dissociation allowed us to establish the composition of methyl-cytidine and methyl-adenine components present in the entire transcriptome of HeLa cells. For this reason, we expect that this technique will benefit not only epitranscriptomic studies requiring the determination of identity and expression levels of RNA modifications, but also metabolomics investigations involving the analysis of natural extracts that may possibly contain subsets of isomeric/isobaric species.

Positive-sense RNA viruses reveal the complexity and dynamics of the cellular and viral epitranscriptomes during infection

More than 140 post-transcriptional modifications (PTMs) are known to decorate cellular RNAs, but their incidence, identity and significance in viral RNA are still largely unknown. We have developed an agnostic analytical approach to comprehensively survey PTMs on viral and cellular RNAs. Specifically, we used mass spectrometry to analyze PTMs on total RNA isolated from cells infected with Zika virus, Dengue virus, hepatitis C virus (HCV), poliovirus and human immunodeficiency virus type 1. All five RNA viruses significantly altered global PTM landscapes. Examination of PTM profiles of individual viral genomes isolated by affinity capture revealed a plethora of PTMs on viral RNAs, which far exceeds the handful of well-characterized modifications. Direct comparison of viral epitranscriptomes identified common and virus-specific PTMs. In particular, specific dimethylcytosine modifications were only present in total RNA from virus-infected cells, and in intracellular HCV RNA, and viral RNA from Zika and HCV virions. Moreover, dimethylcytosine abundance during viral infection was modulated by the cellular DEAD-box RNA helicase DDX6. By opening the Pandora's box on viral PTMs, this report presents numerous questions and hypotheses on PTM function and strongly supports PTMs as a new tier of regulation by which RNA viruses subvert the host and evade cellular surveillance systems.

The Versatile Roles of the tRNA Epitranscriptome during Cellular Responses to Toxic Exposures and Environmental Stress

Living organisms respond to environmental changes and xenobiotic exposures by regulating gene expression. While heat shock, unfolded protein, and DNA damage stress responses are well-studied at the levels of the transcriptome and proteome, tRNA-mediated mechanisms are only recently emerging as important modulators of cellular stress responses. Regulation of the stress response by tRNA shows a high functional diversity, ranging from the control of tRNA maturation and translation initiation, to translational enhancement through modification-mediated codon-biased translation of mRNAs encoding stress response proteins, and translational repression by stress-induced tRNA fragments. tRNAs need to be heavily modified post-transcriptionally for full activity, and it is becoming increasingly clear that many aspects of tRNA metabolism and function are regulated through the dynamic introduction and removal of modifications. This review will discuss the many ways that nucleoside modifications confer high functional diversity to tRNAs, with a focus on tRNA modification-mediated regulation of the eukaryotic response to environmental stress and toxicant exposures. Additionally, the potential applications of tRNA modification biology in the development of early biomarkers of pathology will be highlighted.

Lifestyle modifications: coordinating the tRNA epitranscriptome with codon bias to adapt translation during stress responses

Cells adapt to stress by altering gene expression at multiple levels. Here, we propose a new mechanism regulating stress-dependent gene expression at the level of translation, with coordinated interplay between the tRNA epitranscriptome and biased codon usage in families of stress-response genes. In this model, auxiliary genetic information contained in synonymous codon usage enables regulation of codon-biased and functionally related transcripts by dynamic changes in the tRNA epitranscriptome. This model partly explains the association between synchronous stress-dependent epitranscriptomic marks and significant multi-codon usage skewing in families of translationally regulated transcripts. The model also predicts translational adaptation during viral infections.

The Profile and Dynamics of RNA Modifications in Animals

More than a hundred distinct modified nucleosides have been identified in RNA, but little is known about their distribution across different organisms, their dynamic nature and their response to cellular and environmental stress. Mass-spectrometry-based methods have been at the forefront of identifying and quantifying modified nucleosides. However, they often require synthetic reference standards, which do not exist in the case of many modified nucleosides, and this therefore impedes their analysis. Here we use a metabolic labelling approach to achieve rapid generation of bio-isotopologues of the complete Caenorhabditis elegans transcriptome and its modifications and use them as reference standards to characterise the RNA modification profile in this multicellular organism through an untargeted liquid-chromatography tandem high-resolution mass spectrometry (LC-HRMS) approach. We furthermore show that several of these RNA modifications have a dynamic response to environmental stress and that, in particular, changes in the tRNA wobble base modification 5-methoxycarbonylmethyl-2-thiouridine (mcm5 s2 U) lead to codon-biased gene-expression changes in starved animals.

Defective Ribonucleoproteins, Mistakes in RNA Processing, and Diseases

Ribonucleoproteins (RNPs) are vital to many cellular events. To this end, many neurodegenerative diseases and cancers have been linked to RNP malfunction, particularly as this relates to defective processing of cellular RNA. The connection of RNPs and diseases has also propagated a shift of focus onto RNA targeting from traditional protein targeting treatments. However, therapeutic development in this area has been limited by incomplete molecular insight into the specific contributions of RNPs to disease. This review outlines the role of several RNPs in diseases, focusing on molecular defects in processes that affect proper RNA handling in the cell. This work also evaluates the contributions of recently developed methods to understanding RNP association and function. We review progress in this area by focusing on molecular malfunctions of RNPs associated with the onset and progression of several neurodegenerative diseases and cancer and conclude with a brief discussion of RNA-based therapeutic efforts.