Identification of tRNA nucleoside modification genes critical for stress response and development in rice and Arabidopsis

Phosphate Limitation Responses in Marine Green Algae Are Linked to Reprogramming of the tRNA Epitranscriptome and Codon Usage Bias

Marine algae are central to global carbon fixation, and their productivity is dictated largely by resource availability. Reduced nutrient availability is predicted for vast oceanic regions as an outcome of climate change; however, there is much to learn regarding response mechanisms of the tiny picoplankton that thrive in these environments, especially eukaryotic phytoplankton. Here, we investigate responses of the picoeukaryote Micromonas commoda, a green alga found throughout subtropical and tropical oceans. Under shifting phosphate availability scenarios, transcriptomic analyses revealed altered expression of transfer RNA modification enzymes and biased codon usage of transcripts more abundant during phosphate-limiting versus phosphate-replete conditions, consistent with the role of transfer RNA modifications in regulating codon recognition. To associate the observed shift in the expression of the transfer RNA modification enzyme complement with the transfer RNAs encoded by M. commoda, we also determined the transfer RNA repertoire of this alga revealing potential targets of the modification enzymes. Codon usage bias was particularly pronounced in transcripts encoding proteins with direct roles in managing phosphate limitation and photosystem-associated proteins that have ill-characterized putative functions in "light stress." The observed codon usage bias corresponds to a proposed stress response mechanism in which the interplay between stress-induced changes in transfer RNA modifications and skewed codon usage in certain essential response genes drives preferential translation of the encoded proteins. Collectively, we expose a potential underlying mechanism for achieving growth under enhanced nutrient limitation that extends beyond the catalog of up- or downregulated protein-encoding genes to the cell biological controls that underpin acclimation to changing environmental conditions.

m5U-GEPred: prediction of RNA 5-methyluridine sites based on sequence-derived and graph embedding features

5-Methyluridine (m5U) is one of the most common post-transcriptional RNA modifications, which is involved in a variety of important biological processes and disease development. The precise identification of the m5U sites allows for a better understanding of the biological processes of RNA and contributes to the discovery of new RNA functional and therapeutic targets. Here, we present m5U-GEPred, a prediction framework, to combine sequence characteristics and graph embedding-based information for m5U identification. The graph embedding approach was introduced to extract the global information of training data that complemented the local information represented by conventional sequence features, thereby enhancing the prediction performance of m5U identification. m5U-GEPred outperformed the state-of-the-art m5U predictors built on two independent species, with an average AUROC of 0.984 and 0.985 tested on human and yeast transcriptomes, respectively. To further validate the performance of our newly proposed framework, the experimentally validated m5U sites identified from Oxford Nanopore Technology (ONT) were collected as independent testing data, and in this project, m5U-GEPred achieved reasonable prediction performance with ACC of 91.84%. We hope that m5U-GEPred should make a useful computational alternative for m5U identification.

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.

Pyrimidine catabolism is required to prevent the accumulation of 5-methyluridine in RNA

5-Methylated cytosine is a frequent modification in eukaryotic RNA and DNA influencing mRNA stability and gene expression. Here we show that free 5-methylcytidine (5mC) and 5-methyl-2'-deoxycytidine are generated from nucleic acid turnover in Arabidopsis thaliana, and elucidate how these cytidines are degraded, which is unclear in eukaryotes. First CYTIDINE DEAMINASE produces 5-methyluridine (5mU) and thymidine which are subsequently hydrolyzed by NUCLEOSIDE HYDROLASE 1 (NSH1) to thymine and ribose or deoxyribose. Interestingly, far more thymine is generated from RNA than from DNA turnover, and most 5mU is directly released from RNA without a 5mC intermediate, since 5-methylated uridine (m5U) is an abundant RNA modification (m5U/U ∼1%) in Arabidopsis. We show that m5U is introduced mainly by tRNA-SPECIFIC METHYLTRANSFERASE 2A and 2B. Genetic disruption of 5mU degradation in the NSH1 mutant causes m5U to occur in mRNA and results in reduced seedling growth, which is aggravated by external 5mU supplementation, also leading to more m5U in all RNA species. Given the similarities between pyrimidine catabolism in plants, mammals and other eukaryotes, we hypothesize that the removal of 5mU is an important function of pyrimidine degradation in many organisms, which in plants serves to protect RNA from stochastic m5U modification.

A Preliminary Survey of Transfer RNA Modifications and Modifying Enzymes of the Tropical Plant Cocos nucifera L

The coconut (Cocos nucifera L.) is a commercial crop widely distributed among coastal tropical regions. It provides millions of farmers with food, fuel, cosmetics, folk medicine, and building materials. Among these, oil and palm sugar are representative extracts. However, this unique living species of Cocos has only been preliminarily studied at molecular levels. Benefiting from the genomic sequence data published in 2017 and 2021, we investigated the transfer RNA (tRNA) modifications and modifying enzymes of the coconut in this survey. An extraction method for the tRNA pool from coconut flesh was built. In total, 33 species of modified nucleosides and 66 homologous genes of modifying enzymes were confirmed using a nucleoside analysis using high-performance liquid chromatography combined with high-resolution mass spectrometry (HPLC-HRMS) and homologous protein sequence alignment. The positions of tRNA modifications, including pseudouridines, were preliminarily mapped using a oligonucleotide analysis, and the features of their modifying enzymes were summarized. Interestingly, we found that the gene encoding the modifying enzyme of 2'-O-ribosyladenosine at the 64th position of tRNA (Ar(p)64) was uniquely overexpressed under high-salinity stress. In contrast, most other tRNA-modifying enzymes were downregulated with mining transcriptomic sequencing data. According to previous physiological studies of Ar(p)64, the coconut appears to enhance the quality control of the translation process when subjected to high-salinity stress. We hope this survey can help advance research on tRNA modification and scientific studies of the coconut, as well as thinking of the safety and nutritional value of naturally modified nucleosides.

Elucidation of transfer RNAs as stress regulating agents and the experimental strategies to conceive the functional role of tRNA-derived fragments in plants

In plants, the transfer RNAs (tRNAs) exhibit their profound influence in orchestrating diverse physiological activities like cell growth, development, and response to several surrounding stimuli. The tRNAs, which were known to restrict their function solely in deciphering the codons, are now emerging as frontline defenders in stress biology. The plants that are constantly confronted with a huge panoply of stresses rely on tRNA-mediated stress regulation by altering the tRNA abundance, curbing the transport of tRNAs, fragmenting the mature tRNAs during stress. Among them, the studies on the generation of transfer RNA-derived fragments (tRFs) and their biological implication in stress response have attained huge interest. In plants, the tRFs hold stable expression patterns and regulate biological functions under diverse environmental conditions. In this review, we discuss the fate of plant tRNAs upon stress and thereafter how the tRFs are metamorphosed into sharp ammunition to wrestle with stress. We also address the various methods developed to date for uncovering the role of tRFs and their function in plants.

Ten new high-quality genome assemblies for diverse bioenergy sorghum genotypes

INTRODUCTION

Sorghum (Sorghum bicolor (L.) Moench) is an agriculturally and economically important staple crop that has immense potential as a bioenergy feedstock due to its relatively high productivity on marginal lands. To capitalize on and further improve sorghum as a potential source of sustainable biofuel, it is essential to understand the genomic mechanisms underlying complex traits related to yield, composition, and environmental adaptations.

METHODS

Expanding on a recently developed mapping population, we generated de novo genome assemblies for 10 parental genotypes from this population and identified a comprehensive set of over 24 thousand large structural variants (SVs) and over 10.5 million single nucleotide polymorphisms (SNPs).

RESULTS

We show that SVs and nonsynonymous SNPs are enriched in different gene categories, emphasizing the need for long read sequencing in crop species to identify novel variation. Furthermore, we highlight SVs and SNPs occurring in genes and pathways with known associations to critical bioenergy-related phenotypes and characterize the landscape of genetic differences between sweet and cellulosic genotypes.

DISCUSSION

These resources can be integrated into both ongoing and future mapping and trait discovery for sorghum and its myriad uses including food, feed, bioenergy, and increasingly as a carbon dioxide removal mechanism.

Epitranscriptomic mRNA modifications governing plant stress responses: underlying mechanism and potential application

Plants inevitably encounter environmental adversities, including abiotic and biotic stresses, which significantly impede plant growth and reduce crop yield. Thus, fine-tuning the fate and function of stress-responsive RNAs is indispensable for plant survival under such adverse conditions. Recently, post-transcriptional RNA modifications have been studied as a potent route to regulate plant gene expression under stress. Among over 160 mRNA modifications identified to date, N⁶ -methyladenosine (m⁶ A) in mRNAs is notable because of its multifaceted roles in plant development and stress response. Recent transcriptome-wide mapping has revealed the distribution and patterns of m⁶ A in diverse stress-responsive mRNAs in plants, building a foundation for elucidating the molecular link between m⁶ A and stress response. Moreover, the identification and characterization of m⁶ A writers, readers and erasers in Arabidopsis and other model crops have offered insights into the biological roles of m⁶ A in plant abiotic stress responses. Here, we review the recent progress of research on mRNA modifications, particularly m⁶ A, and their dynamics, distribution, regulation and biological functions in plant stress responses. Further, we posit potential strategies for breeding stress-tolerant crops by engineering mRNA modifications and propose the future direction of research on RNA modifications to gain a much deeper understanding of plant stress biology.

Exploring epitranscriptomics for crop improvement and environmental stress tolerance

Climate change and stressful environmental conditions severely hamper crop growth, development and yield. Plants respond to environmental perturbations, through their plasticity provided by key-genes, governed at post-/transcriptional levels. Gene-regulation in plants is a multilevel process controlled by diverse cellular entities that includes transcription factors (TF), epigenetic regulators and non-coding RNAs beside others. There are successful studies confirming the role of epigenetic modifications (DNA-methylation/histone-modifications) in gene expression. Recent years have witnessed emergence of a highly specialized field the "Epitranscriptomics". Epitranscriptomics deals with investigating post-transcriptional RNA chemical-modifications present across the life forms that change structural, functional and biological characters of RNA. However, deeper insights on of epitranscriptomic modifications, with >140 types known so far, are to be understood fully. Researchers have identified epitranscriptome marks (writers, erasers and readers) and mapped the site-specific RNA modifications (m6A, m⁵C, 3' uridylation, etc.) responsible for fine-tuning gene expression in plants. Simultaneous advancement in sequencing platforms, upgraded bioinformatic tools and pipelines along with conventional labelled techniques have further given a statistical picture of these epitranscriptomic modifications leading to their potential applicability in crop improvement and developing climate-smart crops. We present herein the insights on epitranscriptomic machinery in plants and how epitranscriptome and epitranscriptomic modifications underlying plant growth, development and environmental stress responses/adaptations. Third-generation sequencing technology, advanced bioinformatics tools and databases being used in plant epitranscriptomics are also discussed. Emphasis is given on potential exploration of epitranscriptome engineering for crop-improvement and developing environmental stress tolerant plants covering current status, challenges and future directions.

Epitranscriptomics: An Additional Regulatory Layer in Plants' Development and Stress Response

Epitranscriptomics has added a new layer of regulatory machinery to eukaryotes, and the advancement of sequencing technology has revealed more than 170 post-transcriptional modifications in various types of RNAs, including messenger RNA (mRNA), transfer RNA (tRNA), ribosomal RNA (rRNA), and long non-coding RNA (lncRNA). Among these, N6-methyladenosine (m⁶A) and N5-methylcytidine (m⁵C) are the most prevalent internal mRNA modifications. These regulate various aspects of RNA metabolism, mainly mRNA degradation and translation. Recent advances have shown that regulation of RNA fate mediated by these epitranscriptomic marks has pervasive effects on a plant's development and responses to various biotic and abiotic stresses. Recently, it was demonstrated that the removal of human-FTO-mediated m⁶A from transcripts in transgenic rice and potatoes caused a dramatic increase in their yield, and that the m⁶A reader protein mediates stress responses in wheat and apple, indicating that regulation of m⁶A levels could be an efficient strategy for crop improvement. However, changing the overall m⁶A levels might have unpredictable effects; therefore, the identification of precise m⁶A levels at a single-base resolution is essential. In this review, we emphasize the roles of epitranscriptomic modifications in modulating molecular, physiological, and stress responses in plants, and provide an outlook on epitranscriptome engineering as a promising tool to ensure food security by editing specific m⁶A and m⁵C sites through robust genome-editing technology.

TRM61 is essential for Arabidopsis embryo and endosperm development

Post-transcriptional modifications of tRNA molecules play crucial roles in gene expression and protein biosynthesis. Across the genera, methylation of tRNAs at N¹ of adenosine 58 (A58) by AtTRM61/AtTRM6 complex plays a critical role in maintaining the stability of initiator methionyl-tRNA (tRNA_i^Met). Recently, it was shown that mutation in AtTRM61 or AtTRM6 leads to seed abortion. However, a detailed study about the AtTRM61/AtTRM6 function in plants remains vague. Here, we found that AtTRM61 has a conserved functional structure and possesses conserved binding motifs for cofactor S-adenosyl-L-methionine (AdoMet). Mutations of the complex subunits AtTRM61/AtTRM6 result in embryo and endosperm developmental defects. The endosperm and embryo developmental defects were conditionally complemented by Attrm61-1/ + FIS2pro::AtTRM61 and Attrm61-1/ + ABI3pro::AtTRM61 indicating that AtTRM61 is required for early embryo and endosperm development. Besides, the rescue of the fertility defects in trm61/ + by overexpression of initiator tRNA suggests that AtTRM61 mutation could diminish tRNA_i^Met stability. Moreover, using yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays, we showed that AtMPK4 physically interacts with AtTRM61. The data presented here suggest that AtTRM61 has a conserved structure and function in Arabidopsis. Also, AtTRM61 may be required for tRNA_i^Met stability, embryo and endosperm development.

RNA methylation in chloroplasts or mitochondria in plants

Recent advances in our understanding of epitranscriptomic RNA methylation have expanded the complexity of gene expression regulation beyond epigenetic regulation involving DNA methylation and histone modifications. The instalment, removal, and interpretation of methylation marks on RNAs are carried out by writers (methyltransferases), erasers (demethylases), and readers (RNA-binding proteins), respectively. Contrary to an emerging body of evidence demonstrating the importance of RNA methylation in the diverse fates of RNA molecules, including splicing, export, translation, and decay in the nucleus and cytoplasm, their roles in plant organelles remain largely unclear and are only now being discovered. In particular, extremely high levels of methylation marks in chloroplast and mitochondrial RNAs suggest that RNA methylation plays essential roles in organellar biogenesis and functions in plants that are crucial for plant development and responses to environmental stimuli. Thus, unveiling the cellular components involved in RNA methylation in cell organelles is essential to better understand plant biology.

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.

RsmD, a Chloroplast rRNA m2G Methyltransferase, Plays a Role in Cold Stress Tolerance by Possibly Affecting Chloroplast Translation in Arabidopsis

Ribosomal RNA (rRNA) methylation is a pivotal process in the assembly and activity of ribosomes, which in turn play vital roles in the growth, development and stress responses of plants. Although few methyltransferases responsible for rRNA methylation have been identified in plant chloroplasts, the nature and function of these enzymes in chloroplasts remain largely unknown. In this study, we characterized ArabidopsisRsmD (At3g28460), an ortholog of the methyltransferase responsible for N2-methylguanosine (m2G) modification of 16S rRNA in Escherichia coli. Confocal microscopic analysis of an RsmD- green fluorescent protein fusion protein revealed that RsmD is localized to chloroplasts. Primer extension analysis indicated that RsmD is responsible for m2G methylation at position 915 in the 16S rRNA of Arabidopsis chloroplasts. Under cold stress, rsmd mutant plants exhibited retarded growth, i.e. had shorter roots, lower fresh weight and pale-green leaves, compared with wild-type (WT) plants. However, these phenotypes were not detected in response to drought or salt stress. Notably, the rsmd mutant was hypersensitive to erythromycin or lincomycin and accumulated fewer chloroplast proteins compared with the WT, suggesting that RsmD influences translation in chloroplasts. Complementation lines expressing RsmD in the rsmd mutant background recovered WT phenotypes. Importantly, RsmD harbored RNA methyltransferase activity. Collectively, the findings of this study indicate that RsmD is a chloroplast 16S rRNA methyltransferase responsible for m2G915 modification that plays a role in the adaptation of Arabidopsisto cold stress.

Transfer RNA-derived non-coding RNAs (tncRNAs): Hidden regulation of plants' transcriptional regulatory circuits

The emergence of distinct classes of non-coding RNAs has led to better insights into the eukaryotic gene regulatory networks. Amongst them, the existence of transfer RNA (tRNA)-derived non-coding RNAs (tncRNAs) demands exploration in the plant kingdom. We have designed a methodology to uncover the entire perspective of tncRNAome in plants. Using this pipeline, we have identified diverse tncRNAs with a size ranging from 14 to 50 nucleotides (nt) by utilizing 2448 small RNA-seq samples from six angiosperms, and studied their various features, including length, codon-usage, cleavage pattern, and modified tRNA nucleosides. Codon-dependent generation of tncRNAs suggests that the tRNA cleavage is highly specific rather than random tRNA degradation. The nucleotide composition analysis of tncRNA cleavage positions indicates that they are generated through precise endoribonucleolytic cleavage machinery. Certain nucleoside modifications detected on tncRNAs were found to be conserved across the plants, and hence may influence tRNA cleavage, as well as tncRNA functions. Pathway enrichment analysis revealed that common tncRNA targets are majorly enriched during metabolic and developmental processes. Further distinct tissue-specific tncRNA clusters highlight their role in plant development. Significant number of tncRNAs differentially expressed under abiotic and biotic stresses highlights their potential role in stress resistance. In summary, this study has developed a platform that will help in the understanding of tncRNAs and their involvement in growth, development, and response to various stresses. The workflow, software package, and results are freely available at http://nipgr.ac.in/tncRNA.

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.

Identification of the enzymes responsible for m2,2G and acp3U formation on cytosolic tRNA from insects and plants

Posttranscriptional modification of tRNA is critical for efficient protein translation and proper cell growth, and defects in tRNA modifications are often associated with human disease. Although most of the enzymes required for eukaryotic tRNA modifications are known, many of these enzymes have not been identified and characterized in several model multicellular eukaryotes. Here, we present two related approaches to identify the genes required for tRNA modifications in multicellular organisms using primer extension assays with fluorescent oligonucleotides. To demonstrate the utility of these approaches we first use expression of exogenous genes in yeast to experimentally identify two TRM1 orthologs capable of forming N2,N2-dimethylguanosine (m2,2G) on residue 26 of cytosolic tRNA in the model plant Arabidopsis thaliana. We also show that a predicted catalytic aspartate residue is required for function in each of the proteins. We next use RNA interference in cultured Drosophila melanogaster cells to identify the gene required for m2,2G26 formation on cytosolic tRNA. Additionally, using these approaches we experimentally identify D. melanogaster gene CG10050 as the corresponding ortholog of human DTWD2, which encodes the protein required for formation of 3-amino-3-propylcarboxyuridine (acp3U) on residue 20a of cytosolic tRNA. We further show that A. thaliana gene AT2G41750 can form acp3U20b on an A. thaliana tRNA expressed in yeast cells, and that the aspartate and tryptophan residues in the DXTW motif of this protein are required for modification activity. These results demonstrate that these approaches can be used to study tRNA modification enzymes.

Natural variations of SLG1 confer high-temperature tolerance in indica rice

With global warming and climate change, breeding crop plants tolerant to high-temperature stress is of immense significance. tRNA 2-thiolation is a highly conserved form of tRNA modification among living organisms. Here, we report the identification of SLG1 (Slender Guy 1), which encodes the cytosolic tRNA 2-thiolation protein 2 (RCTU2) in rice. SLG1 plays a key role in the response of rice plants to high-temperature stress at both seedling and reproductive stages. Dysfunction of SLG1 results in plants with thermosensitive phenotype, while overexpression of SLG1 enhances the tolerance of plants to high temperature. SLG1 is differentiated between the two Asian cultivated rice subspecies, indica and japonica, and the variations at both promoter and coding regions lead to an increased level of thiolated tRNA and enhanced thermotolerance of indica rice varieties. Our results demonstrate that the allelic differentiation of SLG1 confers indica rice to high-temperature tolerance, and tRNA thiolation pathway might be a potential target in the next generation rice breeding for the warming globe.

Understanding the mechanism of high-temperature tolerance will help to breed crops adaptive to warming climate. Here, the authors show SLG1, a cytosolic tRNA 2-thiolation protein 2 encoding gene, is differentiated between the two Asian cultivated rice subspecies and confers high temperature tolerance of indica rice.

m5UPred: A Web Server for the Prediction of RNA 5-Methyluridine Sites from Sequences

As one of the widely occurring RNA modifications, 5-methyluridine (m⁵U) has recently been shown to play critical roles in various biological functions and disease pathogenesis, such as under stress response and during breast cancer development. Precise identification of m⁵U sites on RNA is vital for the understanding of the regulatory mechanisms of RNA life. We present here m5UPred, the first web server for in silico identification of m⁵U sites from the primary sequences of RNA. Built upon the support vector machine (SVM) algorithm and the biochemical encoding scheme, m5UPred achieved reasonable prediction performance with the area under the receiver operating characteristic curve (AUC) greater than 0.954 by 5-fold cross-validation and independent testing datasets. To critically test and validate the performance of our newly proposed predictor, the experimentally validated m⁵U sites were further separated by high-throughput sequencing techniques (miCLIP-Seq and FICC-Seq) and cell types (HEK293 and HAP1). When tested on cross-technique and cross-cell-type validation using independent datasets, m5UPred achieved an average AUC of 0.922 and 0.926 under mature mRNA mode, respectively, showing reasonable accuracy and reliability. The m5UPred web server is freely accessible now and it should make a useful tool for the researchers who are interested in m⁵U RNA modification.

Correlation between gene expression levels under drought stress and synonymous codon usage in rice plant by in-silico study

We studied the correlation of synonymous codon usage (SCU) on gene expression levels under drought stress in rice. Sixty genes related to drought stress (with high, intermediate and low expression) were selected from rice meta-analysis data and various codon usage indices such as the effective number of codon usage (ENC), codon adaptation index (CAI) and relative synonymous codon usage (RSCU) were calculated. We found that in genes highly expressing under drought 1) GC content was higher, 2) ENC value was lower, 3) the preferred codons of some amino acids changed and 4) the RSCU ratio of GC-end codons relative to AT-end codons for 18 amino acids increased significantly compared with those in other genes. We introduce ARSCU as the Average ratio of RSCUs of GC-end codons to AT-end codons in each gene that could significantly separate high-expression genes under drought from low-expression genes. ARSCU is calculated using the program ARSCU-Calculator developed by our group to help predicting expression level of rice genes under drought. An index above ARSCU threshold is expected to indicate that the gene under study may belong to the "high expression group under drought". This information may be applied for codon optimization of genes for rice genetic engineering. To validate these findings, we further used 60 other genes (randomly selected subset of 43233 genes studied for their response to drought stress). ARSCU value was able to predict the level of expression at 88.33% of the cases. Using third set of 60 genes selected amongst high expressing genes not related to drought, only 31.65% of the genes showed ARSCU value of higher than the set threshold. This indicates that the phenomenon we described in this report may be unique for drought related genes. To justify the observed correlation between CUB and high expressing genes under drought, possible role of tRNA post transcriptional modification and tRFs was hypothesized as possible underlying biological mechanism.

Glyceroglycolipid Metabolism Regulations under Phosphate Starvation Revealed by Transcriptome Analysis in Synechococcus elongatus PCC 7942

Glyceroglycolipids, abundant in cyanobacteria's photosynthetic membranes, present bioactivities and pharmacological activities, and can be widely used in the pharmaceutical industry. Environmental factors could alter the contents and compositions of cyanobacteria glyceroglycolipids, but the regulation mechanism remains unclear. Therefore, the glyceroglycolipids contents and the transcriptome in Synechococcus elongatus PCC 7942 were analyzed under phosphate starvation. Under phosphate starvation, the decrease of monogalactosyl diacylglycerol (MGDG) and increases of digalactosyl diacylglycerol (DGDG) and sulfoquinovosyl diacylglycerol (SQDG) led to a decrease in the MGDG/DGDG ratio, from 4:1 to 5:3, after 12 days of cultivation. However, UDP-sulfoquinovose synthase gene sqdB, and the SQDG synthase gene sqdX, were down-regulated, and the decreased MGDG/DGDG ratio was later increased back to 2:1 after 15 days of cultivation, suggesting the regulation of glyceroglycolipids on day 12 was based on the MGDG/DGDG ratio maintaining glyceroglycolipid homeostasis. There are 12 differentially expressed transcriptional regulators that could be potential candidates related to glyceroglycolipid regulation, according to the transcriptome analysis. The transcriptome analysis also suggested post-transcriptional or post-translational regulations in glyceroglycolipid synthesis. This study provides further insights into glyceroglycolipid metabolism, as well as the scientific basis for glyceroglycolipid synthesis optimization and cyanobacteria glyceroglycolipids utilization via metabolic engineering.

Identification of small RNAs during cold acclimation in Arabidopsis thaliana

BACKGROUND

Cold stress causes dynamic changes in gene expression that are partially caused by small non-coding RNAs since they regulate protein coding transcripts and act in epigenetic gene silencing pathways. Thus, a detailed analysis of transcriptional changes of small RNAs (sRNAs) belonging to all known sRNA classes such as microRNAs (miRNA) and small interfering RNA (siRNAs) in response to cold contributes to an understanding of cold-related transcriptome changes.

RESULT

We subjected A. thaliana plants to cold acclimation conditions (4 °C) and analyzed the sRNA transcriptomes after 3 h, 6 h and 2 d. We found 93 cold responsive differentially expressed miRNAs and only 14 of these were previously shown to be cold responsive. We performed miRNA target prediction for all differentially expressed miRNAs and a GO analysis revealed the overrepresentation of miRNA-targeted transcripts that code for proteins acting in transcriptional regulation. We also identified a large number of differentially expressed cis- and trans-nat-siRNAs, as well as sRNAs that are derived from long non-coding RNAs. By combining the results of sRNA and mRNA profiling with miRNA target predictions and publicly available information on transcription factors, we reconstructed a cold-specific, miRNA and transcription factor dependent gene regulatory network. We verified the validity of links in the network by testing its ability to predict target gene expression under cold acclimation.

CONCLUSION

In A. thaliana, miRNAs and sRNAs derived from cis- and trans-NAT gene pairs and sRNAs derived from lncRNAs play an important role in regulating gene expression in cold acclimation conditions. This study provides a fundamental database to deepen our knowledge and understanding of regulatory networks in cold acclimation.

The Arabidopsis TRM61/TRM6 complex is a bona fide tRNA N1-methyladenosine methyltransferase

tRNA molecules, which contain the most abundant post-transcriptional modifications, are crucial for proper gene expression and protein biosynthesis. Methylation at N1 of adenosine 58 (A58) is critical for maintaining the stability of initiator methionyl-tRNA (tRNAiMet) in bacterial, archaeal, and eukaryotic tRNAs. However, although research has been conducted in yeast and mammals, it remains unclear how A58 in plant tRNAs is modified and involved in development. In this study, we identify the nucleus-localized complex AtTRM61/AtTRM6 in Arabidopsis as tRNA m1A58 methyltransferase. Deficiency or a lack of either AtTRM61 or AtTRM6 leads to embryo arrest and seed abortion. The tRNA m1A level decreases in conditionally complemented Attrm61/LEC1pro::AtTRM61 plants and this is accompanied by reduced levels of tRNAiMet, indicating the importance of the tRNA m1A modification for tRNAiMet stability. Taken together, our results demonstrate that tRNA m1A58 modification is necessary for tRNAiMet stability and is required for embryo development in Arabidopsis.

FTSJ1 regulates tRNA 2'-O-methyladenosine modification and suppresses the malignancy of NSCLC via inhibiting DRAM1 expression

Non-small cell lung cancer (NSCLC) is the leading cause of cancer mortality worldwide. The mechanisms underlying NSCLC tumorigenesis are incompletely understood. Transfer RNA (tRNA) modification is emerging as a novel regulatory mechanism for carcinogenesis. However, the role of tRNA modification in NSCLC remains obscure. In this study, HPLC/MS assay was used to quantify tRNA modification levels in NSCLC tissues and cells. tRNA-modifying enzyme genes were identified by comparative genomics and validated by qRT-PCR analysis. The biological functions of tRNA-modifying gene in NSCLC were investigated in vitro and in vivo. The mechanisms of tRNA-modifying gene in NSCLC were explored by RNA-seq, qRT-PCR, and rescue assays. The results showed that a total of 18 types of tRNA modifications and up to seven tRNA-modifying genes were significantly downregulated in NSCLC tumor tissues compared with that in normal tissues, with the 2'-O-methyladenosine (Am) modification displaying the lowest level in tumor tissues. Loss- and gain-of-function assays revealed that the amount of Am in tRNAs was significantly associated with expression levels of FTSJ1, which was also downregulated in NSCLC tissues and cells. Upregulation of FTSJ1 inhibited proliferation, migration, and promoted apoptosis of NSCLC cells in vitro. Silencing of FTSJ1 resulted in the opposite effects. In vivo assay confirmed that overexpression of FTSJ1 significantly suppressed the growth of NSCLC cells. Mechanistically, overexpression of FTSJ1 led to a decreased expression of DRAM1. Whereas knockdown of FTSJ1 resulted in an increased expression of DRAM1. Furthermore, silencing of DRAM1 substantially augmented the antitumor effect of FTSJ1 on NSCLC cells. Our findings suggested an important mechanism of tRNA modifications in NSCLC and demonstrated novel roles of FTSJ1 as both tRNA Am modifier and tumor suppressor in NSCLC.

Dynamic gene expression changes in response to micronutrient, macronutrient, and multiple stress exposures in soybean

Preserving crop yield is critical for US soybean production and the global economy. Crop species have been selected for increased yield for thousands of years with individual lines selected for improved performance in unique environments, constraints not experienced by model species such as Arabidopsis. This selection likely resulted in novel stress adaptations, unique to crop species. Given that iron deficiency is a perennial problem in the soybean growing regions of the USA and phosphate deficiency looms as a limitation to global agricultural production, nutrient stress studies in crop species are critically important. In this study, we directly compared whole-genome expression responses of leaves and roots to iron (Fe) and phosphate (Pi) deficiency, representing a micronutrient and macronutrient, respectively. Conducting experiments side by side, we observed soybean responds to both nutrient deficiencies within 24 h. While soybean responds largely to -Fe deficiency, it responds strongly to Pi resupply. Though the timing of the responses was different, both nutrient stress signals used the same molecular pathways. Our study is the first to demonstrate the speed and diversity of the soybean stress response to multiple nutrient deficiencies. We also designed the study to examine gene expression changes in response to multiple stress events. We identified 865 and 3375 genes that either altered their direction of expression after a second stress exposure or were only differentially expressed after a second stress event. Understanding the molecular underpinnings of these responses in crop species could have major implications for improving stress tolerance and preserving yield.

Arabidopsis TRM5 encodes a nuclear-localised bifunctional tRNA guanine and inosine-N1-methyltransferase that is important for growth

Modified nucleosides in tRNAs are critical for protein translation. N1-methylguanosine-37 and N1-methylinosine-37 in tRNAs, both located at the 3'-adjacent to the anticodon, are formed by Trm5. Here we describe Arabidopsis thaliana AtTRM5 (At3g56120) as a Trm5 ortholog. Attrm5 mutant plants have overall slower growth as observed by slower leaf initiation rate, delayed flowering and reduced primary root length. In Attrm5 mutants, mRNAs of flowering time genes are less abundant and correlated with delayed flowering. We show that AtTRM5 complements the yeast trm5 mutant, and in vitro methylates tRNA guanosine-37 to produce N1-methylguanosine (m1G). We also show in vitro that AtTRM5 methylates tRNA inosine-37 to produce N1-methylinosine (m1I) and in Attrm5 mutant plants, we show a reduction of both N1-methylguanosine and N1-methylinosine. We also show that AtTRM5 is localized to the nucleus in plant cells. Proteomics data showed that photosynthetic protein abundance is affected in Attrm5 mutant plants. Finally, we show tRNA-Ala aminoacylation is not affected in Attrm5 mutants. However the abundance of tRNA-Ala and tRNA-Asp 5' half cleavage products are deduced. Our findings highlight the bifunctionality of AtTRM5 and the importance of the post-transcriptional tRNA modifications m1G and m1I at tRNA position 37 in general plant growth and development.

tRNA Wobble Modification Affects Leaf Cell Development in Arabidopsis thaliana

The tRNA modification at the wobble position of Lys, Glu and Gln (wobbleU* modification) is responsible for the fine-tuning of protein translation efficiency and translation rate. This modification influences organism function in accordance with growth and environmental changes. However, the effects of wobbleU* modification at the cellular, tissue, or individual level have not yet been elucidated. In this study, we show that sulfur modification of wobbleU* of the tRNAs affects leaf development in Arabidopsis thaliana. The sulfur modification was impaired in the two wobbleU*-modification mutants: the URM1-like protein-defective mutant and the Elongator complex-defective mutants. Analyses of the mutant phenotypes revealed that the deficiency in the wobbleU* modification increased the airspaces in the leaves and the leaf size without affecting the number and the area of palisade mesophyll cells. On the other hand, both mutants exhibited increased number of leaf epidermal pavement cells but with reduced cell size. The deficiency in the wobbleU* modification also delayed the initiation of the endoreduplication processes of mesophyll cells. The phenotype of ASYMMETRIC LEAVES2-defective mutant was enhanced in the Elongator-defective mutants, while it was unchanged in the URM1-like protein-defective mutant. Collectively, the findings of this study suggest that the tRNA wobbleU* modification plays an important role in leaf morphogenesis by balancing the development between epidermal and mesophyll tissues.

Translational Regulation of Plant Response to High Temperature by a Dual-Function tRNAHis Guanylyltransferase in Rice

As sessile organisms, plants have evolved numerous strategies to acclimate to changes in environmental temperature. However, the molecular basis of this acclimation remains largely unclear. In this study we identified a tRNA^His guanylyltransferase, AET1, which contributes to the modification of pre-tRNA^His and is required for normal growth under high-temperature conditions in rice. Interestingly, AET1 possibly interacts with both RACK1A and eIF3h in the endoplasmic reticulum. Notably, AET1 can directly bind to OsARF mRNAs including the uORFs of OsARF19 and OsARF23, indicating that AET1 is associated with translation regulation. Furthermore, polysome profiling assays suggest that the translational status remains unaffected in the aet1 mutant, but that the translational efficiency of OsARF19 and OsARF23 is reduced; moreover, OsARF23 protein levels are obviously decreased in the aet1 mutant under high temperature, implying that AET1 regulates auxin signaling in response to high temperature. Our findings provide new insights into the molecular mechanisms whereby AET1 regulates the environmental temperature response in rice by playing a dual role in tRNA modification and translational control.

Messenger RNA Modifications in Plants

Over 160 distinct RNA modifications are known and collectively termed the epitranscriptome. Some of these modifications have been discovered in mRNA, uncovering a new layer of gene regulation. Transcriptome-wide mapping of epitranscriptomic codes and the discovery of their writers, erasers, and readers that dynamically install, remove, and interpret RNA modifications, respectively, are fundamental to understanding the epitranscriptome. Recent technologies have enabled the transcriptome-wide profiling of several mRNA modifications in Arabidopsis thaliana, providing key insights into regulating these modifications and their effects on plant development. Here we review technological innovations and recent progress in epitranscriptomics, with specific focus on N⁶-methyladenosine (m⁶A), 5-methylcytosine (m⁵C), uridylation, and their roles in multiple aspects of plant development.

Whole-Genome Resequencing of Seven Eggplant (Solanum melongena) and One Wild Relative (S. incanum) Accessions Provides New Insights and Breeding Tools for Eggplant Enhancement

Whole-genome resequencing provides information of great relevance for crop genetics, evolution, and breeding. Here, we present the first whole-genome resequencing study using seven eggplant (Solanum melongena) and one wild relative (Solanum incanum) accessions. These eight accessions were selected for displaying a high phenotypic and genetic diversity and for being the founder parents of an eggplant multiparent advanced generation intercrosses population. By resequencing at an average depth of 19.8× and comparing to the high-quality reference genome "67/3" over 10 million high-reliable polymorphisms were discovered, of which over 9 million (84.5%) were single nucleotide polymorphisms and more than 700,000 (6.5%) InDels. However, while for the S. melongena accessions, the variants identified ranged from 0.8 to 1.3 million, over 9 million were detected for the wild S. incanum. This confirms the narrow genetic diversity of the domesticated eggplant and suggests that introgression breeding using wild relatives can efficiently contribute to broadening the genetic basis of this crop. Differences were observed among accessions for the enrichment in genes regulating important biological processes. By analyzing the distribution of the variants, we identified the potential footprints of old introgressions from wild relatives that can help to unravel the unclear domestication and breeding history. The comprehensive annotation of these eight genomes and the information provided in this study represents a landmark in eggplant genomics and allows the development of tools for eggplant genetics and breeding.

Correction to: Identification of tRNA nucleoside modification genes critical for stress response and development in rice and Arabidopsis

Following publication of the original article [1], the authors reported that there was a mistake in the presentation of their funding information. The sentence "This study was supported by the National Natural Science Foundation of China (31,100,268 to Peng Chen, 31,270,658 to Bo Zheng);" should instead read "This study was supported by the National Natural Science Foundation of China (31100268 to Peng Chen, 31370604 to Bo Zheng);".

Fine-Tuning of Gene Expression by tRNA-Derived Fragments during Abiotic Stress Signal Transduction

When plants are subjected to unfavorable environmental conditions, overall gene expression in stressed cells is altered from a programmed pattern for normal development to an adaptive pattern for survival. Rapid changes in plant gene expression include production of stress responsive proteins for protection as well as reduction of irrelevant proteins to minimize energy consumption during growth. In addition to the many established mechanisms known to modulate gene expression in eukaryotes, a novel strategy involving tRNA-derived fragments (tRFs) was recently reported to control gene expression. In animals, tRFs are shown to play a certain role in infected or cancer cells. However, tRFs are expected to function in the regulation of gene expression against abiotic stress conditions in plants. Moreover, the underlying mechanism linking up-regulation of tRFs under stress conditions with the stress tolerant response remains unknown. In this review, the biogenesis and putative function of diverse tRFs in abiotic stress signaling are discussed with a focus on tRFs as a transcriptional/post-transcriptional/translational regulator.