Improved TGIRT-seq methods for comprehensive transcriptome profiling with decreased adapter dimer formation and bias correction

RNA methylation: A new promising biomaker in cancer liquid biopsy

RNA methylation is a vital epigenetic modification that regulates gene expression by influencing RNA processes such as transcription, degradation, translation, and transport. Aberrant methylation, including modifications like m6A, m5C, m1A, m7G, and m3C, is closely linked to tumorigenesis and progression. Liquid biopsy, a non-invasive technique analyzing tumor markers in body fluids, offers significant potential for early diagnosis and dynamic monitoring. In this context, RNA methylation, due to its tumor-specific properties, is emerging as a valuable marker. However, significant challenges remain in its clinical application. This review explores the roles of RNA methylation in cancer, recent advances in detection technologies, and its potential as a liquid biopsy marker in tumor management. It highlights its promising applications in cancer diagnosis, prognosis, and personalized treatment in the era of precision oncology.

Optimized identification and characterization of small RNAs with PANDORA-seq

Small noncoding RNAs (sncRNAs) are a diverse group of RNAs including small interfering RNAs, microRNAs, PIWI-interacting RNAs and RNAs derived from structured RNAs such as transfer RNAs, ribosomal RNAs and others. These sncRNAs have varied termini and RNA modifications, which can interfere with adaptor ligation and reverse transcription during cDNA library construction, hindering detection of many types of sncRNA by standard small RNA sequencing methods. To address this limitation, PANDORA sequencing introduces a refined methodology. The procedure includes sequential enzymatic treatments of size-selected RNAs with T4PNK and AlkB, which effectively circumvent the challenges presented by the ligation-blocking termini and reverse transcription-blocking RNA modifications, followed by tailored small RNA library construction protocols and deep sequencing. The obtained datasets are analyzed with the SPORTS pipeline, which can comprehensively analyze various types of sncRNA beyond the traditionally studied classes, to include those derived from various parental RNAs (for example, from transfer RNA and ribosomal RNA), as well as output the locations on the parental RNA from which these sncRNAs are derived. The entire protocol takes ~7 d, depending on the sample size and sequencing turnaround time. PANDORA sequencing provides a transformative tool to further our understanding of the expanding small RNA universe and to explore the uncharted functions of sncRNAs.

All the sites we cannot see: Sources and mitigation of false negatives in RNA modification studies

RNA modifications are essential for human health - too much or too little of them leads to serious illnesses ranging from neurodevelopmental disorders to cancer. Technical advances in RNA modification sequencing are beginning to uncover the RNA targets of diverse RNA-modifying enzymes that are dysregulated in disease. However, the emerging transcriptome-wide maps of modified nucleosides installed by these enzymes should be considered as first drafts. In particular, a range of technical artefacts lead to false negatives - modified sites that are overlooked owing to technique-dependent, and often sequence-context-specific, 'blind spots'. In this Review, we discuss potential sources of false negatives in sequencing-based RNA modification maps, propose mitigation strategies and suggest guidelines for transparent reporting of sensitivity to detect modified sites in profiling studies. Important considerations for recognition and avoidance of false negatives include assessment and reporting of position-specific sequencing depth, identification of protocol-dependent RNA capture biases and applying controls for false negatives as well as for false positives. Despite their limitations, emerging maps of RNA modifications reveal exciting and largely uncharted potential for post-transcriptional control of all aspects of RNA function.

Sources of non-uniform coverage in short-read RNA-Seq data

The origin of several normal cellular functions and related abnormalities can be traced back to RNA splicing. As such, RNA splicing is currently the focus of a vast array of studies. To quantify the transcriptome, short-read RNA-Seq remains the standard assay. The primary technical artifact of RNASeq library prep, which severely interferes with analysis, is extreme non-uniformity in coverage across transcripts. This non-uniformity is present in both bulk and single-cell RNA-Seq and is observed even when the sample contains only full-length transcripts. This issue dramatically affects the accuracy of isoform-level quantification of multi-isoform genes. Understanding the sources of this non-uniformity is critical to developing improved protocols and analysis methods. Here, we explore eight potential sources of non-uniformity. We demonstrate that it cannot be explained by one factor alone. We performed targeted experiments to investigate the effect of fragment length, PCR ramp rate, and ribosomal depletion. We assessed existing data sets with varying sample quality, PCR cycle number, reverse transcriptase, and technical or biological replicates. We found evidence that interference of reverse transcription by secondary structure is unlikely to be the major contributing factor, that rRNA pull-down methods do not cause non-uniformity, that PCR ramp rate does not substantially impact non-uniformity, and that shorter fragments do not reduce non-uniformity. All these findings contradict prior publications or recommendations.

Improved precision, sensitivity, and adaptability of ordered two-template relay cDNA library preparation for RNA sequencing

Sequencing RNAs that are biologically processed or degraded to less than ∼100 nt typically involves multistep, low-yield protocols with bias and information loss inherent to ligation and/or polynucleotide tailing. We recently introduced ordered two-template relay (OTTR), a method that captures obligatorily end-to-end sequences of input molecules and, in the same reverse transcription step, also appends 5' and 3' sequencing adapters of choice. OTTR has been thoroughly benchmarked for optimal production of microRNA, tRNA and tRNA fragments, and ribosome-protected mRNA footprint libraries. Here we sought to characterize, quantify, and ameliorate any remaining bias or imprecision in the end-to-end capture of RNA sequences. We introduce new metrics for the evaluation of sequence capture and use them to optimize reaction buffers, reverse transcriptase sequence, adapter oligonucleotides, and overall workflow. Modifications of the reverse transcriptase and adapter oligonucleotides increased the 3' and 5' end-precision of sequence capture and minimized overall library bias. Improvements in recombinant expression and purification of the truncated Bombyx mori R2 reverse transcriptase used in OTTR reduced nonproductive sequencing reads by minimizing bacterial nucleic acids that compete with low-input RNA molecules for cDNA synthesis, such that with miRNA input of 3 pg (<1 fmol), fewer than 10% of sequencing reads are bacterial nucleic acid contaminants. We also introduce a rapid, automation-compatible OTTR protocol that enables gel-free, length-agnostic enrichment of cDNA duplexes from unwanted adapter-only side products. Overall, this work informs considerations for unbiased end-to-end capture and annotation of RNAs independent of their sequence, structure, or posttranscriptional modifications.

tRNA regulation and amino acid usage bias reflect a coordinated metabolic adaptation in Plasmodium falciparum

An adaptive feature of malaria-causing parasites is the digestion of host hemoglobin (HB) to acquire amino acids (AAs). Here, we describe a link between nutrient availability and translation dependent regulation of gene expression as an adaptive strategy. We show that tRNA expression in Plasmodium falciparum does not match the decoding need expected for optimal translation. A subset of tRNAs decoding AAs that are insufficiently provided by HB are lowly expressed, wherein the abundance of a protein-coding transcript is negatively correlated with the decoding requirement of these tRNAs. Proliferation-related genes have evolved a high requirement of these tRNAs, thereby proliferation can be modulated by repressing protein synthesis of these genes during nutrient stress. We conclude that the parasite modulates translation elongation by maintaining a discordant tRNA profile to exploit variations in AA-composition among genes as an adaptation strategy. This study exemplifies metabolic adaptation as an important driving force for protein evolution.

Improved precision, sensitivity, and adaptability of Ordered Two-Template Relay cDNA library preparation for RNA sequencing

Sequencing RNAs that are biologically processed or degraded to less than ~100 nucleotides typically involves multi-step, low-yield protocols with bias and information loss inherent to ligation and/or polynucleotide tailing. We recently introduced Ordered Two-Template Relay (OTTR), a method that captures obligatorily end-to-end sequences of input molecules and, in the same reverse transcription step, also appends 5' and 3' sequencing adapters of choice. OTTR has been thoroughly benchmarked for optimal production of microRNA, tRNA and tRNA fragments, and ribosome-protected mRNA footprint libraries. Here we sought to characterize, quantify, and ameliorate any remaining bias or imprecision in the end-to-end capture of RNA sequences. We introduce new metrics for the evaluation of sequence capture and use them to optimize reaction buffers, reverse transcriptase sequence, adapter oligonucleotides, and overall workflow. Modifications of the reverse transcriptase and adapter oligonucleotides increased the 3' and 5' end-precision of sequence capture and minimized overall library bias. Improvements in recombinant expression and purification of the truncated Bombyx mori R2 reverse transcriptase used in OTTR reduced non-productive sequencing reads by minimizing bacterial nucleic acids that compete with low-input RNA molecules for cDNA synthesis, such that with miRNA input of 3 picograms (less than 1 fmol), fewer than 10% of sequencing reads are bacterial nucleic acid contaminants. We also introduce a rapid, automation-compatible OTTR protocol that enables gel-free, length-agnostic enrichment of cDNA duplexes from unwanted adapter-only side products. Overall, this work informs considerations for unbiased end-to-end capture and annotation of RNAs independent of their sequence, structure, or post-transcriptional modifications.

Biochemical and structural insights into a 5' to 3' RNA ligase reveal a potential role in tRNA ligation

ATP-grasp superfamily enzymes contain a hand-like ATP-binding fold and catalyze a variety of reactions using a similar catalytic mechanism. More than 30 protein families are categorized in this superfamily, and they are involved in a plethora of cellular processes and human diseases. Here, we identify C12orf29 (RLIG1) as an atypical ATP-grasp enzyme that ligates RNA. Human RLIG1 and its homologs autoadenylate on an active site Lys residue as part of a reaction intermediate that specifically ligates RNA halves containing a 5'-phosphate and a 3'-hydroxyl. RLIG1 binds tRNA in cells and can ligate tRNA within the anticodon loop in vitro. Transcriptomic analyses of Rlig1 knockout mice revealed significant alterations in global tRNA levels in the brains of female mice, but not in those of male mice. Furthermore, crystal structures of a RLIG1 homolog from Yasminevirus bound to nucleotides revealed a minimal and atypical RNA ligase fold with a conserved active site architecture that participates in catalysis. Collectively, our results identify RLIG1 as an RNA ligase and suggest its involvement in tRNA biology.

A ligation-independent sequencing method reveals tRNA-derived RNAs with blocked 3' termini

Despite the numerous sequencing methods available, the diversity in RNA size and chemical modification makes it difficult to capture all RNAs in a cell. We developed a method that combines quasi-random priming with template switching to construct sequencing libraries from RNA molecules of any length and with any type of 3' modifications, allowing for the sequencing of virtually all RNA species. Our ligation-independent detection of all types of RNA (LIDAR) is a simple, effective tool to identify and quantify all classes of coding and non-coding RNAs. With LIDAR, we comprehensively characterized the transcriptomes of mouse embryonic stem cells, neural progenitor cells, mouse tissues, and sperm. LIDAR detected a much larger variety of tRNA-derived RNAs (tDRs) compared with traditional ligation-dependent sequencing methods and uncovered tDRs with blocked 3' ends that had previously escaped detection. Therefore, LIDAR can capture all RNAs in a sample and uncover RNA species with potential regulatory functions.

Human cells contain myriad excised linear intron RNAs with links to gene regulation and potential utility as biomarkers

A previous study using Thermostable Group II Intron Reverse Transcriptase sequencing (TGIRT-seq) found human plasma contains short (≤300 nt) structured full-length excised linear intron (FLEXI) RNAs with potential to serve as blood-based biomarkers. Here, TGIRT-seq identified >9,000 different FLEXI RNAs in human cell lines, including relatively abundant FLEXIs with cell-type-specific expression patterns. Analysis of public CLIP-seq datasets identified 126 RNA-binding proteins (RBPs) that have binding sites within the region corresponding to the FLEXI or overlapping FLEXI splice sites in pre-mRNAs, including 53 RBPs with binding sites for ≥30 different FLEXIs. These included splicing factors, transcription factors, a chromatin remodeling protein, cellular growth regulators, and proteins with cytoplasmic functions. Analysis of ENCODE datasets identified subsets of these RBPs whose knockdown impacted FLEXI host gene mRNA levels or proximate alternative splicing, indicating functional interactions. Hierarchical clustering identified six subsets of RBPs whose FLEXI binding sites were co-enriched in six subsets of functionally related host genes: AGO1-4 and DICER, including but not limited to agotrons or mirtron pre-miRNAs; DKC1, NOLC1, SMNDC1, and AATF (Apoptosis Antagonizing Transcription Factor), including but not limited to snoRNA-encoding FLEXIs; two subsets of alternative splicing factors; and two subsets that included RBPs with cytoplasmic functions (e.g., LARP4, PABPC4, METAP2, and ZNF622) together with regulatory proteins. Cell fractionation experiments showed cytoplasmic enrichment of FLEXI RNAs with binding sites for RBPs with cytoplasmic functions. The subsets of host genes encoding FLEXIs with binding sites for different subsets of RBPs were co-enriched with non-FLEXI other short and long introns with binding sites for the same RBPs, suggesting overarching mechanisms for coordinately regulating expression of functionally related genes. Our findings identify FLEXIs as a previously unrecognized large class of cellular RNAs and provide a comprehensive roadmap for further analyzing their biological functions and the relationship of their RBPs to cellular regulatory mechanisms.

Senataxin deficiency disrupts proteostasis through nucleolar ncRNA-driven protein aggregation

Senataxin is an evolutionarily conserved RNA-DNA helicase involved in DNA repair and transcription termination that is associated with human neurodegenerative disorders. Here, we investigated whether Senataxin loss affects protein homeostasis based on previous work showing R-loop-driven accumulation of DNA damage and protein aggregates in human cells. We find that Senataxin loss results in the accumulation of insoluble proteins, including many factors known to be prone to aggregation in neurodegenerative disorders. These aggregates are located primarily in the nucleolus and are promoted by upregulation of non-coding RNAs expressed from the intergenic spacer region of ribosomal DNA. We also map sites of R-loop accumulation in human cells lacking Senataxin and find higher RNA-DNA hybrids within the ribosomal DNA, peri-centromeric regions, and other intergenic sites but not at annotated protein-coding genes. These findings indicate that Senataxin loss affects the solubility of the proteome through the regulation of transcription-dependent lesions in the nucleus and the nucleolus.

Full-length tRNAs lacking a functional CCA tail are selectively sorted into the lumen of extracellular vesicles

Small extracellular vesicles (sEVs) are heterogenous lipid membrane particles typically less than 200 nm in size and secreted by most cell types either constitutively or upon activation signals. sEVs isolated from biofluids contain RNAs, including small non-coding RNAs (ncRNAs), that can be either encapsulated within the EV lumen or bound to the EV surface. EV-associated microRNAs (miRNAs) are, despite a relatively low abundance, extensively investigated for their selective incorporation and their role in cell-cell communication. In contrast, the sorting of highly-structured ncRNA species is understudied, mainly due to technical limitations of traditional small RNA sequencing protocols. Here, we adapted ALL-tRNAseq to profile the relative abundance of highly structured and potentially methylated small ncRNA species, including transfer RNAs (tRNAs), small nucleolar RNAs (snoRNAs), and Y RNAs in bulk EV preparations. We determined that full-length tRNAs, typically 75 to 90 nucleotides in length, were the dominant small ncRNA species (>60% of all reads in the 18-120 nucleotides size-range) in all cell culture-derived EVs, as well as in human plasma-derived EV samples, vastly outnumbering 21 nucleotides-long miRNAs. Nearly all EV-associated tRNAs were protected from external RNAse treatment, indicating a location within the EV lumen. Strikingly, the vast majority of luminal-sorted, full-length, nucleobase modification-containing EV-tRNA sequences, harbored a dysfunctional 3' CCA tail, 1 to 3 nucleotides truncated, rendering them incompetent for amino acid loading. In contrast, in non-EV associated extracellular particle fractions (NVEPs), tRNAs appeared almost exclusively fragmented or 'nicked' into tRNA-derived small RNAs (tsRNAs) with lengths between 18 to 35 nucleotides. We propose that in mammalian cells, tRNAs that lack a functional 3' CCA tail are selectively sorted into EVs and shuttled out of the producing cell, offering a new perspective into the physiological role of secreted EVs and luminal cargo-selection.

The tRNAVal half: A strong endogenous Toll-like receptor 7 ligand with a 5'-terminal universal sequence signature

Toll-like receptors (TLRs) are crucial components of the innate immune system. Endosomal TLR7 recognizes single-stranded RNAs, yet its endogenous ssRNA ligands are not fully understood. We previously showed that extracellular (ex-) 5'-half molecules of tRNAHisGUG (the 5'-tRNAHisGUG half) in extracellular vesicles (EVs) of human macrophages activate TLR7 when delivered into endosomes of recipient macrophages. Here, we fully explored immunostimulatory ex-5'-tRNA half molecules and identified the 5'-tRNAValCAC/AAC half, the most abundant tRNA-derived RNA in macrophage EVs, as another 5'-tRNA half molecule with strong TLR7 activation capacity. Levels of the ex-5'-tRNAValCAC/AAC half were highly up-regulated in macrophage EVs upon exposure to lipopolysaccharide and in the plasma of patients infected with Mycobacterium tuberculosis. The 5'-tRNAValCAC/AAC half-mediated activation of TLR7 effectively eradicated bacteria infected in macrophages. Mutation analyses of the 5'-tRNAValCAC/AAC half identified the terminal GUUU sequence as a determinant for TLR7 activation. We confirmed that GUUU is the optimal ratio of guanosine and uridine for TLR7 activation; microRNAs or other RNAs with the terminal GUUU motif can indeed stimulate TLR7, establishing the motif as a universal signature for TLR7 activation. These results advance our understanding of endogenous ssRNA ligands of TLR7 and offer insights into diverse TLR7-involved pathologies and their therapeutic strategies.

Untacking small RNA profiling and RNA fragment footprinting: Approaches and challenges in library construction

Small RNAs (sRNAs) with sizes ranging from 15 to 50 nucleotides (nt) are critical regulators of gene expression control. Prior studies have shown that sRNAs are involved in a broad range of biological processes, such as organ development, tumorigenesis, and epigenomic regulation; however, emerging evidence unveils a hidden layer of diversity and complexity of endogenously encoded sRNAs profile in eukaryotic organisms, including novel types of sRNAs and the previously unknown post-transcriptional RNA modifications. This underscores the importance for accurate, unbiased detection of sRNAs in various cellular contexts. A multitude of high-throughput methods based on next-generation sequencing (NGS) are developed to decipher the sRNA expression and their modifications. Nonetheless, distinct from mRNA sequencing, the data from sRNA sequencing suffer frequent inconsistencies and high variations emanating from the adapter contaminations and RNA modifications, which overall skew the sRNA libraries. Here, we summarize the sRNA-sequencing approaches, and discuss the considerations and challenges for the strategies and methods of sRNA library construction. The pros and cons of sRNA sequencing have significant implications for implementing RNA fragment footprinting approaches, including CLIP-seq and Ribo-seq. We envision that this review can inspire novel improvements in small RNA sequencing and RNA fragment footprinting in future. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Processing of Small RNAs Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

Biochemical and structural insights into a 5' to 3' RNA ligase reveal a potential role in tRNA ligation

ATP-grasp superfamily enzymes contain a hand-like ATP-binding fold and catalyze a variety of reactions using a similar catalytic mechanism. More than 30 protein families are categorized in this superfamily, and they are involved in a plethora of cellular processes and human diseases. Here we identify C12orf29 as an atypical ATP-grasp enzyme that ligates RNA. Human C12orf29 and its homologs auto-adenylate on an active site Lys residue as part of a reaction intermediate that specifically ligates RNA halves containing a 5'-phosphate and a 3'-hydroxyl. C12orf29 binds tRNA in cells and can ligate tRNA within the anticodon loop in vitro. Genetic depletion of c12orf29 in female mice alters global tRNA levels in brain. Furthermore, crystal structures of a C12orf29 homolog from Yasminevirus bound to nucleotides reveal a minimal and atypical RNA ligase fold with a unique active site architecture that participates in catalysis. Collectively, our results identify C12orf29 as an RNA ligase and suggest its involvement in tRNA biology.

Mechanisms used for cDNA synthesis and site-specific integration of RNA into DNA genomes by a reverse transcriptase-Cas1 fusion protein

Reverse transcriptase-Cas1 (RT-Cas1) fusion proteins found in some CRISPR systems enable spacer acquisition from both RNA and DNA, but the mechanism of RNA spacer acquisition has remained unclear. Here, we found that Marinomonas mediterranea RT-Cas1/Cas2 adds short 3'-DNA (dN) tails to RNA protospacers, enabling their direct integration into CRISPR arrays as 3'-dN-RNAs or 3'-dN-RNA/cDNA duplexes at rates comparable to similarly configured DNAs. Reverse transcription of RNA protospacers is initiated at 3' proximal sites by multiple mechanisms, including recently described de novo initiation, protein priming with any dNTP, and use of short exogenous or synthesized DNA oligomer primers, enabling synthesis of near full-length cDNAs of diverse RNAs without fixed sequence requirements. The integration of 3'-dN-RNAs or single-stranded DNAs (ssDNAs) is favored over duplexes at higher protospacer concentrations, potentially relevant to spacer acquisition from abundant pathogen RNAs or ssDNA fragments generated by phage defense nucleases. Our findings reveal mechanisms for site-specifically integrating RNA into DNA genomes with potential biotechnological applications.

Small RNA structural biochemistry in a post-sequencing era

High-throughput sequencing has had an enormous impact on small RNA research during the past decade. However, sequencing only offers a one-dimensional view of the transcriptome and is often highly biased. Additionally, the 'sequence, map and annotate' approach, used widely in small RNA research, can lead to flawed interpretations of the data, lacking biological plausibility, due in part to database issues. Even in the absence of technical biases, the loss of three-dimensional information is a major limitation to understanding RNA stability, turnover and function. For example, noncoding RNA-derived fragments seem to exist mainly as dimers, tetramers or as nicked forms of their parental RNAs, contrary to widespread assumptions. In this perspective, we will discuss main sources of bias during small RNA-sequencing, present several useful bias-reducing strategies and provide guidance on the interpretation of small RNA-sequencing results, with emphasis on RNA fragmentomics. As sequencing offers a one-dimensional projection of a four-dimensional reality, prior structure-level knowledge is often needed to make sense of the data. Consequently, while less-biased sequencing methods are welcomed, integration of orthologous experimental techniques is also strongly recommended.

Precision run-on sequencing (PRO-seq) for microbiome transcriptomics

Bacteria respond to environmental stimuli through precise regulation of transcription initiation and elongation. Bulk RNA sequencing primarily characterizes mature transcripts, so to identify actively transcribed loci we need to capture RNA polymerase (RNAP) complexed with nascent RNA. However, such capture methods have only previously been applied to culturable, genetically tractable organisms such as E. coli and B. subtilis. Here we apply precision run-on sequencing (PRO-seq) to profile nascent transcription in cultured E. coli and diverse uncultured bacteria. We demonstrate that PRO-seq can characterize the transcription of small, structured, or post-transcriptionally modified RNAs, which are often absent from bulk RNA-seq libraries. Applying PRO-seq to the human microbiome highlights taxon-specific RNAP pause motifs and pause-site distributions across non-coding RNA loci that reflect structure-coincident pausing. We also uncover concurrent transcription and cleavage of CRISPR guide RNAs and transfer RNAs. We demonstrate the utility of PRO-seq for exploring transcriptional dynamics in diverse microbial communities.

Emerging functional principles of tRNA-derived small RNAs and other regulatory small RNAs

Recent advancements in small RNA sequencing have unveiled a previously hidden world of regulatory small noncoding RNAs (sncRNAs) that extend beyond the well-studied small interfering RNAs, microRNAs, and piwi-interacting RNAs. This exploration, starting with tRNA-derived small RNAs, has led to the discovery of a diverse universe of sncRNAs derived from various longer structured RNAs such as rRNAs, small nucleolar RNAs, small nuclear RNAs, Y RNAs, and vault RNAs, with exciting uncharted functional possibilities. In this perspective, we discuss the emerging functional principles of sncRNAs beyond the well-known RNAi-like mechanisms, focusing on those that operate independent of linear sequence complementarity but rather function in an aptamer-like fashion. Aptamers use 3D structure for specific interactions with ligands and are modulated by RNA modifications and subcellular environments. Given that aptamer-like sncRNA functions are widespread and present in species lacking RNAi, they may represent an ancient functional principle that predates RNAi. We propose a rethinking of the origin of RNAi and its relationship with these aptamer-like functions in sncRNAs and how these complementary mechanisms shape biological processes. Lastly, the aptamer-like function of sncRNAs highlights the need for caution in using small RNA mimics in research and therapeutics, as their specificity is not restricted solely to linear sequence.

Mechanisms used for cDNA synthesis and site-specific integration of RNA into DNA genomes by a reverse transcriptase-Cas1 fusion protein

Reverse transcriptase-Cas1 (RT-Cas1) fusion proteins found in some CRISPR systems enable spacer acquisition from both RNA and DNA, but the mechanism of RNA spacer acquisition has remained unclear. Here, we found Marinomonas mediterranea RT-Cas1/Cas2 adds short 3'-DNA (dN) tails to RNA protospacers enabling their direct integration into CRISPR arrays as 3'-dN-RNA/cDNA duplexes or 3'-dN-RNAs at rates comparable to similarly configured DNAs. Reverse transcription of RNA protospacers occurs by multiple mechanisms, including recently described de novo initiation, protein priming with any dNTP, and use of short exogenous or synthesized DNA oligomer primers, enabling synthesis of cDNAs from diverse RNAs without fixed sequence requirements. The integration of 3'-dN-RNAs or single-stranded (ss) DNAs is favored over duplexes at higher protospacer concentrations, potentially relevant to spacer acquisition from abundant pathogen RNAs or ssDNA fragments generated by phage-defense nucleases. Our findings reveal novel mechanisms for site-specifically integrating RNA into DNA genomes with potential biotechnological applications.

Intronic small nucleolar RNAs regulate host gene splicing through base pairing with their adjacent intronic sequences

BACKGROUND

Small nucleolar RNAs (snoRNAs) are abundant noncoding RNAs best known for their involvement in ribosomal RNA maturation. In mammals, most expressed snoRNAs are embedded in introns of longer genes and produced through transcription and splicing of their host. Intronic snoRNAs were long viewed as inert passengers with little effect on host expression. However, a recent study reported a snoRNA influencing the splicing and ultimate output of its host gene. Overall, the general contribution of intronic snoRNAs to host expression remains unclear.

RESULTS

Computational analysis of large-scale human RNA-RNA interaction datasets indicates that 30% of detected snoRNAs interact with their host transcripts. Many snoRNA-host duplexes are located near alternatively spliced exons and display high sequence conservation suggesting a possible role in splicing regulation. The study of the model SNORD2-EIF4A2 duplex indicates that the snoRNA interaction with the host intronic sequence conceals the branch point leading to decreased inclusion of the adjacent alternative exon. Extended SNORD2 sequence containing the interacting intronic region accumulates in sequencing datasets in a cell-type-specific manner. Antisense oligonucleotides and mutations that disrupt the formation of the snoRNA-intron structure promote the splicing of the alternative exon, shifting the EIF4A2 transcript ratio away from nonsense-mediated decay.

CONCLUSIONS

Many snoRNAs form RNA duplexes near alternative exons of their host transcripts, placing them in optimal positions to control host output as shown for the SNORD2-EIF4A2 model system. Overall, our study supports a more widespread role for intronic snoRNAs in the regulation of their host transcript maturation.

Artifacts and biases of the reverse transcription reaction in RNA sequencing

RNA sequencing has spurred a significant number of research areas in recent years. Most protocols rely on synthesizing a more stable complementary DNA (cDNA) copy of the RNA molecule during the reverse transcription reaction. The resulting cDNA pool is often wrongfully assumed to be quantitatively and molecularly similar to the original RNA input. Sadly, biases and artifacts confound the resulting cDNA mixture. These issues are often overlooked or ignored in the literature by those that rely on the reverse transcription process. In this review, we confront the reader with intra- and intersample biases and artifacts caused by the reverse transcription reaction during RNA sequencing experiments. To fight the reader's despair, we also provide solutions to most issues and inform on good RNA sequencing practices. We hope the reader can use this review to their advantage, thereby contributing to scientifically sound RNA studies.

A ligation-independent sequencing method reveals tRNA-derived RNAs with blocked 3' termini

Despite the numerous sequencing methods available, the vast diversity in size and chemical modifications of RNA molecules makes the capture of the full spectrum of cellular RNAs a difficult task. By combining quasi-random hexamer priming with a custom template switching strategy, we developed a method to construct sequencing libraries from RNA molecules of any length and with any type of 3' terminal modification, allowing the sequencing and analysis of virtually all RNA species. Ligation-independent detection of all types of RNA (LIDAR) is a simple, effective tool to comprehensively characterize changes in small non-coding RNAs and mRNAs simultaneously, with performance comparable to separate dedicated methods. With LIDAR, we comprehensively characterized the coding and non-coding transcriptome of mouse embryonic stem cells, neural progenitor cells, and sperm. LIDAR detected a much larger variety of tRNA-derived RNAs (tDRs) compared to traditional ligation-dependent sequencing methods, and uncovered the presence of tDRs with blocked 3' ends that had previously escaped detection. Our findings highlight the potential of LIDAR to systematically detect all RNAs in a sample and uncover new RNA species with potential regulatory functions.

Improved Bacterial Single-Cell RNA-Seq through Automated MATQ-Seq and Cas9-Based Removal of rRNA Reads

Bulk RNA sequencing technologies have provided invaluable insights into host and bacterial gene expression and associated regulatory networks. Nevertheless, the majority of these approaches report average expression across cell populations, hiding the true underlying expression patterns that are often heterogeneous in nature. Due to technical advances, single-cell transcriptomics in bacteria has recently become reality, allowing exploration of these heterogeneous populations, which are often the result of environmental changes and stressors. In this work, we have improved our previously published bacterial single-cell RNA sequencing (scRNA-seq) protocol that is based on multiple annealing and deoxycytidine (dC) tailing-based quantitative scRNA-seq (MATQ-seq), achieving a higher throughput through the integration of automation. We also selected a more efficient reverse transcriptase, which led to reduced cell loss and higher workflow robustness. Moreover, we successfully implemented a Cas9-based rRNA depletion protocol into the MATQ-seq workflow. Applying our improved protocol on a large set of single Salmonella cells sampled over different growth conditions revealed improved gene coverage and a higher gene detection limit compared to our original protocol and allowed us to detect the expression of small regulatory RNAs, such as GcvB or CsrB at a single-cell level. In addition, we confirmed previously described phenotypic heterogeneity in Salmonella in regard to expression of pathogenicity-associated genes. Overall, the low percentage of cell loss and high gene detection limit makes the improved MATQ-seq protocol particularly well suited for studies with limited input material, such as analysis of small bacterial populations in host niches or intracellular bacteria. IMPORTANCE Gene expression heterogeneity among isogenic bacteria is linked to clinically relevant scenarios, like biofilm formation and antibiotic tolerance. The recent development of bacterial single-cell RNA sequencing (scRNA-seq) enables the study of cell-to-cell variability in bacterial populations and the mechanisms underlying these phenomena. Here, we report a scRNA-seq workflow based on MATQ-seq with increased robustness, reduced cell loss, and improved transcript capture rate and gene coverage. Use of a more efficient reverse transcriptase and the integration of an rRNA depletion step, which can be adapted to other bacterial single-cell workflows, was instrumental for these improvements. Applying the protocol to the foodborne pathogen Salmonella, we confirmed transcriptional heterogeneity across and within different growth phases and demonstrated that our workflow captures small regulatory RNAs at a single-cell level. Due to low cell loss and high transcript capture rates, this protocol is uniquely suited for experimental settings in which the starting material is limited, such as infected tissues.

Arg-tRNA synthetase links inflammatory metabolism to RNA splicing and nuclear trafficking via SRRM2

Cells respond to perturbations such as inflammation by sensing changes in metabolite levels. Especially prominent is arginine, which has known connections to the inflammatory response. Aminoacyl-tRNA synthetases, enzymes that catalyse the first step of protein synthesis, can also mediate cell signalling. Here we show that depletion of arginine during inflammation decreased levels of nuclear-localized arginyl-tRNA synthetase (ArgRS). Surprisingly, we found that nuclear ArgRS interacts and co-localizes with serine/arginine repetitive matrix protein 2 (SRRM2), a spliceosomal and nuclear speckle protein, and that decreased levels of nuclear ArgRS correlated with changes in condensate-like nuclear trafficking of SRRM2 and splice-site usage in certain genes. These splice-site usage changes cumulated in the synthesis of different protein isoforms that altered cellular metabolism and peptide presentation to immune cells. Our findings uncover a mechanism whereby an aminoacyl-tRNA synthetase cognate to a key amino acid that is metabolically controlled during inflammation modulates the splicing machinery.

A New EBS2b-IBS2b Base Paring (A-8/T-8) Improved the Gene-Targeting Efficiency of Thermotargetron in Escherichia coli

Thermophilic group II intron is one type of retrotransposon composed of intron RNA and intron-encoded protein (IEP), which can be utilized in gene targeting by harnessing their novel ribozyme-based DNA integration mechanism termed "retrohoming." It is mediated by a ribonucleoprotein (RNP) complex that contains the excised intron lariat RNA and an IEP with reverse transcriptase (RT) activity. The RNP recognizes targeting sites by exon-binding sequences 2 (EBS2)/intron-binding sequences 2 (IBS2), EBS1/IBS1, and EBS3/IBS3 bases pairing. Previously, we developed the TeI3c/4c intron as a thermophilic gene targeting system-Thermotargetron (TMT). However, we found that the targeting efficiency of TMT varies significantly at different targeting sites, which leads to a relatively low success rate. To further improve the success rate and gene-targeting efficiency of TMT, we constructed a Random Gene-targeting Plasmids Pool (RGPP) to analyze the sequence recognition preference of TMT. A new base pairing, located at the -8 site between EBS2/IBS2 and EBS1/IBS1 (named EBS2b-IBS2b), increased the success rate (2.45- to 5.07-fold) and significantly improved gene-targeting efficiency of TMT. A computer algorithm (TMT 1.0), based on the newly discovered sequence recognition roles, was also developed to facilitate the design of TMT gene-targeting primers. The present work could essentially expand the practicalities of TMT in the genome engineering of heat-tolerance mesophilic and thermophilic bacteria. IMPORTANCE The randomized base pairing in the interval of IBS2 and IBS1 of Tel3c/4c intron (-8 and -7 sites) in Thermotargetron (TMT) results in a low success rate and gene-targeting efficiency in bacteria. In the present work, we constructed a randomized gene-targeting plasmids pool (RGPP) to study whether there is a base preference in target sequences. Among all the successful "retrohoming" targets, we found that a new EBS2b-IBS2b base paring (A-8/T-8) significantly increased TMT's gene-targeting efficiency, and the concept is also applicable to other gene targets in redesigned gene-targeting plasmids pool in E. coli. The improved TMT is a promising tool for the genetic engineering of bacteria and could promote metabolic engineering and synthetic biology research in valuable microbes that recalcitrance for genetic manipulation.

Pitfalls of using sequence databases for heterologous expression studies - a technical review

Synthesis of DNA fragments based on gene sequences that are available in public resources has become an efficient and affordable method that has gradually replaced traditional cloning efforts such as PCR cloning from cDNA. However, database entries based on genome sequencing results are prone to errors which can lead to false sequence information and, ultimately, errors in functional characterisation of proteins such as ion channels and transporters in heterologous expression systems. We have identified five common problems that repeatedly appear in public resources: (1) Not every gene has yet been annotated; (2) not all gene annotations are necessarily correct; (3) transcripts may contain automated corrections; (4) there are mismatches between gene, mRNA and protein sequences; and (5) splicing patterns often lack experimental validation. This technical review highlights and provides a strategy to bypass these issues in order to avoid critical mistakes that could impact future studies of any gene/protein of interest in heterologous expression systems.

Nano3P-seq: transcriptome-wide analysis of gene expression and tail dynamics using end-capture nanopore cDNA sequencing

RNA polyadenylation plays a central role in RNA maturation, fate, and stability. In response to developmental cues, polyA tail lengths can vary, affecting the translation efficiency and stability of mRNAs. Here we develop Nanopore 3' end-capture sequencing (Nano3P-seq), a method that relies on nanopore cDNA sequencing to simultaneously quantify RNA abundance, tail composition, and tail length dynamics at per-read resolution. By employing a template-switching-based sequencing protocol, Nano3P-seq can sequence RNA molecule from its 3' end, regardless of its polyadenylation status, without the need for PCR amplification or ligation of RNA adapters. We demonstrate that Nano3P-seq provides quantitative estimates of RNA abundance and tail lengths, and captures a wide diversity of RNA biotypes. We find that, in addition to mRNA and long non-coding RNA, polyA tails can be identified in 16S mitochondrial ribosomal RNA in both mouse and zebrafish models. Moreover, we show that mRNA tail lengths are dynamically regulated during vertebrate embryogenesis at an isoform-specific level, correlating with mRNA decay. Finally, we demonstrate the ability of Nano3P-seq in capturing non-A bases within polyA tails of various lengths, and reveal their distribution during vertebrate embryogenesis. Overall, Nano3P-seq is a simple and robust method for accurately estimating transcript levels, tail lengths, and tail composition heterogeneity in individual reads, with minimal library preparation biases, both in the coding and non-coding transcriptome.

MiRNAs and snoRNAs in Bone Metastasis: Functional Roles and Clinical Potential

Bone is a frequent site of metastasis. Bone metastasis is associated with a short-term prognosis in cancer patients, and current treatments aim to slow its growth, but are rarely curative. Thus, revealing molecular mechanisms that explain why metastatic cells are attracted to the bone micro-environment, and how they successfully settle in the bone marrow-taking advantage over bone resident cells-and grow into macro-metastasis, is essential to propose new therapeutic approaches. MicroRNAs and snoRNAs are two classes of small non-coding RNAs that post-transcriptionally regulate gene expression. Recently, microRNAs and snoRNAs have been pointed out as important players in bone metastasis by (i) preparing the pre-metastatic niche, directly and indirectly affecting the activities of osteoclasts and osteoblasts, (ii) promoting metastatic properties within cancer cells, and (iii) acting as mediators within cells to support cancer cell growth in bone. This review aims to highlight the importance of microRNAs and snoRNAs in metastasis, specifically in bone, and how their roles can be linked together. We then discuss how microRNAs and snoRNAs are secreted by cancer cells and be found as extracellular vesicle cargo. Finally, we provide evidence of how microRNAs and snoRNAs can be potential therapeutic targets, at least in pre-clinical settings, and how their detection in liquid biopsies can be a useful diagnostic and/or prognostic biomarker to predict the risk of relapse in cancer patients.

RNA modifications: an overview of select web-based tools

The field of epitranscriptomics has expanded dramatically in recent years, both in the number of identified RNA modifications and the number of researchers studying them. As knowledge of post-transcriptional modifications continues to expand, numerous new methods have been developed to detect these modifications. Additionally, modifications are being extended to therapeutic settings, such as with recent mRNA vaccines. With this increase in knowledge and use, the community is recognizing the necessity for user-friendly databases to (i) store information from both high- and low-throughput studies and (ii) provide prediction software on how RNA modifications contribute to RNA function and disease. This mini-review highlights select RNA modification databases and their key attributes with the aim of providing a resource to researchers in the field of epitranscriptomics.

i-tRAP (individual tRNA acylation PCR): A convenient method for selective quantification of tRNA charging

Each transfer RNA (tRNA) is aminoacylated (charged) with a genetic codon-specific amino acid at its 3' end. Charged tRNAs are primarily used for translation, whereas fluctuations in charged tRNA fractions are known to reflect cellular response to stress. Here we report the development of individual tRNA-acylation using PCR (i-tRAP), a convenient PCR-based method that can specifically quantify individual tRNA charging ratio. In this i-tRAP method, demethylases remove base methylations which are problematic for reverse transcription reaction, and β-elimination reaction specifically removes the 3' end of adenine residue in uncharged tRNA. Subsequent TaqMan MGB qRT-PCR can distinguish between cDNA of charged tRNA and uncharged tRNA. By using this method, we revealed that the charging ratio of tRNAGln(CUG) was changed in response to amino acid starvation and also the charging ratio of tRNAGln(CUG) in senescent cells was lower than in young cells under starvation conditions. i-tRAP can be applicable to the quantification of charging ratio of various tRNAs, and provides a simple and convenient method for analyzing tRNA charging.

Group II intron-like reverse transcriptases function in double-strand break repair

Bacteria encode reverse transcriptases (RTs) of unknown function that are closely related to group II intron-encoded RTs. We found that a Pseudomonas aeruginosa group II intron-like RT (G2L4 RT) with YIDD instead of YADD at its active site functions in DNA repair in its native host and when expressed in Escherichia coli. G2L4 RT has biochemical activities strikingly similar to those of human DNA repair polymerase θ and uses them for translesion DNA synthesis and double-strand break repair (DSBR) via microhomology-mediated end-joining (MMEJ). We also found that a group II intron RT can function similarly in DNA repair, with reciprocal active-site substitutions showing isoleucine favors MMEJ and alanine favors primer extension in both enzymes. These DNA repair functions utilize conserved structural features of non-LTR-retroelement RTs, including human LINE-1 and other eukaryotic non-LTR-retrotransposon RTs, suggesting such enzymes may have inherent ability to function in DSBR in a wide range of organisms.

Emerging drug targets for triple-negative breast cancer: a guided tour of the preclinical landscape

INTRODUCTION

Triple-negative breast cancer (TNBC) is the most fatal molecular subtype of breast cancer because of its aggressiveness and resistance to chemotherapy. FDA-approved therapies for TNBC are limited to poly(ADP-ribose) polymerase inhibitors, immune checkpoint inhibitors, and trophoblast cell surface antigen 2-targeted antibody-drug conjugate. Therefore, developing a novel effective targeted therapy for TNBC is an urgent unmet need.

AREAS COVERED

In this narrative review, we discuss emerging targets for TNBC treatment discovered in early translational studies. We focus on cancer cell membrane molecules, hyperactive intracellular signaling pathways, and the tumor microenvironment (TME) based on their druggability, therapeutic potency, specificity to TNBC, and application in immunotherapy.

EXPERT OPINION

The significant challenges in the identification and validation of TNBC-associated targets are 1) application of appropriate genetic, molecular, and immunological approaches for modulating the target, 2) establishment of a proper mouse model that accurately represents the human immune TME, 3) TNBC molecular heterogeneity, and 4) failure translation of preclinical findings to clinical practice. To overcome those difficulties, future research needs to apply novel technology, such as single-cell RNA sequencing, thermostable group II intron reverse transcriptase sequencing, and humanized mouse models. Further, combination treatment targeting multiple pathways in both the TNBC tumor and its TME is essential for effective disease control.

TRMT6/61A-dependent base methylation of tRNA-derived fragments regulates gene-silencing activity and the unfolded protein response in bladder cancer

RNA modifications are important regulatory elements of RNA functions. However, most genome-wide mapping of RNA modifications has focused on messenger RNAs and transfer RNAs, but such datasets have been lacking for small RNAs. Here we mapped N1-methyladenosine (m1A) in the cellular small RNA space. Benchmarked with synthetic m1A RNAs, our workflow identified specific groups of m1A-containing small RNAs, which are otherwise disproportionally under-represented. In particular, 22-nucleotides long 3' tRNA-fragments are highly enriched for TRMT6/61A-dependent m1A located within the seed region. TRMT6/61A-dependent m1A negatively affects gene silencing by tRF-3s. In urothelial carcinoma of the bladder, where TRMT6/61A is over-expressed, higher m1A modification on tRFs is detected, correlated with a dysregulation of tRF targetome. Lastly, TRMT6/61A regulates tRF-3 targets involved in unfolded protein response. Together, our results reveal a mechanism of regulating gene expression via base modification of small RNA.

Exploring the expanding universe of small RNAs

The world of small noncoding RNAs (sncRNAs) is ever-expanding, from small interfering RNA, microRNA and Piwi-interacting RNA to the recently emerging non-canonical sncRNAs derived from longer structured RNAs (for example, transfer, ribosomal, Y, small nucleolar, small nuclear and vault RNAs), showing distinct biogenesis and functional principles. Here we discuss recent tools for sncRNA identification, caveats in sncRNA expression analysis and emerging methods for direct sequencing of sncRNAs and systematic mapping of RNA modifications that are integral to their function.

A Novel Tissue Atlas and Online Tool for the Interrogation of Small RNA Expression in Human Tissues and Biofluids

One promising goal for utilizing the molecular information circulating in biofluids is the discovery of clinically useful biomarkers. Extracellular RNAs (exRNAs) are one of the most diverse classes of molecular cargo, easily assayed by sequencing and with expressions that rapidly change in response to subject status. Despite diverse exRNA cargo, most evaluations from biofluids have focused on small RNA sequencing and analysis, specifically on microRNAs (miRNAs). Another goal of characterizing circulating molecular information, is to correlate expression to injuries associated with specific tissues of origin. Biomarker candidates are often described as being specific, enriched in a particular tissue or associated with a disease process. Likewise, miRNA data is often reported to be specific, enriched for a tissue, without rigorous testing to support the claim. Here we provide a tissue atlas of small RNAs from 30 different tissues and three different blood cell types. We analyzed the tissues for enrichment of small RNA sequences and assessed their expression in biofluids: plasma, cerebrospinal fluid, urine, and saliva. We employed published data sets representing physiological (resting vs. acute exercise) and pathologic states (early- vs. late-stage liver fibrosis, and differential subtypes of stroke) to determine differential tissue-enriched small RNAs. We also developed an online tool that provides information about exRNA sequences found in different biofluids and tissues. The data can be used to better understand the various types of small RNA sequences in different tissues as well as their potential release into biofluids, which should help in the validation or design of biomarker studies.

Role of Fasciola hepatica Small RNAs in the Interaction With the Mammalian Host

MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression being involved in many different biological processes and play a key role in developmental timing. Additionally, recent studies have shown that miRNAs released from parasites are capable of regulating the expression of host genes. In the present work, we studied the expression patterns of ncRNAs of various intra-mammalian life-cycle stages of the liver fluke, Fasciola hepatica, as well as those packaged into extracellular vesicles and shed by the adult fluke. The miRNA expression profile of the intra-mammalian stages shows important variations, despite a set of predominant miRNAs that are highly expressed across all stages. No substantial variations in miRNA expression between dormant and activated metacercariae were detected, suggesting that they might not be central players in regulating fluke gene expression during this crucial step in the invasion of the definitive host. We generated a curated pipeline for the prediction of putative target genes that reports only sites conserved between three different prediction approaches. This pipeline was tested against an iso-seq curated database of the 3’ UTR regions of F. hepatica genes to detect miRNA regulation networks within liver fluke. Several functions related to the host immune response or modulation were enriched among the targets of the most highly expressed parasite miRNAs, stressing that they might be key players during the establishment and maintenance of infection. Additionally, we detected fragments derived from the processing of tRNAs, in all developmental stages analyzed, and documented the presence of novel long tRNA fragments enriched in vesicles. We confirmed the presence of at least 5 putative vault RNAs (vtRNAs), that are expressed across different stages and enriched in vesicles. The presence of tRNA fragments and vtRNAs in vesicles raise the possibility that they could be involved in the host-parasite interaction.

A method for simultaneous detection of small and long RNA biotypes by ribodepleted RNA-Seq

RNA sequencing offers unprecedented access to the transcriptome. Key to this is the identification and quantification of many different species of RNA from the same sample at the same time. In this study we describe a novel protocol for simultaneous detection of coding and non-coding transcripts using modifications to the Ion Total RNA-Seq kit v2 protocol, with integration of QIASeq FastSelect rRNA removal kit. We report highly consistent sequencing libraries can be produced from both frozen high integrity mouse hippocampal tissue and the more challenging post-mortem human tissue. Removal of rRNA using FastSelect was extremely efficient, resulting in less than 1.5% rRNA content in the final library. We identified > 30,000 unique transcripts from all samples, including protein-coding genes and many species of non-coding RNA, in biologically-relevant proportions. Furthermore, the normalized sequencing read count for select genes significantly negatively correlated with Ct values from qRT-PCR analysis from the same samples. These results indicate that this protocol accurately and consistently identifies and quantifies a wide variety of transcripts simultaneously. The highly efficient rRNA depletion, coupled with minimized sample handling and without complicated and high-loss size selection protocols, makes this protocol useful to researchers wishing to investigate whole transcriptomes.

TGIRT-seq Protocol for the Comprehensive Profiling of Coding and Non-coding RNA Biotypes in Cellular, Extracellular Vesicle, and Plasma RNAs

High-throughput RNA sequencing (RNA-seq) has extraordinarily advanced our understanding of gene expression and disease etiology, and is a powerful tool for the identification of biomarkers in a wide range of organisms. However, most RNA-seq methods rely on retroviral reverse transcriptases (RTs), enzymes that have inherently low fidelity and processivity, to convert RNAs into cDNAs for sequencing. Here, we describe an RNA-seq protocol using Thermostable Group II Intron Reverse Transcriptases (TGIRTs), which have high fidelity, processivity, and strand-displacement activity, as well as a proficient template-switching activity that enables efficient and seamless RNA-seq adapter addition. By combining these activities, TGIRT-seq enables the simultaneous profiling of all RNA biotypes from small amounts of starting material, with superior RNA-seq metrics, and unprecedented ability to sequence structured RNAs. The TGIRT-seq protocol for Illumina sequencing consists of three steps: (i) addition of a 3' RNA-seq adapter, coupled to the initiation of cDNA synthesis at the 3' end of a target RNA, via template switching from a synthetic adapter RNA/DNA starter duplex; (ii) addition of a 5' RNA-seq adapter, by using thermostable 5' App DNA/RNA ligase to ligate an adapter oligonucleotide to the 3' end of the completed cDNA; (iii) minimal PCR amplification, to add capture sites and indices for Illumina sequencing. TGIRT-seq for the Illumina sequencing platform has been used for comprehensive profiling of coding and non-coding RNAs in ribodepleted, chemically fragmented cellular RNAs, and for the analysis of intact (non-chemically fragmented) cellular, extracellular vesicle (EV), and plasma RNAs, where it yields continuous full-length end-to-end sequences of structured small non-coding RNAs (sncRNAs), including tRNAs, snoRNAs, snRNAs, pre-miRNAs, and full-length excised linear intron (FLEXI) RNAs. Graphic abstract: Figure 1.Overview of the TGIRT-seq protocol for Illumina sequencing.Major steps are: (1) Template switching from a synthetic R2 RNA/R2R DNA starter duplex with a 1-nt 3' DNA overhang (a mixture of A, C, G, and T residues, denoted N) that base pairs to the 3' nucleotide of a target RNA, and upon initiating reverse transcription by adding dNTPs, seamlessly links an R2R adapter to the 5' end of the resulting cDNA; (2) Ligation of an R1R adapter to the 3' end of the completed cDNA; and (3) Minimal PCR amplification with primers that add Illumina capture sites (P5 and P7) and barcode sequences (indices 5 and 7). The index 7 barcode is required, while the index 5 barcode is optional, to provide unique dual indices (UDIs).

Cell type-specific analysis by single-cell profiling identifies a stable mammalian tRNA-mRNA interface and increased translation efficiency in neurons

The correlation between codon and anticodon pools influences the efficiency of translation, but whether differences exist in these pools across individual cells is unknown. We determined that codon usage and amino acid demand are highly stable across different cell types using available mouse and human single-cell RNA-sequencing atlases. After showing the robustness of ATAC-sequencing measurements for the analysis of tRNA gene usage, we quantified anticodon usage and amino acid supply in both mouse and human single-cell ATAC-seq atlases. We found that tRNA gene usage is overall coordinated across cell types, except in neurons, which clustered separately from other cell types. Integration of these data sets revealed a strong and statistically significant correlation between amino acid supply and demand across almost all cell types. Neurons have an enhanced translation efficiency over other cell types, driven by an increased supply of tRNA^Ala (AGC) anticodons. This results in faster decoding of the Ala-GCC codon, as determined by cell type–specific ribosome profiling, suggesting that the reduction of tRNA^Ala (AGC) anticodon pools may be implicated in neurological pathologies. This study, the first such examination of codon usage, anticodon usage, and translation efficiency resolved at the cell-type level with single-cell information, identifies a conserved landscape of translation elongation across mammalian cellular diversity and evolution.

5´XP sRNA-seq: efficient identification of transcripts with and without 5´ phosphorylation reveals evolutionary conserved small RNA

Small RNA (sRNA) sequencing has been critical for our understanding of many cellular processes, including gene regulation. Nonetheless, the varying biochemical properties of sRNA, such as 5´ nucleotide modifications, make many sRNA subspecies incompatible with common protocols for sRNA sequencing. Here we describe 5XP-seq that outlines a novel strategy that captures a more complete picture of sRNA. By tagging 5´P sRNA during library preparation, 5XP-seq combines an open approach that includes all types of 5'-terminal modifications (5´X), with a selective approach for 5-phosphorylated sRNA (5´P). We show that 5XP-seq not only enriches phosphorylated miRNA and piRNA but successfully discriminates these sRNA from all other sRNA species. We further demonstrate the importance of this strategy by successful inter-species validation of sRNAs that would have otherwise failed, including human to insect translation of several tRNA (tRFs) and rRNA (rRFs) fragments. By combining 5´ insensitive library strategies with 5´ sensitive tagging, we have successfully tackled an intrinsic bias in modern sRNA sequencing that will help us reveal the true complexity and the evolutionary significance of the sRNA world.

Low-bias ncRNA libraries using ordered two-template relay: Serial template jumping by a modified retroelement reverse transcriptase

Retrotransposons are noninfectious, mobile genetic elements that proliferate in host genomes via an RNA intermediate that is copied into DNA by a reverse transcriptase (RT) enzyme. RTs are important for biotechnological applications involving information capture from RNA since RNA is first converted into complementary DNA for detection or sequencing. Here, we biochemically characterized RTs from two retroelements and uncovered several activities that allowed us to design a streamlined, efficient workflow for determining the inventory of RNA sequences in processed RNA pools. The unique properties of nonretroviral RT activities obviate many technical issues associated with current methods of RNA sequence analysis, with wide applications in research, biotechnology, and diagnostics.

Selfish, non-long terminal repeat (non-LTR) retroelements and mobile group II introns encode reverse transcriptases (RTs) that can initiate DNA synthesis without substantial base pairing of primer and template. Biochemical characterization of these enzymes has been limited by recombinant expression challenges, hampering understanding of their properties and the possible exploitation of their properties for research and biotechnology. We investigated the activities of representative RTs using a modified non-LTR RT from Bombyx mori and a group II intron RT from Eubacterium rectale. Only the non-LTR RT supported robust and serial template jumping, producing one complementary DNA (cDNA) from several templates each copied end to end. We also discovered an unexpected terminal deoxynucleotidyl transferase activity of the RTs that adds nucleotide(s) of choice to 3′ ends of single- and/or double-stranded RNA or DNA. Combining these two types of activity with additional insights about nontemplated nucleotide additions to duplexed cDNA product, we developed a streamlined protocol for fusion of next-generation sequencing adaptors to both cDNA ends in a single RT reaction. When benchmarked using a reference pool of microRNAs (miRNAs), library production by Ordered Two-Template Relay (OTTR) using recombinant non-LTR retroelement RT outperformed all commercially available kits and rivaled the low bias of technically demanding home-brew protocols. We applied OTTR to inventory RNAs purified from extracellular vesicles, identifying miRNAs as well as myriad other noncoding RNAs (ncRNAs) and ncRNA fragments. Our results establish the utility of OTTR for automation-friendly, low-bias, end-to-end RNA sequence inventories of complex ncRNA samples.

Quantitative mapping of the cellular small RNA landscape with AQRNA-seq

Current next-generation RNA-sequencing (RNA-seq) methods do not provide accurate quantification of small RNAs within a sample, due to sequence-dependent biases in capture, ligation and amplification during library preparation. We present a method, absolute quantification RNA-sequencing (AQRNA-seq), that minimizes biases and provides a direct, linear correlation between sequencing read count and copy number for all small RNAs in a sample. Library preparation and data processing were optimized and validated using a 963-member microRNA reference library, oligonucleotide standards of varying length, and RNA blots. Application of AQRNA-seq to a panel of human cancer cells revealed >800 detectable miRNAs that varied during cancer progression, while application to bacterial transfer RNA pools, with the challenges of secondary structure and abundant modifications, revealed 80-fold variation in tRNA isoacceptor levels, stress-induced site-specific tRNA fragmentation, quantitative modification maps, and evidence for stress-induced, tRNA-driven, codon-biased translation. AQRNA-seq thus provides a versatile means to quantitatively map the small RNA landscape in cells.

Structural basis for template switching by a group II intron-encoded non-LTR-retroelement reverse transcriptase

Reverse transcriptases (RTs) can switch template strands during complementary DNA synthesis, enabling them to join discontinuous nucleic acid sequences. Template switching (TS) plays crucial roles in retroviral replication and recombination, is used for adapter addition in RNA-Seq, and may contribute to retroelement fitness by increasing evolutionary diversity and enabling continuous complementary DNA synthesis on damaged templates. Here, we determined an X-ray crystal structure of a TS complex of a group II intron RT bound simultaneously to an acceptor RNA and donor RNA template–DNA primer heteroduplex with a 1-nt 3′-DNA overhang. The structure showed that the 3′ end of the acceptor RNA binds in a pocket formed by an N-terminal extension present in non–long terminal repeat–retroelement RTs and the RT fingertips loop, with the 3′ nucleotide of the acceptor base paired to the 1-nt 3′-DNA overhang and its penultimate nucleotide base paired to the incoming dNTP at the RT active site. Analysis of structure-guided mutations identified amino acids that contribute to acceptor RNA binding and a phenylalanine residue near the RT active site that mediates nontemplated nucleotide addition. Mutation of the latter residue decreased multiple sequential template switches in RNA-Seq. Our results provide new insights into the mechanisms of TS and nontemplated nucleotide addition by RTs, suggest how these reactions could be improved for RNA-Seq, and reveal common structural features for TS by non–long terminal repeat–retroelement RTs and viral RNA–dependent RNA polymerases.

SnoRNA copy regulation affects family size, genomic location and family abundance levels

Background

Small nucleolar RNAs (snoRNAs) are an abundant class of noncoding RNAs present in all eukaryotes and best known for their involvement in ribosome biogenesis. In mammalian genomes, many snoRNAs exist in multiple copies, resulting from recombination and retrotransposition from an ancestral snoRNA. To gain insight into snoRNA copy regulation, we used Rfam classification and normal human tissue expression datasets generated using low structure bias RNA-seq to characterize snoRNA families.

Results

We found that although box H/ACA families are on average larger than box C/D families, the number of expressed members is similar for both types. Family members can cover a wide range of average abundance values, but importantly, expression variability of individual members of a family is preferred over the total variability of the family, especially for box H/ACA snoRNAs, suggesting that while members are likely differentially regulated, mechanisms exist to ensure uniformity of the total family abundance across tissues. Box C/D snoRNA family members are mostly embedded in the same host gene while box H/ACA family members tend to be encoded in more than one different host, supporting a model in which box C/D snoRNA duplication occurred mostly by cis recombination while box H/ACA snoRNA families have gained copy members through retrotransposition. And unexpectedly, snoRNAs encoded in the same host gene can be regulated independently, as some snoRNAs within the same family vary in abundance in a divergent way between tissues.

Conclusions

SnoRNA copy regulation affects family sizes, genomic location of the members and controls simultaneously member and total family abundance to respond to the needs of individual tissues.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07757-1.

Changes of the tRNA Modification Pattern during the Development of Dictyostelium discoideum

Dictyostelium discoideum is a social amoeba, which on starvation develops from a single-cell state to a multicellular fruiting body. This developmental process is accompanied by massive changes in gene expression, which also affect non-coding RNAs. Here, we investigate how tRNAs as key regulators of the translation process are affected by this transition. To this end, we used LOTTE-seq to sequence the tRNA pool of D. discoideum at different developmental time points and analyzed both tRNA composition and tRNA modification patterns. We developed a workflow for the specific detection of modifications from reverse transcriptase signatures in chemically untreated RNA-seq data at single-nucleotide resolution. It avoids the comparison of treated and untreated RNA-seq data using reverse transcription arrest patterns at nucleotides in the neighborhood of a putative modification site as internal control. We find that nucleotide modification sites in D. discoideum tRNAs largely conform to the modification patterns observed throughout the eukaroytes. However, there are also previously undescribed modification sites. We observe substantial dynamic changes of both expression levels and modification patterns of certain tRNA types during fruiting body development. Beyond the specific application to D. discoideum our results demonstrate that the developmental variability of tRNA expression and modification can be traced efficiently with LOTTE-seq.

Tau aggregates are RNA-protein assemblies that mislocalize multiple nuclear speckle components

Tau aggregates contribute to neurodegenerative diseases, including frontotemporal dementia and Alzheimer's disease (AD). Although RNA promotes tau aggregation in vitro, whether tau aggregates in cells contain RNA is unknown. We demonstrate, in cell culture and mouse brains, that cytosolic and nuclear tau aggregates contain RNA with enrichment for small nuclear RNAs (snRNAs) and small nucleolar RNAs (snoRNAs). Nuclear tau aggregates colocalize with and alter the composition, dynamics, and organization of nuclear speckles, membraneless organelles involved in pre-mRNA splicing. Moreover, several nuclear speckle components, including SRRM2, mislocalize to cytosolic tau aggregates in cells, mouse brains, and brains of individuals with AD, frontotemporal dementia (FTD), and corticobasal degeneration (CBD). Consistent with these alterations, we observe that the presence of tau aggregates is sufficient to alter pre-mRNA splicing. This work identifies tau alteration of nuclear speckles as a feature of tau aggregation that may contribute to the pathology of tau aggregates.

High-resolution quantitative profiling of tRNA abundance and modification status in eukaryotes by mim-tRNAseq

Measurements of cellular tRNA abundance are hampered by pervasive blocks to cDNA synthesis at modified nucleosides and the extensive similarity among tRNA genes. We overcome these limitations with modification-induced misincorporation tRNA sequencing (mim-tRNAseq), which combines a workflow for full-length cDNA library construction from endogenously modified tRNA with a comprehensive and user-friendly computational analysis toolkit. Our method accurately captures tRNA abundance and modification status in yeast, fly, and human cells and is applicable to any organism with a known genome. We applied mim-tRNAseq to discover a dramatic heterogeneity of tRNA isodecoder pools among diverse human cell lines and a surprising interdependence of modifications at distinct sites within the same tRNA transcript.

•

mim-tRNAseq overcomes experimental and computational hurdles to tRNA quantitation

•

mim-tRNAseq includes a comprehensive computational toolkit for tRNA read analysis

•

tRNA abundance, aminoacylation, and modification status quantified in one reaction

•

mim-tRNAseq reveals an interdependence of modifications at distinct tRNA positions

Identification of protein-protected mRNA fragments and structured excised intron RNAs in human plasma by TGIRT-seq peak calling

Human plasma contains > 40,000 different coding and non-coding RNAs that are potential biomarkers for human diseases. Here, we used thermostable group II intron reverse transcriptase sequencing (TGIRT-seq) combined with peak calling to simultaneously profile all RNA biotypes in apheresis-prepared human plasma pooled from healthy individuals. Extending previous TGIRT-seq analysis, we found that human plasma contains largely fragmented mRNAs from > 19,000 protein-coding genes, abundant full-length, mature tRNAs and other structured small non-coding RNAs, and less abundant tRNA fragments and mature and pre-miRNAs. Many of the mRNA fragments identified by peak calling correspond to annotated protein-binding sites and/or have stable predicted secondary structures that could afford protection from plasma nucleases. Peak calling also identified novel repeat RNAs, miRNA-sized RNAs, and putatively structured intron RNAs of potential biological, evolutionary, and biomarker significance, including a family of full-length excised intron RNAs, subsets of which correspond to mirtron pre-miRNAs or agotrons.

Processing by RNase 1 forms tRNA halves and distinct Y RNA fragments in the extracellular environment

Extracellular RNAs participate in intercellular communication, and are being studied as promising minimally invasive diagnostic markers. Several studies in recent years showed that tRNA halves and distinct Y RNA fragments are abundant in the extracellular space, including in biofluids. While their regulatory and diagnostic potential has gained a substantial amount of attention, the biogenesis of these extracellular RNA fragments remains largely unexplored. Here, we demonstrate that these fragments are produced by RNase 1, a highly active secreted nuclease. We use RNA sequencing to investigate the effect of a null mutation of RNase 1 on the levels of tRNA halves and Y RNA fragments in the extracellular environment of cultured human cells. We complement and extend our RNA sequencing results with northern blots, showing that tRNAs and Y RNAs in the non-vesicular extracellular compartment are released from cells as full-length precursors and are subsequently cleaved to distinct fragments. In support of these results, formation of tRNA halves is recapitulated by recombinant human RNase 1 in our in vitro assay. These findings assign a novel function for RNase 1, and position it as a strong candidate for generation of tRNA halves and Y RNA fragments in biofluids.