Combining tRNA sequencing methods to characterize plant tRNA expression and post-transcriptional modification

Large-scale duplication events underpin population-level flexibility in tRNA gene copy number in Pseudomonas fluorescens SBW25

The complement of tRNA genes within a genome is typically considered to be a (relatively) stable characteristic of an organism. Here, we demonstrate that bacterial tRNA gene set composition can be more flexible than previously appreciated, particularly regarding tRNA gene copy number. We report the high-rate occurrence of spontaneous, large-scale, tandem duplication events in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The identified duplications are up to ∼1 Mb in size (∼15% of the wildtype genome) and are predicted to change the copy number of up to 917 genes, including several tRNA genes. The observed duplications are inherently unstable: they occur, and are subsequently lost, at extremely high rates. We propose that this unusually plastic type of mutation provides a mechanism by which tRNA gene set diversity can be rapidly generated, while simultaneously preserving the underlying tRNA gene set in the absence of continued selection. That is, if a tRNA set variant provides no fitness advantage, then high-rate segregation of the duplication ensures the maintenance of the original tRNA gene set. However, if a tRNA gene set variant is beneficial, the underlying duplication fragment(s) may persist for longer and provide raw material for further, more stable, evolutionary change.

A robust method for measuring aminoacylation through tRNA-Seq

Current methods to quantify the fraction of aminoacylated tRNAs, also known as the tRNA charge, are limited by issues with either low throughput, precision, and/or accuracy. Here, we present an optimized charge tRNA-Seq method that combines previous developments with newly described approaches to establish a protocol for precise and accurate tRNA charge measurements. We verify that this protocol provides robust quantification of tRNA aminoacylation and we provide an end-to-end method that scales to hundreds of samples including software for data processing. Additionally, we show that this method supports measurements of relative tRNA expression levels and can be used to infer tRNA modifications through reverse transcription misincorporations, thereby supporting multipurpose applications in tRNA biology.

Advances in methods for tRNA sequencing and quantification

In the past decade tRNA sequencing (tRNA-seq) has attracted considerable attention as an important tool for the development of novel approaches to quantify highly modified tRNA species and to propel tRNA research aimed at understanding the cellular physiology and disease and development of tRNA-based therapeutics. Many methods are available to quantify tRNA abundance while accounting for modifications and tRNA charging/acylation. Advances in both library preparation methods and bioinformatic workflows have enabled developments in next-generation sequencing (NGS) workflows. Other approaches forgo NGS applications in favor of hybridization-based approaches. In this review we provide a brief comparative overview of various tRNA quantification approaches, focusing on the advantages and disadvantages of these methods, which together facilitate reliable tRNA quantification.

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.

ALL-tRNAseq enables robust tRNA profiling in tissue samples

Transfer RNAs (tRNAs) are small adaptor RNAs essential for mRNA translation. Alterations in the cellular tRNA population can directly affect mRNA decoding rates and translational efficiency during cancer development and progression. To evaluate changes in the composition of the tRNA pool, multiple sequencing approaches have been developed to overcome reverse transcription blocks caused by the stable structures of these molecules and their numerous base modifications. However, it remains unclear whether current sequencing protocols faithfully capture tRNAs existing in cells or tissues. This is specifically challenging for clinical tissue samples that often present variable RNA qualities. For this reason, we developed ALL-tRNAseq, which combines the highly processive MarathonRT and RNA demethylation for the robust assessment of tRNA expression, together with a randomized adapter ligation strategy prior to reverse transcription to assess tRNA fragmentation levels in both cell lines and tissues. Incorporation of tRNA fragments not only informed on sample integrity but also significantly improved tRNA profiling of tissue samples. Our data showed that our profiling strategy effectively improves classification of oncogenic signatures in glioblastoma and diffuse large B-cell lymphoma tissues, particularly for samples presenting higher levels of RNA fragmentation, further highlighting the utility of ALL-tRNAseq for translational research.

The role of post-transcriptional modification on a new tRNAIle(GAU) identified from Ganoderma lucidum in its fragments' cytotoxicity on cancer cells

Ganoderma lucidum (Ganoderma) is a famous Chinese herbal medicine which has been used clinically for thousands of years in China. Despite numerous studies on triterpenes and polysaccharides, the bioactivity of RNAs abundant in Ganoderma remains unknown. Here, based on LC-MS techniques, dihydrouracil, 5-methyluridine (m⁵U) and pseudouridine were identified at position 19, 52 and 53 of a new tRNA^Ile(GAU) which was isolated as the most abundant tRNA species in Ganoderma, and is the first purified tRNA from fungus. Cytotoxic screening of tRNA-half (t-half) and tRNA fragment (tRF) derived from this tRNA, as well as their mimics (t-half or tRF as antisense strand), demonstrated that the double-stranded form, i.e., tRF and t-halve mimics, exhibited stronger cytotoxicity than their single-stranded form, and the cytotoxicity of t-half mimic is significantly stronger than that of tRF mimic. Notably, the cytotoxicity of 3'-t-half mimic is not only much more potent than that of taxol, but also is much more potent than that of ganoderic acids, the major bioactive components in Ganoderma. Furthermore, 3'-t-half mimic_M2 (m⁵U modified) exhibited significantly stronger cytotoxicity than unmodified 3'-t-half mimic, which is consistent with the computational simulation showing that m⁵U modification enhances the stability of the tertiary structure of 3'-t-half mimic. Overall, the present study not only indicates t-halves are bioactive components in Ganoderma which should not be neglected, but also reveals an important role of post-transcriptional modification on tRNA in its fragments' cytotoxicity against cancer cells, which benefits the design and development of RNAi drugs from natural resource.

PltRNAdb: Plant transfer RNA database

Transfer RNAs (tRNAs) are intermediate-sized non-coding RNAs found in all organisms that help translate messenger RNA into protein. Recently, the number of sequenced plant genomes has increased dramatically. The availability of this extensive data greatly accelerates the study of tRNAs on a large scale. Here, 8,768,261 scaffolds/chromosomes containing 229,093 giga-base pairs representing whole-genome sequences of 256 plant species were analyzed to identify tRNA genes. As a result, 331,242 nuclear, 3,216 chloroplast, and 1,467 mitochondrial tRNA genes were identified. The nuclear tRNA genes include 275,134 tRNAs decoding 20 standard amino acids, 1,325 suppressor tRNAs, 6,273 tRNAs with unknown isotypes, 48,475 predicted pseudogenes, and 37,873 tRNAs with introns. Efforts also extended to the creation of PltRNAdb (https://bioinformatics.um6p.ma/PltRNAdb/index.php), a data source for tRNA genes from 256 plant species. PltRNAdb website allows researchers to search, browse, visualize, BLAST, and download predicted tRNA genes. PltRNAdb will help improve our understanding of plant tRNAs and open the door to discovering the unknown regulatory roles of tRNAs in plant genomes.

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.

Advances in Non-Coding RNA Sequencing

Non-coding RNAs (ncRNAs) comprise a set of abundant and functionally diverse RNA molecules. Since the discovery of the first ncRNA in the 1960s, ncRNAs have been shown to be involved in nearly all steps of the central dogma of molecular biology. In recent years, the pace of discovery of novel ncRNAs and their cellular roles has been greatly accelerated by high-throughput sequencing. Advances in sequencing technology, library preparation protocols as well as computational biology helped to greatly expand our knowledge of which ncRNAs exist throughout the kingdoms of life. Moreover, RNA sequencing revealed crucial roles of many ncRNAs in human health and disease. In this review, we discuss the most recent methodological advancements in the rapidly evolving field of high-throughput sequencing and how it has greatly expanded our understanding of ncRNA biology across a large number of different organisms.

Rapid shifts in mitochondrial tRNA import in a plant lineage with extensive mitochondrial tRNA gene loss

In most eukaryotes, transfer RNAs (tRNAs) are one of the very few classes of genes remaining in the mitochondrial genome, but some mitochondria have lost these vestiges of their prokaryotic ancestry. Sequencing of mitogenomes from the flowering plant genus Silene previously revealed a large range in tRNA gene content, suggesting rapid and ongoing gene loss/replacement. Here, we use this system to test longstanding hypotheses about how mitochondrial tRNA genes are replaced by importing nuclear-encoded tRNAs. We traced the evolutionary history of these gene loss events by sequencing mitochondrial genomes from key outgroups (Agrostemma githago and Silene [=Lychnis] chalcedonica). We then performed the first global sequencing of purified plant mitochondrial tRNA populations to characterize the expression of mitochondrial-encoded tRNAs and the identity of imported nuclear-encoded tRNAs. We also confirmed the utility of high-throughput sequencing methods for the detection of tRNA import by sequencing mitochondrial tRNA populations in a species (Solanum tuberosum) with known tRNA trafficking patterns. Mitochondrial tRNA sequencing in Silene revealed substantial shifts in the abundance of some nuclear-encoded tRNAs in conjunction with their recent history of mt-tRNA gene loss and surprising cases where tRNAs with anticodons still encoded in the mitochondrial genome also appeared to be imported. These data suggest that nuclear-encoded counterparts are likely replacing mitochondrial tRNAs even in systems with recent mitochondrial tRNA gene loss, and the redundant import of a nuclear-encoded tRNA may provide a mechanism for functional replacement between translation systems separated by billions of years of evolutionary divergence.

Plant transfer RNA-derived fragments: Biogenesis and functions

Processing of mature transfer RNAs (tRNAs) produces complex populations of tRNA-derived fragments (tRFs). Emerging evidence shows that tRFs have important functions in bacteria, animals, and plants. Here, we review recent advances in understanding plant tRFs, focusing on their biological and cellular functions, such as regulating stress responses, mediating plant-pathogen interactions, and modulating post-transcriptional gene silencing and translation. We also review sequencing strategies and bioinformatics resources for studying tRFs in plants. Finally, we discuss future directions for plant tRF research, which will expand our knowledge of plant non-coding RNAs.

tRNA-derived fragments: A new class of non-coding RNA with key roles in nervous system function and dysfunction

tRNA-derived small RNAs (tsRNA) are a recently identified family of non-coding RNA that have been associated with a variety of cellular functions including the regulation of protein translation and gene expression. Recent sequencing and bioinformatic studies have identified the broad spectrum of tsRNA in the nervous system and demonstrated that this new class of non-coding RNA is produced from tRNA by specific cleavage events catalysed by ribonucleases such as angiogenin and dicer. Evidence is also accumulating that production of tsRNA is increased during disease processes where they regulate stress responses, proteostasis, and neuronal survival. Mutations to tRNA cleaving and modifying enzymes have been implicated in several neurodegenerative disorders, and tsRNA levels in the blood are advancing as biomarkers for neurological disease. In this review we summarize the physiological importance of tsRNA in the central nervous system and their relevance to neurological disease.

Comparative Mitogenomic Analysis Reveals Gene and Intron Dynamics in Rubiaceae and Intra-Specific Diversification in Damnacanthus indicus

The dynamic evolution of mitochondrial gene and intron content has been reported across the angiosperms. However, a reference mitochondrial genome (mitogenome) is not available in Rubiaceae. The phylogenetic utility of mitogenome data at a species level is rarely assessed. Here, we assembled mitogenomes of six Damnacanthus indicus (Rubiaceae, Rubioideae) representing two varieties (var. indicus and var. microphyllus). The gene and intron content of D. indicus was compared with mitogenomes from representative angiosperm species and mitochondrial contigs from the other Rubiaceae species. Mitogenome structural rearrangement and sequence divergence in D. indicus were analyzed in six individuals. The size of the mitogenome in D. indicus varied from 417,661 to 419,435 bp. Comparing the number of intact mitochondrial protein-coding genes in other Gentianales taxa (38), D. indicus included 32 genes representing several losses. The intron analysis revealed a shift from cis to trans splicing of a nad1 intron (nad1i728) in D. indicus and it is a shared character with the other four Rubioideae taxa. Two distinct mitogenome structures (type A and B) were identified. Two-step direct repeat-mediated recombination was proposed to explain structural changes between type A and B mitogenomes. The five individuals from two varieties in D. indicus diverged well in the whole mitogenome-level comparison with one exception. Collectively, our study elucidated the mitogenome evolution in Rubiaceae along with D. indicus and showed the reliable phylogenetic utility of the whole mitogenome data at a species-level evolution.

Targeted genome editing of plants and plant cells for biomanufacturing

Plants have provided humans with useful products since antiquity, but in the last 30 years they have also been developed as production platforms for small molecules and recombinant proteins. This initially niche area has blossomed with the growth of the global bioeconomy, and now includes chemical building blocks, polymers and renewable energy. All these applications can be described as “plant molecular farming” (PMF). Despite its potential to increase the sustainability of biologics manufacturing, PMF has yet to be embraced broadly by industry. This reflects a combination of regulatory uncertainty, limited information on process cost structures, and the absence of trained staff and suitable manufacturing capacity. However, the limited adaptation of plants and plant cells to the requirements of industry-scale manufacturing is an equally important hurdle. For example, the targeted genetic manipulation of yeast has been common practice since the 1980s, whereas reliable site-directed mutagenesis in most plants has only become available with the advent of CRISPR/Cas9 and similar genome editing technologies since around 2010. Here we summarize the applications of new genetic engineering technologies to improve plants as biomanufacturing platforms. We start by identifying current bottlenecks in manufacturing, then illustrate the progress that has already been made and discuss the potential for improvement at the molecular, cellular and organism levels. We discuss the effects of metabolic optimization, adaptation of the endomembrane system, modified glycosylation profiles, programmable growth and senescence, protease inactivation, and the expression of enzymes that promote biodegradation. We outline strategies to achieve these modifications by targeted gene modification, considering case-by-case examples of individual improvements and the combined modifications needed to generate a new general-purpose “chassis” for PMF.

Hopeful monsters: unintended sequencing of famously malformed mite mitochondrial tRNAs reveals widespread expression and processing of sense-antisense pairs

Although tRNA structure is one of the most conserved and recognizable shapes in molecular biology, aberrant tRNAs are frequently found in the mitochondrial genomes of metazoans. The extremely degenerate structures of several mitochondrial tRNAs (mt-tRNAs) have led to doubts about their expression and function. Mites from the arachnid superorder Acariformes are predicted to have some of the shortest mt-tRNAs, with a complete loss of cloverleaf-like shape. While performing mitochondrial isolations and recently developed tRNA-seq methods in plant tissue, we inadvertently sequenced the mt-tRNAs from a common plant pest, the acariform mite Tetranychus urticae, to a high enough coverage to detect all previously annotated T. urticae tRNA regions. The results not only confirm expression, CCA-tailing and post-transcriptional base modification of these highly divergent tRNAs, but also revealed paired sense and antisense expression of multiple T. urticae mt-tRNAs. Mirrored expression of mt-tRNA genes has been hypothesized but not previously demonstrated to be common in any system. We discuss the functional roles that these divergent tRNAs could have as both decoding molecules in translation and processing signals in transcript maturation pathways, as well as how sense–antisense pairs add another dimension to the bizarre tRNA biology of mitochondrial genomes.

Non-Redundant tRNA Reference Sequences for Deep Sequencing Analysis of tRNA Abundance and Epitranscriptomic RNA Modifications

Analysis of RNA by deep-sequencing approaches has found widespread application in modern biology. In addition to measurements of RNA abundance under various physiological conditions, such techniques are now widely used for mapping and quantification of RNA modifications. Transfer RNA (tRNA) molecules are among the frequent targets of such investigation, since they contain multiple modified residues. However, the major challenge in tRNA examination is related to a large number of duplicated and point-mutated genes encoding those RNA molecules. Moreover, the existence of multiple isoacceptors/isodecoders complicates both the analysis and read mapping. Existing databases for tRNA sequencing provide near exhaustive listings of tRNA genes, but the use of such highly redundant reference sequences in RNA-seq analyses leads to a large number of ambiguously mapped sequencing reads. Here we describe a relatively simple computational strategy for semi-automatic collapsing of highly redundant tRNA datasets into a non-redundant collection of reference tRNA sequences. The relevance of the approach was validated by analysis of experimentally obtained tRNA-sequencing datasets for different prokaryotic and eukaryotic model organisms. The data demonstrate that non-redundant tRNA reference sequences allow improving unambiguous mapping of deep sequencing data.

The birth of a bacterial tRNA gene by large-scale, tandem duplication events

Organisms differ in the types and numbers of tRNA genes that they carry. While the evolutionary mechanisms behind tRNA gene set evolution have been investigated theoretically and computationally, direct observations of tRNA gene set evolution remain rare. Here, we report the evolution of a tRNA gene set in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The growth defect caused by deleting the single-copy tRNA gene, serCGA, is rapidly compensated by large-scale (45–290 kb) duplications in the chromosome. Each duplication encompasses a second, compensatory tRNA gene (serTGA) and is associated with a rise in tRNA-Ser(UGA) in the mature tRNA pool. We postulate that tRNA-Ser(CGA) elimination increases the translational demand for tRNA-Ser(UGA), a pressure relieved by increasing serTGA copy number. This work demonstrates that tRNA gene sets can evolve through duplication of existing tRNA genes, a phenomenon that may contribute to the presence of multiple, identical tRNA gene copies within genomes.