Transfer RNA detection by small RNA deep sequencing and disease association with myelodysplastic syndromes

Disruption of mitochondrial energy metabolism is a putative pathogenesis of Diamond-Blackfan anemia

Energy metabolism in the context of erythropoiesis and related diseases remains largely unexplored. Here, we developed a primary cell model by differentiating hematopoietic stem progenitor cells toward the erythroid lineage and suppressing the mitochondrial oxidative phosphorylation (OXPHOS) pathway. OXPHOS suppression led to differentiation failure of erythroid progenitors and defects in ribosome biogenesis. Ran GTPase-activating protein 1 (RanGAP1) was identified as a target of mitochondrial OXPHOS for ribosomal defects during erythropoiesis. Overexpression of RanGAP1 largely alleviated erythroid defects resulting from OXPHOS suppression. Coenzyme Q10, an activator of OXPHOS, largely rescued erythroid defects and increased RanGAP1 expression. Patients with Diamond-Blackfan anemia (DBA) exhibited OXPHOS suppression and a concomitant suppression of ribosome biogenesis. RNA-seq analysis implied that the substantial mutation (approximately 10%) in OXPHOS genes accounts for OXPHOS suppression in these patients. Conclusively, OXPHOS disruption and the associated disruptive mitochondrial energy metabolism are linked to the pathogenesis of DBA.

The Role of Non-Coding RNAs in Myelodysplastic Neoplasms

Myelodysplastic syndromes or neoplasms (MDS) are a heterogeneous group of myeloid clonal disorders characterized by peripheral blood cytopenias, blood and marrow cell dysplasia, and increased risk of evolution to acute myeloid leukemia (AML). Non-coding RNAs, especially microRNAs and long non-coding RNAs, serve as regulators of normal and malignant hematopoiesis and have been implicated in carcinogenesis. This review presents a comprehensive summary of the biology and role of non-coding RNAs, including the less studied circRNA, siRNA, piRNA, and snoRNA as potential prognostic and/or predictive biomarkers or therapeutic targets in MDS.

Pseudouridine-modified tRNA fragments repress aberrant protein synthesis and predict leukaemic progression in myelodysplastic syndrome

Transfer RNA-derived fragments (tRFs) are emerging small noncoding RNAs that, although commonly altered in cancer, have poorly defined roles in tumorigenesis¹. Here we show that pseudouridylation (Ψ) of a stem cell-enriched tRF subtype², mini tRFs containing a 5′ terminal oligoguanine (mTOG), selectively inhibits aberrant protein synthesis programmes, thereby promoting engraftment and differentiation of haematopoietic stem and progenitor cells (HSPCs) in patients with myelodysplastic syndrome (MDS). Building on evidence that mTOG-Ψ targets polyadenylate-binding protein cytoplasmic 1 (PABPC1), we employed isotope exchange proteomics to reveal critical interactions between mTOG and functional RNA-recognition motif (RRM) domains of PABPC1. Mechanistically, this hinders the recruitment of translational co-activator PABPC1-interacting protein 1 (PAIP1)³ and strongly represses the translation of transcripts sharing pyrimidine-enriched sequences (PES) at the 5′ untranslated region (UTR), including 5′ terminal oligopyrimidine tracts (TOP) that encode protein machinery components and are frequently altered in cancer⁴. Significantly, mTOG dysregulation leads to aberrantly increased translation of 5′ PES messenger RNA (mRNA) in malignant MDS-HSPCs and is clinically associated with leukaemic transformation and reduced patient survival. These findings define a critical role for tRFs and Ψ in difficult-to-treat subsets of MDS characterized by high risk of progression to acute myeloid leukaemia (AML).

Bellodi, Dimitriou and colleagues report that pseudouridine-modified transfer-RNA fragments modulate the translation of transcripts sharing pyrimidine-enriched sequences at their 5′ untranslated regions and their dysregulation impacts myelodysplastic syndrome pathogenesis.

Riddle of the Sphinx: Emerging Role of Transfer RNAs in Human Cancer

The dysregulation of transfer RNA (tRNA) expression contributes to the diversity of proteomics, heterogeneity of cell populations, and instability of the genome, which may be related to human cancer susceptibility. However, the relationship between tRNA dysregulation and cancer susceptibility remains elusive because the landscape of cancer-associated tRNAs has not been portrayed yet. Furthermore, the molecular mechanisms of tRNAs involved in tumorigenesis and cancer progression have not been systematically understood. In this review, we detail current knowledge of cancer-related tRNAs and comprehensively summarize the basic characteristics and functions of these tRNAs, with a special focus on their role and involvement in human cancer. This review bridges the gap between tRNAs and cancer and broadens our understanding of their relationship, thus providing new insights and strategies to improve the potential clinical applications of tRNAs for cancer diagnosis and therapy.

EV Translational Horizons as Viewed Across the Complex Landscape of Liquid Biopsies

Extracellular Vesicle (EV)-based diagnostic and therapeutic tools are an area of intensive study and substantial promise, but EVs as liquid biopsies have advanced years ahead of EVs as therapeutic tools. EVs are emerging as a promising approach for detecting tumors, evaluating the molecular profiles of known disease, and monitoring treatment responses. Although correlative assays based on liquid biopsies are already having an impact on translational studies and clinical practice, much remains to be learned before these assays will be optimized for clinical correlations, functional biological studies, and therapeutic use. What follows is an overview of current evidence supporting the investigation and use of liquid biopsies, organized by specific liquid biopsy components available for analysis, along with a summary of what challenges must be overcome before these assays will provide functional biological insights into the pathogenesis and treatment of disease. The same challenges must also be overcome before it will be feasible to measure and monitor the dosing, distribution, pharmacokinetics, and delivery of EV therapeutics and their cargo in complex biofluids where EVs and circulate with and are co-isolated with a number of other nanoscale materials, including lipoproteins (LPPs), ribonucleoprotein complexes (RNPs), and cell free nucleic acids (cfNA).

DBtRend: A Web-Server of tRNA Expression Profiles from Small RNA Sequencing Data in Humans

Transfer RNA (tRNA), a key component of the translation machinery, plays critical roles in stress conditions and various diseases. While knowledge regarding the importance of tRNA function is increasing, its biological roles are still not well understood. There is currently no comprehensive database or web server providing the expression landscape of tRNAs across a variety of human tissues and diseases. Here, we constructed a user-friendly and interactive database, DBtRend, which provides a profile of mature tRNA expression across various biological conditions by reanalyzing the small RNA or microRNA sequencing data from the Cancer Genome Atlas (TCGA) and NCBI's Gene Expression Omnibus (GEO) in humans. Users can explore not only the expression values of mature individual tRNAs in the human genome, but also those of isodecoders and isoacceptors based on our specific pipelines. DBtRend provides the expressed patterns of tRNAs, the differentially expressed tRNAs in different biological conditions, and the information of samples or patients, tissue types, and molecular subtype of cancers. The database is expected to help researchers interested in functional discoveries of tRNAs.

Emerging Functions for snoRNAs and snoRNA-Derived Fragments

The widespread implementation of mass sequencing has revealed a diverse landscape of small RNAs derived from larger precursors. Whilst many of these are likely to be byproducts of degradation, there are nevertheless metabolically stable fragments derived from tRNAs, rRNAs, snoRNAs, and other non-coding RNA, with a number of examples of the production of such fragments being conserved across species. Coupled with specific interactions to RNA-binding proteins and a growing number of experimentally reported examples suggesting function, a case is emerging whereby the biological significance of small non-coding RNAs extends far beyond miRNAs and piRNAs. Related to this, a similarly complex picture is emerging of non-canonical roles for the non-coding precursors, such as for snoRNAs that are also implicated in such areas as the silencing of gene expression and the regulation of alternative splicing. This is in addition to a body of literature describing snoRNAs as an additional source of miRNA-like regulators. This review seeks to highlight emerging roles for such non-coding RNA, focusing specifically on “new” roles for snoRNAs and the small fragments derived from them.

Transfer RNA-derived fragment tRF-28-QSZ34KRQ590K in plasma exosomes may be a potential biomarker for atopic dermatitis in pediatric patients

Atopic dermatitis (AD) is a common chronic relapsing inflammatory disease. There is substantial evidence suggesting that noncoding RNAs have indispensable roles in the pathogenesis of AD. Exosomal transfer RNA-derived fragments (tRFs) have been identified as potential biomarkers for various disorders. However, the role of tRFs in AD has remained to be elucidated, which was thus the aim of the present study. Plasma samples from 23 pediatric patients with AD and 23 healthy controls were collected. Exosomes were successfully isolated from plasma according to the manufacturer's protocol. Small RNA sequencing was performed in 3 patients with AD and 3 controls, and 135 significantly differentially expressed plasma exosomal tRFs were identified, including 36 upregulated and 99 downregulated tRFs. Using the miRanda and RNAhybrid databases, 58,227 target genes of these 135 differentially expressed tRFs were predicted. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses suggested that these target genes of tRFs are involved in multiple functions and pathways associated with AD. Downregulation of tRF-28-QSZ34KRQ590K and tRF-33-Q99P9P9NH57SD3 were validated in 20 patients with AD and 20 controls by reverse transcription-quantitative PCR and tRF-28-QSZ34KRQ590K exhibited significance in the receiver operating characteristic curve analysis. The present study was the first to provide a tRF profile in AD and implied that plasma exosomal tRF-28-QSZ34KRQ590K may be a potential biomarker for pediatric patients with AD. The present study enhanced the understanding of the pathogenesis of AD and provided a novel marker for the diagnosis and targeted treatment of AD.

tRFTars: predicting the targets of tRNA-derived fragments

Background

tRNA-derived fragments (tRFs) are 14–40-nucleotide-long, small non-coding RNAs derived from specific tRNA cleavage events with key regulatory functions in many biological processes. Many studies have shown that tRFs are associated with Argonaute (AGO) complexes and inhibit gene expression in the same manner as miRNAs. However, there are currently no tools for accurately predicting tRF target genes.

Methods

We used tRF-mRNA pairs identified by crosslinking, ligation, and sequencing of hybrids (CLASH) and covalent ligation of endogenous AGO-bound RNAs (CLEAR)-CLIP to assess features that may participate in tRF targeting, including the sequence context of each site and tRF-mRNA interactions. We applied genetic algorithm (GA) to select key features and support vector machine (SVM) to construct tRF prediction models.

Results

We first identified features that globally influenced tRF targeting. Among these features, the most significant were the minimum free folding energy (MFE), position 8 match, number of bases paired in the tRF-mRNA duplex, and length of the tRF, which were consistent with previous findings. Our constructed model yielded an area under the receiver operating characteristic (ROC) curve (AUC) = 0.980 (0.977–0.983) in the training process and an AUC = 0.847 (0.83–0.861) in the test process. The model was applied to all the sites with perfect Watson–Crick complementarity to the seed in the 3′ untranslated region (3′-UTR) of the human genome. Seven of nine target/nontarget genes of tRFs confirmed by reporter assay were predicted. We also validated the predictions via quantitative real-time PCR (qRT-PCR). Thirteen potential target genes from the top of the predictions were significantly down-regulated at the mRNA levels by overexpression of the tRFs (tRF-3001a, tRF-3003a or tRF-3009a).

Conclusions

Predictions can be obtained online, tRFTars, freely available at http://trftars.cmuzhenninglab.org:3838/tar/, which is the first tool to predict targets of tRFs in humans with a user-friendly interface.

On the expanding roles of tRNA fragments in modulating cell behavior

The fragments that derive from transfer RNAs (tRNAs) are an emerging category of regulatory RNAs. Known as tRFs, these fragments were reported for the first time only a decade ago, making them a relatively recent addition to the ever-expanding pantheon of non-coding RNAs. tRFs are short, 16-35 nucleotides (nts) in length, and produced through cleavage of mature and precursor tRNAs at various positions. Both cleavage positions and relative tRF abundance depend strongly on context, including the tissue type, tissue state, and disease, as well as the sex, population of origin, and race/ethnicity of an individual. These dependencies increase the urgency to understand the regulatory roles of tRFs. Such efforts are gaining momentum, and comprise experimental and computational approaches. System-level studies across many tissues and thousands of samples have produced strong evidence that tRFs have important and multi-faceted roles. Here, we review the relevant literature on tRF biology in higher organisms, single cell eukaryotes, and prokaryotes.

Circulating Small Noncoding RNAs Have Specific Expression Patterns in Plasma and Extracellular Vesicles in Myelodysplastic Syndromes and Are Predictive of Patient Outcome

Myelodysplastic syndromes (MDS) are hematopoietic stem cell disorders with large heterogeneity at the clinical and molecular levels. As diagnostic procedures shift from bone marrow biopsies towards less invasive techniques, circulating small noncoding RNAs (sncRNAs) have become of particular interest as potential novel noninvasive biomarkers of the disease. We aimed to characterize the expression profiles of circulating sncRNAs of MDS patients and to search for specific RNAs applicable as potential biomarkers. We performed small RNA-seq in paired samples of total plasma and plasma-derived extracellular vesicles (EVs) obtained from 42 patients and 17 healthy controls and analyzed the data with respect to the stage of the disease, patient survival, response to azacitidine, mutational status, and RNA editing. Significantly higher amounts of RNA material and a striking imbalance in RNA content between plasma and EVs (more than 400 significantly deregulated sncRNAs) were found in MDS patients compared to healthy controls. Moreover, the RNA content of EV cargo was more homogeneous than that of total plasma, and different RNAs were deregulated in these two types of material. Differential expression analyses identified that many hematopoiesis-related miRNAs (e.g., miR-34a, miR-125a, and miR-150) were significantly increased in MDS and that miRNAs clustered on 14q32 were specifically increased in early MDS. Only low numbers of circulating sncRNAs were significantly associated with somatic mutations in the SF3B1 or DNMT3A genes. Survival analysis defined a signature of four sncRNAs (miR-1237-3p, U33, hsa_piR_019420, and miR-548av-5p measured in EVs) as the most significantly associated with overall survival (HR = 5.866, p < 0.001). In total plasma, we identified five circulating miRNAs (miR-423-5p, miR-126-3p, miR-151a-3p, miR-125a-5p, and miR-199a-3p) whose combined expression levels could predict the response to azacitidine treatment. In conclusion, our data demonstrate that circulating sncRNAs show specific patterns in MDS and that their expression changes during disease progression, providing a rationale for the potential clinical usefulness of circulating sncRNAs in MDS prognosis. However, monitoring sncRNA levels in total plasma or in the EV fraction does not reflect one another, instead, they seem to represent distinctive snapshots of the disease and the data should be interpreted circumspectly with respect to the type of material analyzed.

Emerging roles of novel small non-coding regulatory RNAs in immunity and cancer

The term small non-coding RNAs (ncRNAs) refers to all those RNAs that even without encoding for a protein, can play important functional roles. Transfer RNA and ribosomal RNA-derived fragments (tRFs and rRFs, respectively) are an emerging class of ncRNAs originally considered as simple degradation products, which though play important roles in stress responses, signalling, or gene expression. They control all levels of gene expression regulating transcription and translation and affecting RNA processing and maturation. They have been linked to pivotal cellular processes such as self-renewal, differentiation, and proliferation. For this reason, mis-regulation of this novel class of ncRNAs can lead to various pathological processes such as neurodegenerative and development diseases, metabolism and immune system disorders, and cancer. In this review, we summarise the classification, biogenesis, and functions of tRFs and rRFs with a special focus on their role in immunity and cancer.

Translational efficiency across healthy and tumor tissues is proliferation-related

Different tissues express genes with particular codon usage and anticodon tRNA repertoires. However, the codon-anticodon co-adaptation in humans is not completely understood, nor is its effect on tissue-specific protein levels. Here, we first validated the accuracy of small RNA-seq for tRNA quantification across five human cell lines. We then analyzed the tRNA abundance of more than 8,000 tumor samples from TCGA, together with their paired mRNA-seq and proteomics data, to determine the Supply-to-Demand Adaptation. We thereby elucidate that the dynamic adaptation of the tRNA pool is largely related to the proliferative state across tissues. The distribution of such tRNA pools over the whole cellular translatome affects the subsequent translational efficiency, which functionally determines a condition-specific expression program both in healthy and tumor states. Furthermore, the aberrant translational efficiency of some codons in cancer, exemplified by ProCCA and GlyGGT, is associated with poor patient survival. The regulation of these tRNA profiles is partly explained by the tRNA gene copy numbers and their promoter DNA methylation.

tRic: a user-friendly data portal to explore the expression landscape of tRNAs in human cancers

Transfer RNAs (tRNAs) play critical roles in human cancer. Currently, no database provides the expression landscape and clinical relevance of tRNAs across a variety of human cancers. Utilizing miRNA-seq data from The Cancer Genome Atlas, we quantified the relative expression of tRNA genes and merged them into the codon level and amino level across 31 cancer types. The expression of tRNAs is associated with clinical features of patient smoking history and overall survival, and disease stage, subtype, and grade. We further analysed codon frequency and amino acid frequency for each protein coding gene and linked alterations of tRNA expression with protein translational efficiency. We include these data resources in a user-friendly data portal, tRic (tRNA in cancer, https://hanlab.uth.edu/tRic/ or http://bioinfo.life.hust.edu.cn/tRic/), which can be of significant interest to the research community.

Integrated Analysis of Mouse and Human Gastric Neoplasms Identifies Conserved microRNA Networks in Gastric Carcinogenesis

BACKGROUND & AIMS

microRNAs (miRNAs) are small noncoding RNAs that bind to the 3' untranslated regions of mRNAs to promote their degradation or block their translation. Mice with disruption of the trefoil factor 1 gene (Tff1) develop gastric neoplasms. We studied these mice to identify conserved miRNA networks involved in gastric carcinogenesis.

METHODS

We performed next-generation miRNA sequencing analysis of normal gastric tissues (based on histology) from patients without evidence of gastric neoplasm (n = 64) and from TFF1-knockout mice (n = 22). We validated our findings using 270 normal gastric tissues (including 61 samples from patients without evidence of neoplastic lesions) and 234 gastric tumor tissues from 3 separate cohorts of patients and from mice. We performed molecular and functional assays using cell lines (MKN28, MKN45, STKM2, and AGS cells), gastric organoids, and mice with xenograft tumors.

RESULTS

We identified 117 miRNAs that were significantly deregulated in mouse and human gastric tumor tissues compared with nontumor tissues. We validated changes in levels of 6 miRNAs by quantitative real-time polymerase chain reaction analyses of neoplastic gastric tissues from mice (n = 39) and 3 independent patient cohorts (n = 332 patients total). We found levels of MIR135B-5p, MIR196B-5p, and MIR92A-5p to be increased in tumor tissues, whereas levels of MIR143-3p, MIR204-5p, and MIR133-3p were decreased in tumor tissues. Levels of MIR143-3p were reduced not only in gastric cancer tissues but also in normal tissues adjacent to tumors in humans and low-grade dysplasia in mice. Transgenic expression of MIR143-3p in gastric cancer cell lines reduced their proliferation and restored their sensitivity to cisplatin. AGS cells with stable transgenic expression of MIR143-3p grew more slowly as xenograft tumors in mice than control AGS cells; tumor growth from AGS cells that expressed MIR143-3p, but not control cells, was sensitive to cisplatin. We identified and validated bromodomain containing 2 (BRD2) as a direct target of MIR143-3p; increased levels of BRD2 in gastric tumors was associated with shorter survival times for patients.

CONCLUSIONS

In an analysis of miRNA profiles of gastric tumors from mice and human patients, we identified a conserved signature associated with the early stages of gastric tumorigenesis. Strategies to restore MIR143-3p or inhibit BRD2 might be developed for treatment of gastric cancer.

Altered splicing and cytoplasmic levels of tRNA synthetases in SF3B1-mutant myelodysplastic syndromes as a therapeutic vulnerability

Myelodysplastic syndromes (MDS) are haematopoietic malignancies that are characterised by a heterogeneous clinical course. In recent years, sequencing efforts have uncovered recurrent somatic mutations within RNA splicing factors, including SF3B1, SRSF2, U2AF1 and ZRSR2. The most frequently mutated gene is SF3B1, mutated in 17% of MDS patients. While SF3B1 mutations and their effects on splicing have been well characterised, much remains to be explored about their more far-reaching effects on cellular homeostasis. Given that mRNA splicing and nuclear export are coordinated processes, we hypothesised that SF3B1 mutation might also affect export of certain mRNAs and that this may represent a targetable pathway for the treatment of SF3B1-mutant MDS. We used CRISPR/Cas9-genome editing to create isogenic cellular models. Comprehensive transcriptome and proteome profiling of these cells identified alterations in the splicing and export of components of the translational machinery, primarily tRNA synthetases, in response to the SF3B1 K700E mutation. While steady-state protein synthesis was unaffected, SF3B1 mutant cells were more sensitive to the clinically-relevant purine analogue, 8-azaguanine. In this study, we also demonstrated that 8-azaguanine affects splicing. Our results suggest that the simultaneous targeting of RNA metabolism and splicing by 8-azaguanine represents a therapeutic opportunity for SF3B1-mutant myelodysplastic syndromes.

Global analysis of tRNA and translation factor expression reveals a dynamic landscape of translational regulation in human cancers

The protein translational system, including transfer RNAs (tRNAs) and several categories of enzymes, plays a key role in regulating cell proliferation. Translation dysregulation also contributes to cancer development, though relatively little is known about the changes that occur to the translational system in cancer. Here, we present global analyses of tRNAs and three categories of enzymes involved in translational regulation in ~10,000 cancer patients across 31 cancer types from The Cancer Genome Atlas. By analyzing the expression levels of tRNAs at the gene, codon, and amino acid levels, we identified unequal alterations in tRNA expression, likely due to the uneven distribution of tRNAs decoding different codons. We find that overexpression of tRNAs recognizing codons with a low observed-over-expected ratio may overcome the translational bottleneck in tumorigenesis. We further observed overall overexpression and amplification of tRNA modification enzymes, aminoacyl-tRNA synthetases, and translation factors, which may play synergistic roles with overexpression of tRNAs to activate the translational systems across multiple cancer types.

Non-coding RNAs Function as Immune Regulators in Teleost Fish

Non-coding RNAs (ncRNAs) are functional RNA molecules that are transcribed from DNA but not translated into proteins. ncRNAs function as key regulators of gene expression and chromatin modification. Recently, the functional role of ncRNAs in teleost fish has been extensively studied. Teleost fish are a highly diverse group among the vertebrate lineage. Fish are also important in terms of aquatic ecosystem, food production and human life, being the source of animal proteins worldwide and models of biomedical research. However, teleost fish always suffer from the invasion of infectious pathogens including viruses and bacteria, which has resulted in a tremendous economic loss to the fishing industry worldwide. Emerging evidence suggests that ncRNAs, especially miRNAs and lncRNAs, may serve as important regulators in cytokine and chemokine signaling, antigen presentation, complement and coagulation cascades, and T cell response in teleost fish. In this review, we summarize current knowledge and understanding of the roles of both miRNAs and lncRNAs in immune regulation in teleost fish. Molecular mechanism insights into the function of ncRNAs in fish immune response may contribute to the development of potential biomarkers and therapeutic targets for the prevention and treatment of fish diseases.

Multi-perspective quality control of Illumina RNA sequencing data analysis

Quality control (QC) is a critical step in RNA sequencing (RNA-seq). Yet, it is often ignored or conducted on a limited basis. Here, we present a multi-perspective strategy for QC of RNA-seq experiments. The QC of RNA-seq can be divided into four related stages: (1) RNA quality, (2) raw read data (FASTQ), (3) alignment and (4) gene expression. We illustrate the importance of conducting QC at each stage of an RNA-seq experiment and demonstrate our recommended RNA-seq QC strategy. Furthermore, we discuss the major and often neglected quality issues associated with the three major types of RNA-seq: mRNA, total RNA and small RNA. This RNA-seq QC overview provides comprehensive guidance for researchers who conduct RNA-seq experiments.

The Landscape of Small Non-Coding RNAs in Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) is an operational term for breast cancers lacking targetable estrogen receptor expression and HER2 amplifications. TNBC is, therefore, inherently heterogeneous, and is associated with worse prognosis, greater rates of metastasis, and earlier onset. TNBC displays mutational and transcriptional diversity, and distinct mRNA transcriptional subtypes exhibiting unique biology. High-throughput sequencing has extended cancer research far beyond protein coding regions that include non-coding small RNAs, such as miRNA, isomiR, tRNA, snoRNAs, snRNA, yRNA, 7SL, and 7SK. In this study, we performed small RNA profiling of 26 TNBC cell lines, and compared the abundance of non-coding RNAs among the transcriptional subtypes of triple negative breast cancer. We also examined their co-expression pattern with corresponding mRNAs. This study provides a detailed description of small RNA expression in triple-negative breast cancer cell lines that can aid in the development of future biomarker and novel targeted therapies.

Current Research on Non-Coding Ribonucleic Acid (RNA)

Non-coding ribonucleic acid (RNA) has without a doubt captured the interest of biomedical researchers. The ability to screen the entire human genome with high-throughput sequencing technology has greatly enhanced the identification, annotation and prediction of the functionality of non-coding RNAs. In this review, we discuss the current landscape of non-coding RNA research and quantitative analysis. Non-coding RNA will be categorized into two major groups by size: long non-coding RNAs and small RNAs. In long non-coding RNA, we discuss regular long non-coding RNA, pseudogenes and circular RNA. In small RNA, we discuss miRNA, transfer RNA, piwi-interacting RNA, small nucleolar RNA, small nuclear RNA, Y RNA, single recognition particle RNA, and 7SK RNA. We elaborate on the origin, detection method, and potential association with disease, putative functional mechanisms, and public resources for these non-coding RNAs. We aim to provide readers with a complete overview of non-coding RNAs and incite additional interest in non-coding RNA research.

Discovery and function of transfer RNA-derived fragments and their role in disease

Until recently, transfer RNAs (tRNAs) were thought to function in protein translation only. However, recent findings demonstrate that both pre- and mature tRNAs can undergo endonucleolytic cleavage by different ribonucleases originating different types of small non-coding RNAs, known as tRNA-derived fragments (tRFs). tRFs are classified according to their origin and are implicated in various cellular processes, namely apoptosis, protein synthesis control, and RNA interference. Although their functions are still poorly understood, their mechanisms of action vary according to the tRF sub-type. Several tRFs have been associated with cancer, neurodegenerative disorders, and viral infections and growing evidence shows that they may constitute novel molecular targets for modulating pathological processes. Here, we recapitulate the current knowledge of tRF biology, highlight the known functions and mechanisms of action of the different sub-classes of tRFs and discuss their implications in human disease. WIREs RNA 2017, 8:e1423. doi: 10.1002/wrna.1423 For further resources related to this article, please visit the WIREs website.

Techniques for detecting chromosomal aberrations in myelodysplastic syndromes

Myelodysplastic syndromes (MDS) are a group of heterogeneous hematologic diseases. Chromosomal aberrations are important for the initiation, development, and progression of MDS. Detection of chromosomal abnormalities in MDS is important for categorization, risk stratification, therapeutic selection, and prognosis evaluation of the disease. Recent progress of multiple techniques has brought powerful molecular cytogenetic information to reveal copy number variation, uniparental disomy, and complex chromosomal aberrations in MDS. In this review, we will introduce some common chromosomal aberrations in MDS and their clinical significance. Then we will explain the application, advantages, and limitations of different techniques for detecting chromosomal abnormalities in MDS. The information in this review may be helpful for clinicians to select appropriate methods in patient-related decision making.

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.

MicroRNAs and tRNA-derived fragments predict the transformation of myelodysplastic syndromes to acute myeloid leukemia

Myelodysplastic syndromes (MDS) are clonal hematopoietic disorders of the elderly that carry an increased risk of progression to acute myeloid leukemia (AML). Since small non-coding RNAs (sRNAs), including microRNA (miRNAs), act as regulators of cellular differentiation, we hypothesized that changes to sRNAs might be implicated in the progression of MDS to AML. We conducted sRNA sequencing on three sets of patients: Group A (MDS patients who never progressed to AML); Group B (MDS patients who later progressed to an AML); and Group C (AML patients with myelodysplasia-related changes, including patients with a known preceding diagnosis of MDS). We identified five miRNAs that differentiated Groups A and B, independent of bone marrow blast percentage, including three members of the miR-181 family, as well as differential patterns of miRNA isoforms (isomiRs) and tDRs. Thus, we have identified sRNA biomarkers that predict MDS cases that are likely to progress to AML.

Comprehensive evaluation of extracellular small RNA isolation methods from serum in high throughput sequencing

Background

DNA and RNA fractions from whole blood, serum and plasma are increasingly popular analytes that are currently under investigation for their utility in the diagnosis and staging of disease. Small non-coding ribonucleic acids (sRNAs), specifically microRNAs (miRNAs) and their variant isoforms (isomiRs), and transfer RNA (tRNA)-derived small RNAs (tDRs) comprise a repertoire of molecules particularly promising in this regard.

Results

In this designed study, we compared the performance of various methods and kits for isolating circulating extracellular sRNAs (ex-sRNAs). ex-sRNAs from one healthy individual were isolated using five different isolation kits: Qiagen Circulating Nucleic Acid Kit, ThermoFisher Scientific Ambion TRIzol LS Reagent, Qiagen miRNEasy, QiaSymphony RNA extraction kit and the Exiqon MiRCURY RNA Isolation Kit. Each isolation method was repeated four times. A total of 20 small RNA sequencing (sRNAseq) libraries were constructed, sequenced and compared using a rigorous bioinformatics approach. The Circulating Nucleic Acid Kit had the greatest miRNA isolation variability, but had the lowest isolation variability for other RNA classes (isomiRs, tDRs, and other miscellaneous sRNAs (osRNA). However, the Circulating Nucleic Acid Kit consistently generated the fewest number of reads mapped to the genome, as compared to the best-performing method, Ambion TRIzol, which mapped 10% of the miRNAs, 7.2% of the tDRs and 23.1% of the osRNAs. The other methods performed intermediary, with QiaSymphony mapping 14% of the osRNAs, and miRNEasy mapping 4.6% of the tDRs and 2.9% of the miRNAs, achieving the second best kit performance rating overall.

Conclusions

In summary, each isolation kit displayed different performance characteristics that could be construed as biased or advantageous, depending upon the downstream application and number of samples that require processing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3470-z) contains supplementary material, which is available to authorized users.

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 micro-RNA expression signature for human NAFLD progression

BACKGROUND

The spectrum of nonalcoholic fatty liver disease (NAFLD) describes disease conditions deteriorating from nonalcoholic fatty liver (NAFL) to nonalcoholic steatohepatitis (NASH) to cirrhosis (CIR) to hepatocellular carcinoma (HCC). From a molecular and biochemical perspective, our understanding of the etiology of this disease is limited by the broad spectrum of disease presentations, the lack of a thorough understanding of the factors contributing to disease susceptibility, and ethical concerns related to repeat sampling of the liver. To better understand the factors associated with disease progression, we investigated by next-generation RNA sequencing the altered expression of microRNAs (miRNAs) in liver biopsies of class III obese subjects (body mass index ≥40 kg/m(2)) biopsied at the time of elective bariatric surgery.

METHODS

Clinical characteristics and unbiased RNA expression profiles for 233 miRs, 313 transfer RNAs (tRNAs), and 392 miscellaneous small RNAs (snoRNAs, snRNAs, rRNAs) were compared among 36 liver biopsy specimens stratified by disease severity.

RESULTS

The abundances of 3 miRNAs that were found to be differentially regulated (miR-301a-3p and miR-34a-5p increased and miR-375 decreased) with disease progression were validated by RT-PCR. No tRNAs or miscellaneous RNAs were found to be associated with disease severity. Similar patterns of increased miR-301a and decreased miR-375 expression were observed in 134 hepatocellular carcinoma (HCC) samples deposited in The Cancer Genome Atlas (TCGA).

CONCLUSIONS

Our analytical results suggest that NAFLD severity is associated with a specific pattern of altered hepatic microRNA expression that may drive the hallmark of this disorder: altered lipid and carbohydrate metabolism. The three identified miRNAs can potentially be used as biomarkers to access the severity of NAFLD. The persistence of this miRNA expression pattern in an external validation cohort of HCC samples suggests that specific microRNA expression patterns may permit and/or sustain NAFLD development to HCC.

Practicability of mitochondrial heteroplasmy detection through an Illumina genotyping array

MOTIVATION

High-throughput genomic data often contain unexpected information that can be mined for alternative applications. Despite the rise of high-throughput sequencing, Illumina genotyping arrays remain a driving force in large scale genetic and epidemiology studies. By processing and analyzing genotyping data of over 100,000 samples genotyped on Illumina genotyping arrays, we discovered evidence that indicates that mitochondrial heteroplasmy can be estimated from the fluorescence intensity data of the array. To verify our hypothesis, we conducted a sequencing validation study.

RESULT

Mitochondrial DNA targeted sequencing was performed on three samples that had been genotyped using the Illumina exome genotyping array. In each sample chosen, one heteroplasmy target was identified from the genotyping array, and sequencing data verified all three putative heteroplasmic sites. The estimated heteroplasmy level difference between that estimated from the genotyping fluorescence intensity and that directly measured from sequencing was 3.2% on average. Our analysis showed that an Illumina genotyping array can accurately and reliably estimate high-level heteroplasmy (>40%); however, intensity data from a genotyping array is not suitable for estimating low level heteroplasmy (<25%).

Genome-wide profiling of transfer RNAs and their role as novel prognostic markers for breast cancer

Transfer RNAs (tRNAs, key molecules in protein synthesis) have not been investigated as potential prognostic markers in breast cancer (BC), despite early findings of their dysregulation and diagnostic potential. We aim to comprehensively profile tRNAs from breast tissues and to evaluate their role as prognostic markers (Overall Survival, OS and Recurrence Free Survival, RFS). tRNAs were profiled from 11 normal breast and 104 breast tumor tissues using next generation sequencing. We adopted a Case-control (CC) and Case-Only (CO) association study designs. Risk scores constructed from tRNAs were subjected to univariate and multivariate Cox-proportional hazards regression to investigate their prognostic value. Of the 571 tRNAs profiled, 76 were differentially expressed (DE) and three were significant for OS in the CC approach. We identified an additional 11 tRNAs associated with OS and 14 tRNAs as significant for RFS in the CO approach, indicating that CC alone may not capture all discriminatory tRNAs in prognoses. In both the approaches, the risk scores were significant in the multivariate analysis as independent prognostic factors, and patients belonging to high-risk group were associated with poor prognosis. Our results confirmed global up-regulation of tRNAs in BC and identified tRNAs as potential novel prognostic markers for BC.

Translating RNA sequencing into clinical diagnostics: opportunities and challenges

RNA-based measurements have the potential for application across diverse areas of human health, including disease diagnosis, prognosis and therapeutic selection. Current clinical applications include infectious diseases, cancer, transplant medicine and fetal monitoring.

RNA sequencing (RNA-seq) allows for the detection of a wide variety of RNA species, including mRNA, non-coding RNA, pathogen RNA, chimeric gene fusions, transcript isoforms and splice variants, and provides the capability to quantify known, pre-defined RNA species and rare RNA transcript variants within a sample. In addition to differential expression and detection of novel transcripts, RNA-seq also supports the detection of mutations and germline variation for hundreds to thousands of expressed genetic variants, facilitating assessment of allele-specific expression of these variants.

Circulating RNAs and small regulatory RNAs, such as microRNAs, are very stable. These RNA species are vigorously being tested for their potential as biomarkers. However, there are currently few agreed upon methods for isolation or quantitative measurements and a current lack of quality controls that can be used to test platform accuracy and sample preparation quality.

Analytical, bioinformatic and regulatory challenges exist, and ongoing efforts toward the establishment of benchmark standards, assay optimization for clinical conditions and demonstration of assay reproducibility are required to expand the clinical utility of RNA-seq.

RNA sequencing (RNA-seq) is a powerful approach for comprehensive analyses of transcriptomes. This Review describes the widespread potential applications of RNA-seq in clinical medicine, such as detecting disease-associated mutations and gene expression disruptions, as well as characteristic non-coding RNAs, circulating extracellular RNAs or pathogen RNAs. The authors also highlight the challenges in adopting RNA-seq routinely into clinical practice.

With the emergence of RNA sequencing (RNA-seq) technologies, RNA-based biomolecules hold expanded promise for their diagnostic, prognostic and therapeutic applicability in various diseases, including cancers and infectious diseases. Detection of gene fusions and differential expression of known disease-causing transcripts by RNA-seq represent some of the most immediate opportunities. However, it is the diversity of RNA species detected through RNA-seq that holds new promise for the multi-faceted clinical applicability of RNA-based measures, including the potential of extracellular RNAs as non-invasive diagnostic indicators of disease. Ongoing efforts towards the establishment of benchmark standards, assay optimization for clinical conditions and demonstration of assay reproducibility are required to expand the clinical utility of RNA-seq.