tRNA-Derived Fragments (tRFs): Emerging New Roles for an Ancient RNA in the Regulation of Gene Expression

Transfer RNA-derived fragments and tRNA halves: biogenesis, biological functions and their roles in diseases

The number of studies on non-coding RNAs has increased substantially in recent years owing to their importance in gene regulation. However, the biological functions of small RNAs from abundant species of housekeeping non-coding RNAs (rRNA, tRNA, etc.) remain a highly studied topic. tRNA-derived small RNAs (tsRNAs) refer to the specific cleavage of tRNAs by specific nucleases [e.g., Dicer and angiogenin (ANG)] in particular cells or tissues or under certain conditions such as stress and hypoxia. tsRNAs are a type of non-coding small RNA that are widely found in the prokaryotic and eukaryotic transcriptomes and are generated from mature tRNAs or precursor tRNAs at different sites. There are two main types of tsRNAs, tRNA-derived fragments (tRFs) and tRNA halves. tRFs are 14-30 nucleotides (nt) long and mainly consist of three subclasses: tRF-5, tRF-3, and tRF-1. tRNA halves, which are 31-40 nt long, are generated by specific cleavage in the anticodon loops of mature tRNAs. There are two types of tRNA halves, 5'-tRNA halves and 3'-tRNA halves. tsRNAs have multiple biological functions including acting as signaling molecules in stress responses and as regulators of gene expression. Additionally, they have been considered to be involved in RNA processing, cell proliferation, translation suppression, the modulation of DNA damage response, and neurodegeneration. More importantly, they are closely related to the occurrence of many human diseases such as tumors, infectious diseases, metabolic diseases, and neurological diseases. Moreover, tsRNAs have the potential to become new biomarkers for disease diagnosis. Continuous investigations will help us to understand their generation and regulatory mechanisms as well as the possible roles of tRFs and tRNA halves.

An integrative piRNA analysis of mouse gametes and zygotes reveals new potential origins and gene regulatory roles

Piwi-interacting RNAs (piRNAs) are a subclass of the small non-coding RNAs (sncRNAs). Their main reported function was to exert control over transposable elements (TEs) in mammalian germline. In this study undertaking a deeper bioinformatics analysis of piRNAs present in mouse oocytes, sperm cells and zygotes, we first elaborated a new piRNA database based on sequences identified as piRNAs by immunoprecipitation with PIWI proteins. Our bioinformatics analysis revealed that, at least in gametes and zygotes, piRNAs could encompass multifunctional cell-dependent regulatory molecules. Indeed, genome analysis of the piRNA mapping density (reads/kb) evidenced in all samples an enrichment of intron-derived piRNAs. Further, piRNA population was classified into sequences not associated to TEs or repeats (NRapiRNAs) and associated to repetitive genome elements (RapiRNAs). In oocytes most of the NRapiRNAs mapped to the 5′UTRs of coding mRNAs, while higher proportion of NRapiRNAs was detected in sperm cells associated to the 3′UTRs of mRNAs. This piRNA complementarity to mRNA UTRs suggests key post-transcriptional regulatory roles over mRNAs such as those encoding MHC genes. In addition, a striking association of RapiRNA with long non-coding RNAs (lncRNAs) was identified. piRNAs associated with relevant lncRNAs such as: Rab26os and GAS5 and key mRNAs, were particularly assessed.

Golgi Reassembly and Stacking Protein (GRASP) Participates in Vesicle-Mediated RNA Export in Cryptococcus Neoformans

Golgi reassembly and stacking protein (GRASP) is required for polysaccharide secretion and virulence in Cryptococcus neoformans. In fungal species, extracellular vesicles (EVs) participate in the export of polysaccharides, proteins and RNA. In the present work, we investigated if EV-mediated RNA export is functionally connected with GRASP in C. neoformans using a graspΔ mutant. Since GRASP-mediated unconventional secretion involves autophagosome formation in yeast, we included the atg7Δ mutant with defective autophagic mechanisms in our analysis. All fungal strains exported EVs but deletion of GRASP or ATG7 profoundly affected vesicular dimensions. The mRNA content of the graspΔ EVs differed substantially from that of the other two strains. The transcripts associated to the endoplasmic reticulum were highly abundant transcripts in graspΔ EVs. Among non-coding RNAs (ncRNAs), tRNA fragments were the most abundant in both mutant EVs but graspΔ EVs alone concentrated 22 exclusive sequences. In general, our results showed that the EV RNA content from atg7Δ and WT were more related than the RNA content of graspΔ, suggesting that GRASP, but not the autophagy regulator Atg7, is involved in the EV export of RNA. This is a previously unknown function for a key regulator of unconventional secretion in eukaryotic cells.

tRNA fragments (tRFs) guide Ago to regulate gene expression post-transcriptionally in a Dicer-independent manner

tRNA related RNA fragments (tRFs), also known as tRNA-derived RNAs (tdRNAs), are abundant small RNAs reported to be associated with Argonaute proteins, yet their function is unclear. We show that endogenous 18 nucleotide tRFs derived from the 3' ends of tRNAs (tRF-3) post-transcriptionally repress genes in HEK293T cells in culture. tRF-3 levels increase upon parental tRNA overexpression. This represses target genes with a sequence complementary to the tRF-3 in the 3' UTR. The tRF-3-mediated repression is Dicer-independent, Argonaute-dependent, and the targets are recognized by sequence complementarity. Furthermore, tRF-3:target mRNA pairs in the RNA induced silencing complex associate with GW182 proteins, known to repress translation and promote the degradation of target mRNAs. RNA-seq demonstrates that endogenous target genes are specifically decreased upon tRF-3 induction. Therefore, Dicer-independent tRF-3s, generated upon tRNA overexpression, repress genes post-transcriptionally through an Argonaute-GW182 containing RISC via sequence matches with target mRNAs.

A Novel Class of tRNA-Derived Small Non-Coding RNAs Respond to Myocardial Hypertrophy and Contribute to Intergenerational Inheritance

tRNA-derived fragments (tRFs) are a new class of non-coding RNA that play an important role in regulating cellular RNA processing and protein translation. However, there is currently no study reporting the influence of tRFs on myocardial hypertrophy. In this study, we used an isoproterenol (ISO)-induced myocardial hypertrophy rat model. Small RNA (<40 nts) transcriptome sequencing was used to select differentially expressed tRFs. We also compared the tRFs expression pattern in F0 sperm and the hearts of F1 offspring between the myocardial hypertrophy group (Hyp) and the control group (Con). Isoproterenol successfully induced a typical cardiac hypertrophy model in our study. Small RNA-seq revealed that tRFs were extremely enriched (84%) in the Hyp heart. Overexpression of tRFs1 and tRFs2 both enlarged the surface area of cardiac cells and increased expression of hypertrophic markers (ANF, BNP, and β-MHC). Luciferase reporter assay identified that tRFs1 directly target 3′UTR of Timp3. tRFs1, tRFs2, tRFs3, and tRFs4 were also highly expressed in Hyp F0 sperm and in Hyp F1 offspring hearts, but there was no differential expression of tRFs7, tRFs9, and tRFs10. Compared to Con F1 offspring, Hyp F1 offspring had elevated expression levels of β-MHC and ANP genes, and they had increased fibrosis and apoptosis in their hearts. These results demonstrated that tRFs are involved in regulating the response of myocardial hypertrophy. Besides, tRFs might serve as novel epigenetic factors that contribute to the intergenerational inheritance of cardiac hypertrophy.

Targeting RNA in mammalian systems with small molecules

The recognition of RNA functions beyond canonical protein synthesis has challenged the central dogma of molecular biology. Indeed, RNA is now known to directly regulate many important cellular processes, including transcription, splicing, translation, and epigenetic modifications. The misregulation of these processes in disease has led to an appreciation of RNA as a therapeutic target. This potential was first recognized in bacteria and viruses, but discoveries of new RNA classes following the sequencing of the human genome have invigorated exploration of its disease-related functions in mammals. As stable structure formation is evolving as a hallmark of mammalian RNAs, the prospect of utilizing small molecules to specifically probe the function of RNA structural domains and their interactions is gaining increased recognition. To date, researchers have discovered bioactive small molecules that modulate phenotypes by binding to expanded repeats, microRNAs, G-quadruplex structures, and RNA splice sites in neurological disorders, cancers, and other diseases. The lessons learned from achieving these successes both call for additional studies and encourage exploration of the plethora of mammalian RNAs whose precise mechanisms of action remain to be elucidated. Efforts toward understanding fundamental principles of small molecule-RNA recognition combined with advances in methodology development should pave the way toward targeting emerging RNA classes such as long noncoding RNAs. Together, these endeavors can unlock the full potential of small molecule-based probing of RNA-regulated processes and enable us to discover new biology and underexplored avenues for therapeutic intervention in human disease. This article is categorized under: RNA Methods > RNA Analyses In Vitro and In Silico RNA Interactions with Proteins and Other Molecules > Small Molecule-RNA Interactions RNA in Disease and Development > RNA in Disease.

Effects of Acute Aerobic Exercise on Rats Serum Extracellular Vesicles Diameter, Concentration and Small RNAs Content

Physical exercise stimulates organs, mainly the skeletal muscle, to release a broad range of molecules, recently dubbed exerkines. Among them, RNAs, such as miRNAs, piRNAs, and tRNAs loaded in extracellular vesicles (EVs) have the potential to play a significant role in the way muscle and other organs communicate to translate exercise into health. Low, moderate and high intensity treadmill protocols were applied to rat groups, aiming to investigate the impact of exercise on serum EVs and their associated small RNA molecules. Transmission electron microscopy, resistive pulse sensing, and western blotting were used to investigate EVs morphology, size distribution, concentration and EVs marker proteins. Small RNA libraries from EVs RNA were sequenced. Exercise did not change EVs size, while increased EVs concentration. Twelve miRNAs were found differentially expressed after exercise: rno-miR-128-3p, 103-3p, 330-5p, 148a-3p, 191a-5p, 10b-5p, 93-5p, 25-3p, 142-5p, 3068-3p, 142-3p, and 410-3p. No piRNA was found differentially expressed, and one tRNA, trna8336, was found down-regulated after exercise. The differentially expressed miRNAs were predicted to target genes involved in the MAPK pathway. A single bout of exercise impacts EVs and their small RNA load, reinforcing the need for a more detailed investigation into EVs and their load as mediators of health-promoting exercise.

The dysregulation of tRNAs and tRNA derivatives in cancer

Transfer RNAs (tRNAs), traditionally considered to participate in protein translation, were interspersed in the entire genome. Recent studies suggested that dysregulation was observed in not only tRNAs, but also tRNA derivatives generated by the specific cleavage of pre- and mature tRNAs in the progression of cancer. Accumulating evidence had identified that certain tRNAs and tRNA derivatives were involved in proliferation, metastasis and invasiveness of cancer cell, as well as tumor growth and angiogenesis in several malignant human tumors. This paper reviews the importance of the dysregulation of tRNAs and tRNA derivatives during the development of cancer, such as breast cancer, lung cancer, and melanoma, aiming at a better understanding of the tumorigenesis and providing new ideas for the treatment of these cancers.

RNA Cloaking by Reversible Acylation

We describe a selective and mild chemical approach for controlling RNA hybridization, folding, and enzyme interactions. Reaction of RNAs in aqueous buffer with an azide-substituted acylating agent (100-200 mm) yields several 2'-OH acylations per RNA strand in as little as 10 min. This poly-acylated ("cloaked") RNA is strongly blocked from hybridization with complementary nucleic acids, from cleavage by RNA-processing enzymes, and from folding into active aptamer structures. Importantly, treatment with a water-soluble phosphine triggers a Staudinger reduction of the azide groups, resulting in spontaneous loss of acyl groups ("uncloaking"). This fully restores RNA folding and biochemical activity.

Heavy Chronic Intermittent Ethanol Exposure Alters Small Noncoding RNAs in Mouse Sperm and Epididymosomes

While the risks of maternal alcohol abuse during pregnancy are well-established, several preclinical studies suggest that chronic preconception alcohol consumption by either parent may also have significance consequences for offspring health and development. Notably, since isogenic male mice used in these studies are not involved in gestation or rearing of offspring, the cross-generational effects of paternal alcohol exposure suggest a germline-based epigenetic mechanism. Many recent studies have demonstrated that the effects of paternal environmental exposures such as stress or malnutrition can be transmitted to the next generation via alterations to small noncoding RNAs in sperm. Therefore, we used high throughput sequencing to examine the effect of preconception ethanol on small noncoding RNAs in sperm. We found that chronic intermittent ethanol exposure altered several small noncoding RNAs from three of the major small RNA classes in sperm, tRNA-derived small RNA (tDR), mitochondrial small RNA, and microRNA. Six of the ethanol-responsive small noncoding RNAs were evaluated with RT-qPCR on a separate cohort of mice and five of the six were confirmed to be altered by chronic ethanol exposure, supporting the validity of the sequencing results. In addition to altered sperm RNA abundance, chronic ethanol exposure affected post-transcriptional modifications to sperm small noncoding RNAs, increasing two nucleoside modifications previously identified in mitochondrial tRNA. Furthermore, we found that chronic ethanol reduced epididymal expression of a tRNA methyltransferase, Nsun2, known to directly regulate tDR biogenesis. Finally, ethanol-responsive sperm tDR are similarly altered in extracellular vesicles of the epididymis (i.e., epididymosomes), supporting the hypothesis that alterations to sperm tDR emerge in the epididymis and that epididymosomes are the primary source of small noncoding RNAs in sperm. These results add chronic ethanol to the growing list of paternal exposures that can affect small noncoding RNA abundance and nucleoside modifications in sperm. As small noncoding RNAs in sperm have been shown to causally induce heritable phenotypes in offspring, additional research is warranted to understand the potential effects of ethanol-responsive sperm small noncoding RNAs on offspring health and development.

Multifaceted functional implications of an endogenously expressed tRNA fragment in the vector mosquito Aedes aegypti

The mosquito Aedes aegypti is the primary vector of human arboviral diseases caused by dengue, chikungunya and Zika viruses. Many studies have shown the potential roles of small RNA molecules such as microRNA, small interfering RNA and PIWI-interacting RNA in vector mosquitoes. The function of tRNA fragments (tRF), the newly discovered class of small RNAs, in mosquitoes is not known. In this study, we show that specific tRFs are expressed in significantly differential manner between males and females of Ae. aegypti strains. Specific tRFs also show differential response during developmental transition from larvae to adults, as well as after blood feeding of adult females. The expression pattern of tRFs upon blood feeding varied depending upon if the blood contained dengue virus, and also if the females were treated with antibiotic prior to feeding to cleanse of the gut bacteria. Our findings show that a single tRF derived from the precursor sequences of a tRNA-Gly was differentially expressed between males and females, developmental transitions and also upon blood feeding by females of two laboratory strains that vary in midgut susceptibility to dengue virus infection. The multifaceted functional implications of this specific tRF suggest that biogenesis of small regulatory molecules from a tRNA can have wide ranging effects on key aspects of Ae. aegypti vector biology.

The mosquito Aedes aegypti is a major vector of arboviral diseases in subtropics and tropics. The confounding effects of immature development and adult microbiome on the ability of Ae. aegypti to transmit diseases (vector competence) have gained renewed attention in the recent years. However, the molecular nature of these links/ effects remains unknown. This is major gap in knowledge regarding how vector competence is regulated at molecular level, and how that regulation may be variable among different strains of this mosquito. In this study, we investigated expression of newly discovered class of small RNAs, called tRNA fragments (tRF) in Ae. aegypti strains. Based on small RNA sequencing and bioinformatics analyses, we show that tRFs are expressed in Ae. aegypti, and they are associated with significant changes in expression between males and females, during development stages, and post blood feeding responses. A single tRF showed association with sex-biased expression, developmental regulation and in response to blood meals between Moyo-S and Moyo-R strains that differ in midgut susceptibility to dengue virus. The findings of this study are expected to guide future research efforts directed toward examining detailed regulatory mechanisms of tRFs in vector competence of Ae. aegypti to disease transmission.

Signaling to and from the RNA Polymerase III Transcription and Processing Machinery

RNA polymerase (Pol) III has a specialized role in transcribing the most abundant RNAs in eukaryotic cells, transfer RNAs (tRNAs), along with other ubiquitous small noncoding RNAs, many of which have functions related to the ribosome and protein synthesis. The high energetic cost of producing these RNAs and their central role in protein synthesis underlie the robust regulation of Pol III transcription in response to nutrients and stress by growth regulatory pathways. Downstream of Pol III, signaling impacts posttranscriptional processes affecting tRNA function in translation and tRNA cleavage into smaller fragments that are increasingly attributed with novel cellular activities. In this review, we consider how nutrients and stress control Pol III transcription via its factors and its negative regulator, Maf1. We highlight recent work showing that the composition of the tRNA population and the function of individual tRNAs is dynamically controlled and that unrestrained Pol III transcription can reprogram central metabolic pathways.

MINTbase v2.0: a comprehensive database for tRNA-derived fragments that includes nuclear and mitochondrial fragments from all The Cancer Genome Atlas projects

MINTbase is a repository that comprises nuclear and mitochondrial tRNA-derived fragments ('tRFs') found in multiple human tissues. The original version of MINTbase comprised tRFs obtained from 768 transcriptomic datasets. We used our deterministic and exhaustive tRF mining pipeline to process all of The Cancer Genome Atlas datasets (TCGA). We identified 23 413 tRFs with abundance of ≥ 1.0 reads-per-million (RPM). To facilitate further studies of tRFs by the community, we just released version 2.0 of MINTbase that contains information about 26 531 distinct human tRFs from 11 719 human datasets as of October 2017. Key new elements include: the ability to filter tRFs on-the-fly by minimum abundance thresholding; the ability to filter tRFs by tissue keywords; easy access to information about a tRF's maximum abundance and the datasets that contain it; the ability to generate relative abundance plots for tRFs across cancer types and convert them into embeddable figures; MODOMICS information about modifications of the parental tRNA, etc. Version 2.0 of MINTbase contains 15x more datasets and nearly 4x more distinct tRFs than the original version, yet continues to offer fast, interactive access to its contents. Version 2.0 is available freely at http://cm.jefferson.edu/MINTbase/.

Differentially Expressed tRNA-Derived Small RNAs Co-Sediment Primarily with Non-Polysomal Fractions in Drosophila

Recent studies point to the existence of poorly characterized small regulatory RNAs generated from mRNAs, rRNAs and tRNAs. To explore the subcellular location of tRNA-derived small RNAs, 0–1 and 7–8 h Drosophila embryos were fractionated on sucrose density gradients. Analysis of 12,553,921 deep-sequencing reads from unfractionated and fractionated Drosophila embryos has revealed that tRFs, which are detected mainly from the 5’ends of tRNAs, co-sediment with the non-polysomal fractions. Interestingly, the expression levels of a subset of tRFs change temporally following the maternal-to-zygotic transition in embryos. We detected non-polysomal association of tRFs in S2 cells as well. Differential tRF expression pattern points to developmental significance at the organismal level. These results suggest that tRFs are associated primarily with the non-polysomal complexes in Drosophila embryos and S2 cells.

Protein folding and tRNA biology

Polypeptides can fold into tertiary structures while they are synthesized by the ribosome. In addition to the amino acid sequence, protein folding is determined by several factors within the cell. Among others, the folding pathway of a nascent polypeptide can be affected by transient interactions with other proteins, ligands, or the ribosome, as well as by the translocation through membrane pores. Particularly, the translation machinery and the population of tRNA under different physiological or adaptive responses can dramatically affect protein folding. This review summarizes the scientific evidence describing the role of translation kinetics and tRNA populations on protein folding and addresses current efforts to better understand tRNA biology. It is organized into three main parts, which are focused on: (i) protein folding in the cellular context; (ii) tRNA biology and the complexity of the tRNA population; and (iii) available methods and technical challenges in the characterization of tRNA pools. In this manner, this work illustrates the ways by which functional properties of proteins may be modulated by cellular tRNA populations.

The emerging complexity of the tRNA world: mammalian tRNAs beyond protein synthesis

The discovery of the genetic code and tRNAs as decoders of the code transformed life science. However, after establishing the role of tRNAs in protein synthesis, the field moved to other parts of the RNA world. Now, tRNA research is blooming again, with demonstration of the involvement of tRNAs in various other pathways beyond translation and in adapting translation to environmental cues. These roles are linked to the presence of tRNA sequence variants known as isoacceptors and isodecoders, various tRNA base modifications, the versatility of protein binding partners and tRNA fragmentation events, all of which collectively create an incalculable complexity. This complexity provides a vast repertoire of tRNA species that can serve various functions in cellular homeostasis and in adaptation of cellular functions to changing environments, and it likely arose from the fundamental role of RNAs in early evolution.

unitas: the universal tool for annotation of small RNAs

BACKGROUND

Next generation sequencing is a key technique in small RNA biology research that has led to the discovery of functionally different classes of small non-coding RNAs in the past years. However, reliable annotation of the extensive amounts of small non-coding RNA data produced by high-throughput sequencing is time-consuming and requires robust bioinformatics expertise. Moreover, existing tools have a number of shortcomings including a lack of sensitivity under certain conditions, limited number of supported species or detectable sub-classes of small RNAs.

RESULTS

Here we introduce unitas, an out-of-the-box ready software for complete annotation of small RNA sequence datasets, supporting the wide range of species for which non-coding RNA reference sequences are available in the Ensembl databases (currently more than 800). unitas combines high quality annotation and numerous analysis features in a user-friendly manner. A complete annotation can be started with one simple shell command, making unitas particularly useful for researchers not having access to a bioinformatics facility. Noteworthy, the algorithms implemented in unitas are on par or even outperform comparable existing tools for small RNA annotation that map to publicly available ncRNA databases.

CONCLUSIONS

unitas brings together annotation and analysis features that hitherto required the installation of numerous different bioinformatics tools which can pose a challenge for the non-expert user. With this, unitas overcomes the problem of read normalization. Moreover, the high quality of sequence annotation and analysis, paired with the ease of use, make unitas a valuable tool for researchers in all fields connected to small RNA biology.

The nuclear and organellar tRNA-derived RNA fragment population in Arabidopsis thaliana is highly dynamic

In the expanding repertoire of small noncoding RNAs (ncRNAs), tRNA-derived RNA fragments (tRFs) have been identified in all domains of life. Their existence in plants has been already proven but no detailed analysis has been performed. Here, short tRFs of 19–26 nucleotides were retrieved from Arabidopsis thaliana small RNA libraries obtained from various tissues, plants submitted to abiotic stress or fractions immunoprecipitated with ARGONAUTE 1 (AGO1). Large differences in the tRF populations of each extract were observed. Depending on the tRNA, either tRF-5D (due to a cleavage in the D region) or tRF-3T (via a cleavage in the T region) were found and hot spots of tRNA cleavages have been identified. Interestingly, up to 25% of the tRFs originate from plastid tRNAs and we provide evidence that mitochondrial tRNAs can also be a source of tRFs. Very specific tRF-5D deriving not only from nucleus-encoded but also from plastid-encoded tRNAs are strongly enriched in AGO1 immunoprecipitates. We demonstrate that the organellar tRFs are not found within chloroplasts or mitochondria but rather accumulate outside the organelles. These observations suggest that some organellar tRFs could play regulatory functions within the plant cell and may be part of a signaling pathway.

Transfer RNAs with novel cloverleaf structures

We report the identification of novel tRNA species with 12-base pair amino-acid acceptor branches composed of longer acceptor stem and shorter T-stem. While canonical tRNAs have a 7/5 configuration of the branch, the novel tRNAs have either 8/4 or 9/3 structure. They were found during the search for selenocysteine tRNAs in terabytes of genome, metagenome and metatranscriptome sequences. Certain bacteria and their phages employ the 8/4 structure for serine and histidine tRNAs, while minor cysteine and selenocysteine tRNA species may have a modified 8/4 structure with one bulge nucleotide. In Acidobacteria, tRNAs with 8/4 and 9/3 structures may function as missense and nonsense suppressor tRNAs and/or regulatory noncoding RNAs. In δ-proteobacteria, an additional cysteine tRNA with an 8/4 structure mimics selenocysteine tRNA and may function as opal suppressor. We examined the potential translation function of suppressor tRNA species in Escherichia coli; tRNAs with 8/4 or 9/3 structures efficiently inserted serine, alanine and cysteine in response to stop and sense codons, depending on the identity element and anticodon sequence of the tRNA. These findings expand our view of how tRNA, and possibly the genetic code, is diversified in nature.

LTR-Retrotransposon Control by tRNA-Derived Small RNAs

Transposon reactivation is an inherent danger in cells that lose epigenetic silencing during developmental reprogramming. In the mouse, long terminal repeat (LTR)-retrotransposons, or endogenous retroviruses (ERV), account for most novel insertions and are expressed in the absence of histone H3 lysine 9 trimethylation in preimplantation stem cells. We found abundant 18 nt tRNA-derived small RNA (tRF) in these cells and ubiquitously expressed 22 nt tRFs that include the 3' terminal CCA of mature tRNAs and target the tRNA primer binding site (PBS) essential for ERV reverse transcription. We show that the two most active ERV families, IAP and MusD/ETn, are major targets and are strongly inhibited by tRFs in retrotransposition assays. 22 nt tRFs post-transcriptionally silence coding-competent ERVs, while 18 nt tRFs specifically interfere with reverse transcription and retrotransposon mobility. The PBS offers a unique target to specifically inhibit LTR-retrotransposons, and tRF-targeting is a potentially highly conserved mechanism of small RNA-mediated transposon control.

SMCHD1 regulates a limited set of gene clusters on autosomal chromosomes

BACKGROUND

Facioscapulohumeral muscular dystrophy (FSHD) is in most cases caused by a contraction of the D4Z4 macrosatellite repeat on chromosome 4 (FSHD1) or by mutations in the SMCHD1 or DNMT3B gene (FSHD2). Both situations result in the incomplete epigenetic repression of the D4Z4-encoded retrogene DUX4 in somatic cells, leading to the aberrant expression of DUX4 in the skeletal muscle. In mice, Smchd1 regulates chromatin repression at different loci, having a role in CpG methylation establishment and/or maintenance.

METHODS

To investigate the global effects of harboring heterozygous SMCHD1 mutations on DNA methylation in humans, we combined 450k methylation analysis on mononuclear monocytes from female heterozygous SMCHD1 mutation carriers and unaffected controls with reduced representation bisulfite sequencing (RRBS) on FSHD2 and control myoblast cell lines. Candidate loci were then evaluated for SMCHD1 binding using ChIP-qPCR and expression was evaluated using RT-qPCR.

RESULTS

We identified a limited number of clustered autosomal loci with CpG hypomethylation in SMCHD1 mutation carriers: the protocadherin (PCDH) cluster on chromosome 5, the transfer RNA (tRNA) and 5S rRNA clusters on chromosome 1, the HOXB and HOXD clusters on chromosomes 17 and 2, respectively, and the D4Z4 repeats on chromosomes 4 and 10. Furthermore, minor increases in RNA expression were seen in FSHD2 myoblasts for some of the PCDHβ cluster isoforms, tRNA isoforms, and a HOXB isoform in comparison to controls, in addition to the previously reported effects on DUX4 expression. SMCHD1 was bound at DNAseI hypersensitivity sites known to regulate the PCDHβ cluster and at the chromosome 1 tRNA cluster, with decreased binding in SMCHD1 mutation carriers at the PCDHβ cluster sites.

CONCLUSIONS

Our study is the first to investigate the global methylation effects in humans resulting from heterozygous mutations in SMCHD1. Our results suggest that SMCHD1 acts as a repressor on a limited set of autosomal gene clusters, as an observed reduction in methylation associates with a loss of SMCHD1 binding and increased expression for some of the loci.

Exercise alters mouse sperm small noncoding RNAs and induces a transgenerational modification of male offspring conditioned fear and anxiety

There is growing evidence that the preconceptual lifestyle and other environmental exposures of a father can significantly alter the physiological and behavioral phenotypes of their children. We and others have shown that paternal preconception stress, regardless of whether the stress was experienced during early-life or adulthood, results in offspring with altered anxiety and depression-related behaviors, attributed to hypothalamic-pituitary-adrenal axis dysregulation. The transgenerational response to paternal preconceptual stress is believed to be mediated by sperm-borne small noncoding RNAs, specifically microRNAs. As physical activity confers physical and mental health benefits for the individual, we used a model of voluntary wheel-running and investigated the transgenerational response to paternal exercise. We found that male offspring of runners had suppressed reinstatement of juvenile fear memory, and reduced anxiety in the light-dark apparatus during adulthood. No changes in these affective behaviors were observed in female offspring. We were surprised to find that running had a limited impact on sperm-borne microRNAs. The levels of three unique microRNAs (miR-19b, miR-455 and miR-133a) were found to be altered in the sperm of runners. In addition, we discovered that the levels of two species of tRNA-derived RNAs (tDRs)-tRNA-Gly and tRNA-Pro-were also altered by running. Taken together, we believe this is the first evidence that paternal exercise is associated with an anxiolytic behavioral phenotype of male offspring and altered levels of small noncoding RNAs in sperm. These small noncoding RNAs are known to have an impact on post-transcriptional gene regulation and can thus change the developmental trajectory of offspring brains and associated affective behaviors.

Dissecting miRNA gene repression on single cell level with an advanced fluorescent reporter system

Despite major advances on miRNA profiling and target predictions, functional readouts for endogenous miRNAs are limited and frequently lead to contradicting conclusions. Numerous approaches including functional high-throughput and miRISC complex evaluations suggest that the functional miRNAome differs from the predictions based on quantitative sRNA profiling. To resolve the apparent contradiction of expression versus function, we generated and applied a fluorescence reporter gene assay enabling single cell analysis. This approach integrates and adapts a mathematical model for miRNA-driven gene repression. This model predicts three distinct miRNA-groups with unique repression activities (low, mid and high) governed not just by expression levels but also by miRNA/target-binding capability. Here, we demonstrate the feasibility of the system by applying controlled concentrations of synthetic siRNAs and in parallel, altering target-binding capability on corresponding reporter-constructs. Furthermore, we compared miRNA-profiles with the modeled predictions of 29 individual candidates. We demonstrate that expression levels only partially reflect the miRNA function, fitting to the model-projected groups of different activities. Furthermore, we demonstrate that subcellular localization of miRNAs impacts functionality. Our results imply that miRNA profiling alone cannot define their repression activity. The gene regulatory function is a dynamic and complex process beyond a minimalistic conception of "highly expressed equals high repression".

MINTmap: fast and exhaustive profiling of nuclear and mitochondrial tRNA fragments from short RNA-seq data

Transfer RNA fragments (tRFs) are an established class of constitutive regulatory molecules that arise from precursor and mature tRNAs. RNA deep sequencing (RNA-seq) has greatly facilitated the study of tRFs. However, the repeat nature of the tRNA templates and the idiosyncrasies of tRNA sequences necessitate the development and use of methodologies that differ markedly from those used to analyze RNA-seq data when studying microRNAs (miRNAs) or messenger RNAs (mRNAs). Here we present MINTmap (for MItochondrial and Nuclear TRF mapping), a method and a software package that was developed specifically for the quick, deterministic and exhaustive identification of tRFs in short RNA-seq datasets. In addition to identifying them, MINTmap is able to unambiguously calculate and report both raw and normalized abundances for the discovered tRFs. Furthermore, to ensure specificity, MINTmap identifies the subset of discovered tRFs that could be originating outside of tRNA space and flags them as candidate false positives. Our comparative analysis shows that MINTmap exhibits superior sensitivity and specificity to other available methods while also being exceptionally fast. The MINTmap codes are available through https://github.com/TJU-CMC-Org/MINTmap/ under an open source GNU GPL v3.0 license.

Secretory pathway optimization of CHO producer cells by co-engineering of the mitosRNA-1978 target genes CerS2 and Tbc1D20

Chinese Hamster Ovary (CHO) cells are the most commonly used host for the production of biopharmaceuticals. Although transcription and translation engineering strategies have been employed to generate high-producer cell clones, the secretory pathway still remains a bottleneck in cellular productivity. In this study we show that ectopic expression of a human mitochondrial genome-encoded small RNA (mitosRNA-1978) in an IgG expressing CHO cell line strongly improved specific productivity by functioning in a microRNA-like fashion. By next generation sequencing we identified two endoplasmic reticulum (ER)-localized proteins, Ceramide Synthase 2 (CerS2) and the Rab1 GAP Tbc domain family member 20 (Tbc1D20), as target genes of mitosRNA-1978. Combined transient siRNA-mediated knockdown of CerS2 and Tbc1D20 resulted in increased specific productivity of CHO-IgG cells, thus recapitulating the mitosRNA-1978 phenotype. In support of a function in vesicular trafficking at the level of the ER, we provide evidence for altered cellular ceramide composition upon CerS2 knockdown and increased activity of Rab1 in CHO-IgG cells depleted of Tbc1D20. Importantly, in a fed-batch process, the combined stable knockdown of CerS2 and Tbc1D20 in CHO-IgG cells resulted in dramatically increased antibody production which was accompanied by enhanced cell growth. Thus, by identifying mitosRNA-1978 target genes in combination with an informed shRNA-mediated co-engineering approach we successfully optimized the secretory capacity of CHO producer cells used for the manufacturing of therapeutic proteins.

tRNA-Derived RNA Fragments Associate with Human Multisynthetase Complex (MSC) and Modulate Ribosomal Protein Translation

The functionality of small RNAs from abundant species of "housekeeping" noncoding RNAs (e.g., rRNA, tRNA, snRNA, snoRNA, etc.) remains a highly studied topic. The current state of research on short RNAs derived from transfer RNA (tRNA), called tRNA-derived fragments (tRFs), has been restricted largely to expression studies and limited functional studies. 5' tRFs are known translational inhibitors in mammalian cells, yet little is known about their functionality. Here we report on the first experimental evidence of the tRF protein interactome, identifying the mammalian multisynthetase complex as the primary interactor of the 5' tRF Gln19. We also present proteome-wide SILAC evidence that 5' tRFs increase ribosomal and poly(A)-binding protein translation.

YM500v3: a database for small RNA sequencing in human cancer research

We previously presented the YM500 database, which contains >8000 small RNA sequencing (smRNA-seq) data sets and integrated analysis results for various cancer miRNome studies. In the updated YM500v3 database (http://ngs.ym.edu.tw/ym500/) presented herein, we not only focus on miRNAs but also on other functional small non-coding RNAs (sncRNAs), such as PIWI-interacting RNAs (piRNAs), tRNA-derived fragments (tRFs), small nuclear RNAs (snRNAs) and small nucleolar RNAs (snoRNAs). There is growing knowledge of the role of sncRNAs in gene regulation and tumorigenesis. We have also incorporated >10 000 cancer-related RNA-seq and >3000 more smRNA-seq data sets into the YM500v3 database. Furthermore, there are two main new sections, ‘Survival' and ‘Cancer', in this updated version. The ‘Survival’ section provides the survival analysis results in all cancer types or in a user-defined group of samples for a specific sncRNA. The ‘Cancer’ section provides the results of differential expression analyses, miRNA–gene interactions and cancer miRNA-related pathways. In the ‘Expression’ section, sncRNA expression profiles across cancer and sample types are newly provided. Cancer-related sncRNAs hold potential for both biotech applications and basic research.

A tRNA-derived fragment competes with mRNA for ribosome binding and regulates translation during stress

Posttranscriptional processing of RNA molecules is a common strategy to enlarge the structural and functional repertoire of RNomes observed in all 3 domains of life. Fragmentation of RNA molecules of basically all functional classes has been reported to yield smaller non-protein coding RNAs (ncRNAs) that typically possess different roles compared with their parental transcripts. Here we show that a valine tRNA-derived fragment (Val-tRF) that is produced under certain stress conditions in the halophilic archaeon Haloferax volcanii is capable of binding to the small ribosomal subunit. As a consequence of Val-tRF binding mRNA is displaced from the initiation complex which results in global translation attenuation in vivo and in vitro. The fact that the archaeal Val-tRF also inhibits eukaryal as well as bacterial protein biosynthesis implies a functionally conserved mode of action. While tRFs and tRNA halves have been amply identified in recent RNA-seq project, Val-tRF described herein represents one of the first functionally characterized tRNA processing products to date.

Comprehensive Small RNA-Seq of Adeno-Associated Virus (AAV)-Infected Human Cells Detects Patterns of Novel, Non-Coding AAV RNAs in the Absence of Cellular miRNA Regulation

Most DNA viruses express small regulatory RNAs, which interfere with viral or cellular gene expression. For adeno-associated virus (AAV), a small ssDNA virus with a complex biphasic life cycle miRNAs or other small regulatory RNAs have not yet been described. This is the first comprehensive Illumina-based RNA-Seq analysis of small RNAs expressed by AAV alone or upon co-infection with helper adenovirus or HSV. Several hotspots of AAV-specific small RNAs were detected mostly close to or within the AAV-ITR and apparently transcribed from the newly identified anti-p5 promoter. An additional small RNA hotspot was located downstream of the p40 promoter, from where transcription of non-coding RNAs associated with the inhibition of adenovirus replication were recently described. Parallel detection of known Ad and HSV miRNAs indirectly validated the newly identified small AAV RNA species. The predominant small RNAs were analyzed on Northern blots and by human argonaute protein-mediated co-immunoprecipitation. None of the small AAV RNAs showed characteristics of bona fide miRNAs, but characteristics of alternative RNA processing indicative of differentially regulated AAV promoter-associated small RNAs. Furthermore, the AAV-induced regulation of cellular miRNA levels was analyzed at different time points post infection. In contrast to other virus groups AAV infection had virtually no effect on the expression of cellular miRNA, which underscores the long-established concept that wild-type AAV infection is apathogenic.

Non-canonical roles of tRNAs and tRNA mimics in bacterial cell biology

Transfer RNAs (tRNAs) are the macromolecules that transfer activated amino acids from aminoacyl-tRNA synthetases to the ribosome, where they are used for the mRNA guided synthesis of proteins. Transfer RNAs are ancient molecules, perhaps even predating the existence of the translation machinery. Albeit old, these molecules are tremendously conserved, a characteristic that is well illustrated by the fact that some bacterial tRNAs are efficient and specific substrates of eukaryotic aminoacyl-tRNA synthetases and ribosomes. Considering their ancient origin and high structural conservation, it is not surprising that tRNAs have been hijacked during evolution for functions outside of translation. These roles beyond translation include synthetic, regulatory and information functions within the cell. Here we provide an overview of the non-canonical roles of tRNAs and their mimics in bacteria, and discuss some of the common themes that arise when comparing these different functions.

tRF2Cancer: A web server to detect tRNA-derived small RNA fragments (tRFs) and their expression in multiple cancers

tRNA-derived small RNA fragments (tRFs) are one class of small non-coding RNAs derived from transfer RNAs (tRNAs). tRFs play important roles in cellular processes and are involved in multiple cancers. High-throughput small RNA (sRNA) sequencing experiments can detect all the cellular expressed sRNAs, including tRFs. However, distinguishing genuine tRFs from RNA fragments generated by random degradation remains a major challenge. In this study, we developed an integrated web-based computing system, tRF2Cancer, to accurately identify tRFs from sRNA deep-sequencing data and evaluate their expression in multiple cancers. The binomial test was introduced to evaluate whether reads from a small RNA-seq data set represent tRFs or degraded fragments. A classification method was then used to annotate the types of tRFs based on their sites of origin in pre-tRNA or mature tRNA. We applied the pipeline to analyze 10 991 data sets from 32 types of cancers and identified thousands of expressed tRFs. A tool called ‘tRFinCancer’ was developed to facilitate the users to inspect the expression of tRFs across different types of cancers. Another tool called ‘tRFBrowser’ shows both the sites of origin and the distribution of chemical modification sites in tRFs on their source tRNA. The tRF2Cancer web server is available at http://rna.sysu.edu.cn/tRFfinder/.