Thermostable group II intron reverse transcriptase fusion proteins and their use in cDNA synthesis and next-generation RNA sequencing

Rapid folding of nascent RNA regulates eukaryotic RNA biogenesis

RNA's catalytic, regulatory, or coding potential depends on structure formation. Because base pairing occurs during transcription, early structural states can govern RNA processing events and dictate the formation of functional conformations. These co-transcriptional states remain mostly unknown. Here, we develop co-transcriptional structure tracking (CoSTseq), which detects nascent RNA base pairing within and upon exit from RNA polymerases (Pols) transcriptome wide in living yeast cells. Monitoring each nucleotide's base pairing activity during transcription, CoSTseq reveals predominantly rapid pairing-within 25 bp of transcription after addition to the nascent chain. Moreover, ∼23% of rRNA nucleotides attain their final base pairing state near Pol I, while most other nucleotides must undergo changes in pairing status during later steps of ribosome biogenesis. We show that helicases act immediately to remodel structures across the rDNA locus to facilitate ribosome biogenesis. By contrast, nascent pre-mRNAs attain local structures indistinguishable from mature mRNAs, suggesting that refolding behind elongating ribosomes resembles co-transcriptional folding behind Pol II.

Next-Generation Sequencing Methods to Determine the Accuracy of Retroviral Reverse Transcriptases: Advantages and Limitations

Retroviruses, like other RNA viruses, mutate at very high rates and exist as genetically heterogeneous populations. The error-prone activity of viral reverse transcriptase (RT) is largely responsible for the observed variability, most notably in HIV-1. In addition, RTs are widely used in biotechnology to detect RNAs and to clone expressed genes, among many other applications. The fidelity of retroviral RTs has been traditionally analyzed using enzymatic (gel-based) or reporter-based assays. However, these methods are laborious and have important limitations. The development of next-generation sequencing (NGS) technologies opened the possibility of obtaining reverse transcription error rates from a large number of sequences, although appropriate protocols had to be developed. In this review, we summarize the developments in this field that allowed the determination of RNA-dependent DNA synthesis error rates for different RTs (viral and non-viral), including methods such as PRIMER IDs, REP-SEQ, ARC-SEQ, CIR-SEQ, SMRT-SEQ and ROLL-SEQ. Their advantages and limitations are discussed. Complementary DNA (cDNA) synthesis error rates obtained in different studies, using RTs and RNAs of diverse origins, are presented and compared. Future improvements in methodological pipelines will be needed for the precise identification of mutations in the RNA template, including modified bases.

Genome-wide profiling of tRNA modifications by Induro-tRNAseq reveals coordinated changes

While all native tRNAs undergo extensive post-transcriptional modifications as a mechanism to regulate gene expression, mapping these modifications remains challenging. The critical barrier is the difficulty of readthrough of modifications by reverse transcriptases (RTs). Here we use Induro-a new group-II intron-encoded RT-to map and quantify genome-wide tRNA modifications in Induro-tRNAseq. We show that Induro progressively increases readthrough over time by selectively overcoming RT stops without altering the misincorporation frequency. In a parallel analysis of Induro vs. a related RT, we provide comparative datasets to facilitate the prediction of each modification. We assess tRNA modifications across five human cell lines and three mouse tissues and show that, while the landscape of modifications is highly variable throughout the tRNA sequence framework, it is stabilized for modifications that are required for reading of the genetic code. The coordinated changes have fundamental importance for development of tRNA modifications in protein homeostasis.

Structural basis for the evolution of a domesticated group II intron-like reverse transcriptase to function in host cell DNA repair

A previous study found that a domesticated bacterial group II intron-like reverse transcriptase (G2L4 RT) functions in double-strand break repair (DSBR) via microhomology-mediated end joining (MMEJ) and that a mobile group II intron-encoded RT has a basal DSBR activity that uses conserved structural features of non-LTR-retroelement RTs. Here, we determined G2L4 RT apoenzyme and snap-back DNA synthesis structures revealing novel structural adaptations that optimized its cellular function in DSBR. These included a unique RT3a structure that stabilizes the apoenzyme in an inactive conformation until encountering an appropriate substrate; a longer N-terminal extension/RT0-loop with conserved residues that together with a modified active site favors strand annealing; and a conserved dimer interface that localizes G2L4 RT homodimers to DSBR sites with both monomers positioned for MMEJ. Our findings reveal how a non-LTR-retroelement RT evolved a dedicated cellular function and suggest new ways of optimizing these RTs for genome engineering applications.

WITHDRAWN: Gammaherpesvirus infection alters transfer RNA splicing and triggers tRNA cleavage

The authors have withdrawn this manuscript due to a duplicate posting of manuscript number BIORXIV/2024/592122. Therefore, the authors do not wish this work to be cited as reference for the project. If you have any questions, please contact the corresponding author. The correct preprint can be found at doi: https://doi.org/10.1101/2024.05.01.592122 .

Small Nucleolar RNAs in Head and Neck Squamous Cell Carcinomas

Small nucleolar RNAs (snoRNAs), a distinct class of noncoding RNAs, encompass highly diverse structures and have a range of 60 to 300 nucleotides in length. About 90% of human snoRNAs are intronic and embedded within introns of their host gene transcripts. Most snoRNAs enriched in specific tissue correlate in abundance with their parental host genes. Advancements in high-throughput sequencing have facilitated the discovery of dysregulated snoRNA expression in numerous human malignancies including head and neck squamous cell carcinoma (HNSCC). Hundreds of differentially expressed snoRNAs have been identified in HNSCC tissues. Among 1,524 snoRNA genes in a 567 HNSCC cohort, 113 snoRNAs were found to be survival related. As for snoRNA's roles in HNSCC, based on the available evidence, dysregulated snoRNAs are closely associated with the carcinogenesis and development of HNSCC. Upregulated snoRNAs have been shown to augment the expression of other oncogenes or activate the Wnt/β-catenin signaling pathway, thereby promoting tumor cell viability, glycolysis, migration, and the epithelial-mesenchymal transition while inhibiting apoptosis in vitro. In vivo animal studies have further elucidated the functional roles of snoRNAs. Knockdown of host genes of these snoRNAs suppressed the Wnt/β-catenin signaling pathway and restrained tumor proliferation and aggressiveness in mice. The putative mechanisms underlying these observations are associated with the biological functions of snoRNAs, primarily involving microRNA-like functions through the generation of microRNA-like fragments and regulation of alternative splicing to yield diverse transcripts. While most of the snoRNAs are upregulated in HNSCC, 4 downregulated snoRNAs have been identified and annotated. SNORA36B (implicated in the regulation of DNA templates) and U3 (chr17, influencing cell proliferation) may serve as protective factors associated with prolonged overall survival. This review describes the viable structures of snoRNAs, endeavors to refine snoRNA sequencing technology, and summarizes snoRNAs' expression profile as well as their role in HNSCC progression for potential diagnostic and therapeutic strategies for HNSCC management.

Insights into the human cDNA: A descriptive study using library screening in yeast

The utilization of human cDNA libraries in yeast genetic screens is an approach that has been used to identify novel gene functions and/or genetic and physical interaction partners through forward genetics using yeast two-hybrid (Y2H) and classical cDNA library screens. Here, we summarize several challenges that have been observed during the implementation of human cDNA library screens in Saccharomyces cerevisiae (budding yeast). Upon the utilization of DNA repair deficient-yeast strains to identify novel genes that rescue the toxic effect of DNA-damage inducing drugs, we have observed a wide range of transcripts that could rescue the strains. However, after several rounds of screening, most of these hits turned out to be false positives, most likely due to spontaneous mutations in the yeast strains that arise as a rescue mechanism due to exposure to toxic DNA damage inducing-drugs. The observed transcripts included mitochondrial hits, non-coding RNAs, truncated cDNAs, and transcription products that resulted from the internal priming of genomic regions. We have also noticed that most cDNA transcripts are not fused with the GAL4 activation domain (GAL4AD), rendering them unsuitable for Y2H screening. Consequently, we utilized Sanger sequencing to screen 282 transcripts obtained from either four different yeast screens or through direct fishing from a human kidney cDNA library. The aim was to gain insights into the different transcription products and to highlight the challenges of cDNA screening approaches in the presence of a significant number of undesired transcription products. In summary, this study describes the challenges encountering human cDNA library screening in yeast as a valuable technique that led to the identification of important molecular mechanisms. The results open research venues to further optimize the process and increase its efficiency.

Rapid folding of nascent RNA regulates eukaryotic RNA biogenesis

An RNA's catalytic, regulatory, or coding potential depends on RNA structure formation. Because base pairing occurs during transcription, early structural states can govern RNA processing events and dictate the formation of functional conformations. These co-transcriptional states remain unknown. Here, we develop CoSTseq, which detects nascent RNA base pairing within and upon exit from RNA polymerases (Pols) transcriptome-wide in living yeast cells. By monitoring each nucleotide's base pairing activity during transcription, we identify distinct classes of behaviors. While 47% of rRNA nucleotides remain unpaired, rapid and delayed base pairing - with rates of 48.5 and 13.2 kb -1 of transcribed rDNA, respectively - typically completes when Pol I is only 25 bp downstream. We show that helicases act immediately to remodel structures across the rDNA locus and facilitate ribosome biogenesis. In contrast, nascent pre-mRNAs attain local structures indistinguishable from mature mRNAs, suggesting that refolding behind elongating ribosomes resembles co-transcriptional folding behind Pol II.

TGIRT-seq to profile tRNA-derived RNAs and associated RNA modifications

RNA modifications are key regulators for RNA processes. tRNA-derived RNAs are small RNAs with size between 15 and 50 bases long that are processed from mature or precursor tRNAs. Despite their more recent discovery, tRNA-derived RNAs have been found to play regulatory roles in many cellular processes including gene silencing, protein synthesis, stress response, and transgenerational inheritance. Furthermore, tRNA-derived RNAs are highly abundant in bodily fluids, posing as potential biomarkers. A unique feature of tRNA-derived RNAs is that they are rich in RNA modifications. Many of the RNA modifications on tRNA-derived RNAs disrupt Watson-Crick base pairing and will thus stall reverse transcriptase, such as N1-methyladenosine (m1A), N1-methylguanosine (m1G) and N2, N2-dimethylguanosine (m22G). These RNA modifications add another layer of regulation onto tRNA-derived RNAs' functions and are of interests for future research. However, these RNA modifications could also lead to lower detection of modification-containing RNAs in genome-wide small RNA sequencing analysis due to reverse transcriptase stall. To circumvent this bias, TGIRT (Thermostable Group II Intron Reverse Transcriptase) has been used to readthrough RNA modifications inserting mismatches. These mismatch signatures can then be used to precisely map the modification sites at base resolution. Here we describe the step-by-step experimental protocol to start with purified RNAs from cells or tissues and use TGIRT to make small RNA sequencing library for Illumina sequencing to profile the abundance of tRNA-derived RNAs and the associated RNA modifications.

Challenges in identifying mRNA transcript starts and ends from long-read sequencing data

Long-read sequencing (LRS) technologies have the potential to revolutionize scientific discoveries in RNA biology through the comprehensive identification and quantification of full-length mRNA isoforms. Despite great promise, challenges remain in the widespread implementation of LRS technologies for RNA-based applications, including concerns about low coverage, high sequencing error, and robust computational pipelines. Although much focus has been placed on defining mRNA exon composition and structure with LRS data, less careful characterization has been done of the ability to assess the terminal ends of isoforms, specifically, transcription start and end sites. Such characterization is crucial for completely delineating full mRNA molecules and regulatory consequences. However, there are substantial inconsistencies in both start and end coordinates of LRS reads spanning a gene, such that LRS reads often fail to accurately recapitulate annotated or empirically derived terminal ends of mRNA molecules. Here, we describe the specific challenges of identifying and quantifying mRNA terminal ends with LRS technologies and how these issues influence biological interpretations of LRS data. We then review recent experimental and computational advances designed to alleviate these problems, with ideal use cases for each approach. Finally, we outline anticipated developments and necessary improvements for the characterization of terminal ends from LRS data.

DMS-MapSeq Analysis of Antisense Oligonucleotide Binding to lncRNA PANDA

While various methods exist for examining and visualizing the structures of RNA molecules, dimethyl sulfate-mutational profiling and sequencing (DMS-MaPseq) stands out for its simplicity and versatility. This technique has been proven effective for studying RNA structures both in vitro and in complex biological settings. We present an updated protocol of DMS-MaPseq, as well as methodology that enables it to be used for detection of antisense oligonucleotides (ASOs) binding to RNA. By applying this protocol, we successfully characterized the structural ensemble of the HIV1 Rev Response Element (RRE), along with its two alternative structures. The findings align with previously published research validating the accuracy of the method. We also demonstrate the utility of the DMS-MaPseq protocol by resolving and confirming ASO binding at the complementary sites of the P21-associated noncoding RNA DNA damage-activated (PANDA) long non-coding RNA via decreased DMS reactivity. © 2024 Wiley Periodicals LLC. Basic Protocol 1: DMS-MaPseq on HIV1-RRE Basic Protocol 2: DMS-MaPseq on PANDA with ASO probing.

The dynamics and functional impact of tRNA repertoires during early embryogenesis in zebrafish

Embryogenesis entails dramatic shifts in mRNA translation and turnover that reprogram gene expression during cellular proliferation and differentiation. Codon identity modulates mRNA stability during early vertebrate embryogenesis, but how the composition of tRNA pools is matched to translational demand is unknown. By quantitative profiling of tRNA repertoires in zebrafish embryos during the maternal-to-zygotic transition, we show that zygotic tRNA repertoires are established after the onset of gastrulation, succeeding the major wave of zygotic mRNA transcription. Maternal and zygotic tRNA pools are distinct, but their reprogramming does not result in a better match to the codon content of the zygotic transcriptome. Instead, we find that an increase in global translation at gastrulation sensitizes decoding rates to tRNA supply, thus destabilizing maternal mRNAs enriched in slowly translated codons. Translational activation and zygotic tRNA expression temporally coincide with an increase of TORC1 activity at gastrulation, which phosphorylates and inactivates the RNA polymerase III repressor Maf1a/b. Our data indicate that a switch in global translation, rather than tRNA reprogramming, determines the onset of codon-dependent maternal mRNA decay during zebrafish embryogenesis.

TGIRT-seq of Inflammatory Breast Cancer Tumor and Blood Samples Reveals Widespread Enhanced Transcription Impacting RNA Splicing and Intronic RNAs in Plasma

Inflammatory breast cancer (IBC) is the most aggressive and lethal breast cancer subtype but lacks unequivocal genomic differences or robust biomarkers that differentiate it from non-IBC. Here, Thermostable Group II intron Reverse Transcriptase RNA-sequencing (TGIRT-seq) revealed myriad differences in tumor samples, Peripheral Blood Mononuclear Cells (PBMCs), and plasma that distinguished IBC from non-IBC patients and healthy donors across all tested receptor-based subtypes. These included numerous differentially expressed protein-coding gene and non-coding RNAs in all three sample types, a granulocytic immune response in IBC PBMCs, and over-expression of antisense RNAs, suggesting wide-spread enhanced transcription in both IBC tumors and PBMCs. By using TGIRT-seq to quantitate Intron-exon Depth Ratios (IDRs) and mapping reads to both genome and transcriptome reference sequences, we developed methods for parallel analysis of transcriptional and post-transcriptional gene regulation. This analysis identified numerous differentially and non-differentially expressed protein-coding genes in IBC tumors and PBMCs with high IDRs, the latter reflecting rate-limiting RNA splicing that negatively impacts mRNA production. Mirroring gene expression differences in tumors and PBMCs, over-represented protein-coding gene RNAs in IBC patient plasma were largely intronic RNAs, while those in non-IBC patients and healthy donor plasma were largely mRNA fragments. Potential IBC biomarkers in plasma included T-cell receptor pre-mRNAs and intronic, LINE-1, and antisense RNAs. Our findings provide new insights into IBC and set the stage for monitoring disease progression and response to treatment by liquid biopsy. The methods developed for parallel transcriptional and post-transcriptional gene regulation analysis have potentially broad RNA-seq and clinical applications.

Characterization and implementation of the MarathonRT template-switching reaction to expand the capabilities of RNA-seq

End-to-end RNA-sequencing methods that capture 5'-sequence content without cumbersome library manipulations are of great interest, particularly for analysis of long RNAs. While template-switching methods have been developed for RNA sequencing by distributive short-read RTs, such as the MMLV RTs used in SMART-Seq methods, they have not been adapted to leverage the power of ultraprocessive RTs, such as those derived from group II introns. To facilitate this transition, we dissected the individual processes that guide the enzymatic specificity and efficiency of the multistep template-switching reaction carried out by RTs, in this case, by MarathonRT. Remarkably, this is the first study of its kind, for any RT. First, we characterized the nucleotide specificity of nontemplated addition (NTA) reaction that occurs when the RT extends past the RNA 5'-terminus. We then evaluated the binding specificity of specialized template-switching oligonucleotides, optimizing their sequences and chemical properties to guide efficient template-switching reaction. Having dissected and optimized these individual steps, we then unified them into a procedure for performing RNA sequencing with MarathonRT enzymes, using a well-characterized RNA reference set. The resulting reads span a six-log range in transcript concentration and accurately represent the input RNA identities in both length and composition. We also performed RNA-seq from total human RNA and poly(A)-enriched RNA, with short- and long-read sequencing demonstrating that MarathonRT enhances the discovery of unseen RNA molecules by conventional RT. Altogether, we have generated a new pipeline for rapid, accurate sequencing of complex RNA libraries containing mixtures of long RNA transcripts.

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.

tRNA expression and modification landscapes, and their dynamics during zebrafish embryo development

tRNA genes exist in multiple copies in the genome of all organisms across the three domains of life. Besides the sequence differences across tRNA copies, extensive post-transcriptional modification adds a further layer to tRNA diversification. Whilst the crucial role of tRNAs as adapter molecules in protein translation is well established, whether all tRNAs are actually expressed, and whether the differences across isodecoders play any regulatory role is only recently being uncovered. Here we built upon recent developments in the use of NGS-based methods for RNA modification detection and developed tRAM-seq, an experimental protocol and in silico analysis pipeline to investigate tRNA expression and modification. Using tRAM-seq, we analysed the full ensemble of nucleo-cytoplasmic and mitochondrial tRNAs during embryonic development of the model vertebrate zebrafish. We show that the repertoire of tRNAs changes during development, with an apparent major switch in tRNA isodecoder expression and modification profile taking place around the start of gastrulation. Taken together, our findings suggest the existence of a general reprogramming of the expressed tRNA pool, possibly gearing the translational machinery for distinct stages of the delicate and crucial process of embryo development.

A robust method for measuring aminoacylation through tRNA-Seq

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

The modification landscape of Pseudomonas aeruginosa tRNAs

RNA modifications have a substantial impact on tRNA function, with modifications in the anticodon loop contributing to translational fidelity and modifications in the tRNA core impacting structural stability. In bacteria, tRNA modifications are crucial for responding to stress and regulating the expression of virulence factors. Although tRNA modifications are well-characterized in a few model organisms, our knowledge of tRNA modifications in human pathogens, such as Pseudomonas aeruginosa, remains limited. Here, we leveraged two orthogonal approaches to build a reference landscape of tRNA modifications in Escherichia coli, which enabled us to identify similar modifications in P. aeruginosa Our analysis supports a substantial degree of conservation between the two organisms, while also uncovering potential sites of tRNA modification in P. aeruginosa tRNAs that are not present in E. coli The mutational signature at one of these sites, position 46 of tRNAGln1(UUG) is dependent on the P. aeruginosa homolog of TapT, the enzyme responsible for the 3-(3-amino-3-carboxypropyl) uridine (acp3U) modification. Identifying which modifications are present on different tRNAs will uncover the pathways impacted by the different tRNA-modifying enzymes, some of which play roles in determining virulence and pathogenicity.

Perturbation of METTL1-mediated tRNA N7- methylguanosine modification induces senescence and aging

Cellular senescence is characterized by a decrease in protein synthesis, although the underlying processes are mostly unclear. Chemical modifications to transfer RNAs (tRNAs) frequently influence tRNA activity, which is crucial for translation. We describe how tRNA N7-methylguanosine (m7G46) methylation, catalyzed by METTL1-WDR4, regulates translation and influences senescence phenotypes. Mettl1/Wdr4 and m7G gradually diminish with senescence and aging. A decrease in METTL1 causes a reduction in tRNAs, especially those with the m7G modification, via the rapid tRNA degradation (RTD) pathway. The decreases cause ribosomes to stall at certain codons, impeding the translation of mRNA that is essential in pathways such as Wnt signaling and ribosome biogenesis. Furthermore, chronic ribosome stalling stimulates the ribotoxic and integrative stress responses, which induce senescence-associated secretory phenotype. Moreover, restoring eEF1A protein mitigates senescence phenotypes caused by METTL1 deficiency by reducing RTD. Our findings demonstrate that tRNA m7G modification is essential for preventing premature senescence and aging by enabling efficient mRNA translation.

Diversity of ribosomes at the level of rRNA variation associated with human health and disease

Ribosomal DNA and RNA (rDNA and rRNA) sequences are usually discarded from sequencing analyses. But with hundreds of copies of rDNA genes it is unknown whether they possess sequence variations that form different types of ribosomes that affect human physiology and disease. Here, we developed an algorithm for variant-calling between paralog genes (termed RGA) and compared rDNA variations found in short- and long-read sequencing data from the 1,000 Genomes Project (1KGP) and Genome In A Bottle (GIAB). We additionally developed a novel protocol for long-read sequencing full-length rRNA (RIBO-RT) from actively translating ribosomes. Our analyses identified hundreds of rDNA variants, most of which, surprisingly, are short insertion-deletions (indels) and dozens of highly abundant rRNA variants that are incorporated into translationally active ribosomes. To visualize variant ribosomes at the single cell level, we developed an in-situ rRNA sequencing method (SWITCH-seq) which revealed that variants are co-expressed within individual cells. Strikingly, by analyzing rDNA, we found that variants assemble into distinct ribosome subtypes. We discovered that these subtypes acquire different rRNA structures by successfully employing dimethyl sulfate (DMS) probing of full length rRNA. With this atlas we investigated rRNA variation changes across human tissues and cancer types. This revealed tissue-specific rRNA subtype expression in endoderm/ectoderm-derived tissues. In cancer, low abundant rRNA variants can become highly expressed, which suggests the presence of cancer-specific ribosomes. Together, this study identifies and comprehensively characterizes the diversity of ribosomes at the level of rRNA variants which is dominated by indel variants, their chromosomal location and unique structure as well as the association of ribosome variation with tissue-specific biology and cancer.

Isoform-specific RNA structure determination using Nano-DMS-MaP

RNA structure determination is essential to understand how RNA carries out its diverse biological functions. In cells, RNA isoforms are readily expressed with partial variations within their sequences due, for example, to alternative splicing, heterogeneity in the transcription start site, RNA processing or differential termination/polyadenylation. Nanopore dimethyl sulfate mutational profiling (Nano-DMS-MaP) is a method for in situ isoform-specific RNA structure determination. Unlike similar methods that rely on short sequencing reads, Nano-DMS-MaP employs nanopore sequencing to resolve the structures of long and highly similar RNA molecules to reveal their previously hidden structural differences. This Protocol describes the development and applications of Nano-DMS-MaP and outlines the main considerations for designing and implementing a successful experiment: from bench to data analysis. In cell probing experiments can be carried out by an experienced molecular biologist in 3-4 d. Data analysis requires good knowledge of command line tools and Python scripts and requires a further 3-5 d.

Gammaherpesvirus infection triggers the formation of tRNA fragments from premature tRNAs

Transfer RNAs (tRNAs) are fundamental for both cellular and viral gene expression during viral infection. In addition, mounting evidence supports biological function for tRNA cleavage products, including in the control of gene expression during conditions of stress and infection. We previously reported that infection with the model murine gammaherpesvirus, MHV68, leads to enhanced tRNA transcription. However, whether this has any influence on tRNA transcript processing, viral replication, or the host response is not known. Here, we combined two new approaches, sequencing library preparation by Ordered Two Template Relay (OTTR) and tRNA bioinformatic analysis by tRAX, to quantitatively profile full-length tRNAs and tRNA fragment (tRF) identities during MHV68 infection. We find that MHV68 infection triggers both pre-tRNA and mature tRNA cleavage, resulting in the accumulation of specific tRFs. OTTR-tRAX revealed not only host tRNAome changes, but also the expression patterns of virally-encoded tRNAs (virtRNAs) and virtRFs made from the MHV68 genome, including their base modification signatures. Because the transcript ends of several host tRFs matched tRNA splice junctions, we tested and confirmed the role of tRNA splicing factors TSEN2 and CLP1 in MHV68-induced tRF biogenesis. Further, we show that CLP1 kinase, and by extension tRNA splicing, is required for productive MHV68 infection. Our findings provide new insight into how gammaherpesvirus infection both impacts and relies on tRNA transcription and processing.

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.

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.

NAP-seq reveals multiple classes of structured noncoding RNAs with regulatory functions

Up to 80% of the human genome produces "dark matter" RNAs, most of which are noncapped RNAs (napRNAs) that frequently act as noncoding RNAs (ncRNAs) to modulate gene expression. Here, by developing a method, NAP-seq, to globally profile the full-length sequences of napRNAs with various terminal modifications at single-nucleotide resolution, we reveal diverse classes of structured ncRNAs. We discover stably expressed linear intron RNAs (sliRNAs), a class of snoRNA-intron RNAs (snotrons), a class of RNAs embedded in miRNA spacers (misRNAs) and thousands of previously uncharacterized structured napRNAs in humans and mice. These napRNAs undergo dynamic changes in response to various stimuli and differentiation stages. Importantly, we show that a structured napRNA regulates myoblast differentiation and a napRNA DINAP interacts with dyskerin pseudouridine synthase 1 (DKC1) to promote cell proliferation by maintaining DKC1 protein stability. Our approach establishes a paradigm for discovering various classes of ncRNAs with regulatory functions.

A robust method for measuring aminoacylation through tRNA-Seq

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

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.

The modification landscape of P. aeruginosa tRNAs

RNA modifications have a substantial impact on tRNA function, with modifications in the anticodon loop contributing to translational fidelity and modifications in the tRNA core impacting structural stability. In bacteria, tRNA modifications are crucial for responding to stress and regulating the expression of virulence factors. Although tRNA modifications are well-characterized in a few model organisms, our knowledge of tRNA modifications in human pathogens, such as Pseudomonas aeruginosa, remains limited. Here we leveraged two orthogonal approaches to build a reference landscape of tRNA modifications in E. coli, which enabled us to identify similar modifications in P. aeruginosa. Our analysis revealed a substantial degree of conservation between the two organisms, while also uncovering potential sites of tRNA modification in P. aeruginosa tRNAs that are not present in E. coli. The mutational signature at one of these sites, position 46 of tRNAGln1(UUG) is dependent on the P. aeruginosa homolog of TapT, the enzyme responsible for the 3-(3-amino-3-carboxypropyl) uridine (acp3U) modification. Identifying which modifications are present on different tRNAs will uncover the pathways impacted by the different tRNA modifying enzymes, some of which play roles in determining virulence and pathogenicity.

A Novel Thermostable and Processive Reverse Transcriptase from a Group II Intron of Anoxybacillus flavithermus

Reverse transcriptases (RTs) are a family of enzymes that synthesize DNA using an RNA template and are involved in retrovirus propagation and telomere lengthening. In vitro, RTs are widely applied in various methods, including RNA-seq, RT-PCR, and RT-LAMP. Thermostable RTs from bacterial group II introns are promising tools for biotechnology due to their higher thermostability, fidelity, and processivity compared to commonly used M-MuLV RT and its mutants. However, the diversity of group II intron-encoded RTs is still understudied. In this work, we biochemically characterized a novel RT from a thermophilic bacterium, Anoxybacillus flavithermus, which was isolated from a hot spring in New Zealand and has an optimal growth temperature of around 60 °C. The cloned RT, named Afl RT, retained approximately 40% of the specific activity after a 45 min incubation at 50 °C. The optimal pH was 8.5, the optimal temperature was between 45 and 50 °C, and Mn2+ ions were found to be an optimal cofactor. The processivity analysis with MS2 phage gRNA (3569 b) demonstrated that Afl RT elongated fully up to 36% of the template molecules. In reverse transcription and RT-qLAMP, the enzyme allowed up to 10 copies of MS2 phage genomic RNA to be detected per reaction. Thus, Afl RT holds great potential for a variety of practical applications that require the use of thermostable and processive RTs.

Genome-Wide Profiling of tRNA Using an Unexplored Reverse Transcriptase with High Processivity

Monitoring the dynamic changes of cellular tRNA pools is challenging, due to the extensive post-transcriptional modifications of individual species. The most critical component in tRNAseq is a processive reverse transcriptase (RT) that can read through each modification with high efficiency. Here we show that the recently developed group-II intron RT Induro has the processivity and efficiency necessary to profile tRNA dynamics. Using our Induro-tRNAseq, simpler and more comprehensive than the best methods to date, we show that Induro progressively increases readthrough of tRNA over time and that the mechanism of increase is selective removal of RT stops, without altering the misincorporation frequency. We provide a parallel dataset of the misincorporation profile of Induro relative to the related TGIRT RT to facilitate the prediction of non-annotated modifications. We report an unexpected modification profile among human proline isoacceptors, absent from mouse and lower eukaryotes, that indicates new biology of decoding proline codons.

Charge Engineering of the Nucleic Acid Binding Cleft of a Thermostable HIV-1 Reverse Transcriptase Reveals Key Interactions and a Novel Mechanism of RNase H Inactivation

Coupled with PCR, reverse transcriptases (RTs) have been widely used for RNA detection and gene expression analysis. Increased thermostability and nucleic acid binding affinity are desirable RT properties to improve yields and sensitivity of these applications. The effects of amino acid substitutions in the RT RNase H domain were tested in an engineered HIV-1 group O RT, containing mutations K358R/A359G/S360A and devoid of RNase H activity due to the presence of E478Q (O3MQ RT). Twenty mutant RTs with Lys or Arg at positions interacting with the template-primer (i.e., at positions 473-477, 499-502 and 505) were obtained and characterized. Most of them produced significant amounts of cDNA at 37, 50 and 65 °C, as determined in RT-PCR reactions. However, a big loss of activity was observed with mutants A477K/R, S499K/R, V502K/R and Y505K/R, particularly at 65 °C. Binding affinity experiments confirmed that residues 477, 502 and 505 were less tolerant to mutations. Amino acid substitutions Q500K and Q500R produced a slight increase of cDNA synthesis efficiency at 50 and 65 °C, without altering the KD for model DNA/DNA and RNA/DNA heteroduplexes. Interestingly, molecular dynamics simulations predicted that those mutations inactivate the RNase H activity by altering the geometry of the catalytic site. Proof of this unexpected effect was obtained after introducing Q500K or Q500R in the wild-type HIV-1BH10 RT and mutant K358R/A359G/S360A RT. Our results reveal a novel mechanism of RNase H inactivation that preserves RT DNA binding and polymerization efficiency without substituting RNase H active site residues.

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.

Reverse transcriptases prime DNA synthesis

The discovery of reverse transcriptases (RTs) challenged the central dogma by establishing that genetic information can also flow from RNA to DNA. Although they act as DNA polymerases, RTs are distantly related to replicases that also possess de novo primase activity. Here we identify that CRISPR associated RTs (CARTs) directly prime DNA synthesis on both RNA and DNA. We demonstrate that RT-dependent priming is utilized by some CRISPR-Cas complexes to synthesise new spacers and integrate these into CRISPR arrays. Expanding our analyses, we show that primer synthesis activity is conserved in representatives of other major RT classes, including group II intron RT, telomerase and retroviruses. Together, these findings establish a conserved innate ability of RTs to catalyse de novo DNA primer synthesis, independently of accessory domains or alternative priming mechanisms, which likely plays important roles in a wide variety of biological pathways.

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.

Nano-DMS-MaP allows isoform-specific RNA structure determination

Genome-wide measurements of RNA structure can be obtained using reagents that react with unpaired bases, leading to adducts that can be identified by mutational profiling on next-generation sequencing machines. One drawback of these experiments is that short sequencing reads can rarely be mapped to specific transcript isoforms. Consequently, information is acquired as a population average in regions that are shared between transcripts, thus blurring the underlying structural landscape. Here, we present nanopore dimethylsulfate mutational profiling (Nano-DMS-MaP)-a method that exploits long-read sequencing to provide isoform-resolved structural information of highly similar RNA molecules. We demonstrate the value of Nano-DMS-MaP by resolving the complex structural landscape of human immunodeficiency virus-1 transcripts in infected cells. We show that unspliced and spliced transcripts have distinct structures at the packaging site within the common 5' untranslated region, likely explaining why spliced viral RNAs are excluded from viral particles. Thus, Nano-DMS-MaP is a straightforward method to resolve biologically important transcript-specific RNA structures that were previously hidden in short-read ensemble analyses.

Functional analysis of tRNA modification enzymes using mutational profiling

Transfer RNA (tRNA) delivers amino acids to the ribosome and functions as an essential adapter molecule for decoding codons on the messenger RNA (mRNA) during protein synthesis. Before attaining their proper activity, tRNAs undergo multiple post-transcriptional modifications with highly diversified roles such as stabilization of the tRNA structure, recognition of aminoacyl tRNA synthetases, precise codon-anticodon recognition, support of viral replication and onset of immune responses. The synthesis of the majority of modified nucleosides is catalyzed by a site-specific tRNA modification enzyme. This chapter provides a detailed protocol for using mutational profiling to analyze the enzymatic function of a tRNA methyltransferase in a high-throughput manner. In a previous study, we took tRNA m1A22 methyltransferase TrmK from Geobacillus stearothermophilus as a model tRNA methyltransferase and applied this protocol to gain mechanistic insights into how TrmK recognizes the substrate tRNAs. In theory, this protocol can be used unaltered for studying enzymes that catalyze modifications at the Watson-Crick face such as 1-methyladenosine (m1A), 3-methylcytosine (m3C), 3-methyluridine (m3U), 1-methylguanosine (m1G), and N2,N2-dimethylguanosine (m22G).

Reverse Transcriptases: From Discovery and Applications to Xenobiology

Reverse transcriptases are DNA polymerases that can use RNA as a template for DNA synthesis. They thus catalyze the reverse of transcription. Although discovered in 1970, reverse transcriptases are still of great interest and are constantly being further developed for numerous modern research approaches. They are frequently used in biotechnological and molecular diagnostic applications. In this review, we describe the discovery of these fascinating enzymes and summarize research results and applications ranging from molecular cloning, direct virus detection, and modern sequencing methods to xenobiology.

Engineered CRISPR prime editors with compact, untethered reverse transcriptases

The CRISPR prime editor PE2 consists of a Streptococcus pyogenes Cas9 nickase (nSpCas9) fused at its C-terminus to a Moloney murine leukemia virus reverse transcriptase (MMLV-RT). Here we show that separated nSpCas9 and MMLV-RT proteins function as efficiently as intact PE2 in human cells. We use this Split-PE system to rapidly identify and engineer more compact prime editor architectures that also broaden the types of RTs used for prime editing.

ALL-tRNAseq enables robust tRNA profiling in tissue samples

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

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.

Comparison of Reverse Transcriptase (RT) Activities of Various M-MuLV RTs for RT-LAMP Assays

Reverse transcriptases (RTs) are a family of enzymes synthesizing DNA using RNA as a template and serving as indispensable tools in studies related to RNA. M-MuLV RT and its analogs are the most commonly used RTs. RTs are widely applied in various diagnostics methods, including reverse-transcription loop-mediated isothermal amplification (RT-LAMP). However, the performance of different RTs in LAMP remains relatively unknown. Here, we report on the first direct comparison of various M-MuLV RTs in RT-LAMP, including enzymes with a different number of mutations and fusions with Sto7d. Several parameters were assessed, namely: optimal reaction temperature, enzyme concentration, reverse transcription time, a minimal amount of RNA template, and tolerance to inhibitors. Mutations increased the optimal reaction temperature from 55 °C to 60-65 °C. All of the RTs were suitable for RT-LAMP with RNA templates in the range of 10¹-10⁶ copies per reaction. Highly mutated enzymes were 1.5-3-fold more tolerant to whole blood, blood plasma, and guanidinium, but they were two-fold more sensitive to high concentrations of NaCl. The comparison of different RTs presented here could be helpful for selecting the optimal enzyme when developing novel LAMP-based diagnostic tests.

Two isoleucyl tRNAs that decode synonymous codons divergently regulate breast cancer metastatic growth by controlling translation of proliferation-regulating genes

The human genome contains 61 codons encoding 20 amino acids. Synonymous codons representing a given amino acid are decoded by a set of transfer RNAs (tRNAs) called isoacceptors. We report the surprising observation that two isoacceptor tRNAs that decode synonymous codons become modulated in opposing directions during breast cancer progression. Specifically, tRNAIleUAU became upregulated, whereas tRNAIleGAU became repressed as breast cancer cells attained enhanced metastatic capacity. Functionally, tRNAIleUAU promoted and tRNAIleGAU suppressed metastatic colonization in mouse xenograft models. These tRNAs mediated opposing effects on codon-dependent translation of growth-promoting genes, consistent with genomic enrichment or depletion of their cognate codons in mitotic genes. Our findings uncover a specific isoacceptor tRNA pair that act in opposition, divergently impacting growth-regulating genes and a disease phenotype. Degeneracy of the genetic code can thus be biologically exploited by human cancer cells via tRNA isoacceptor shifts that causally facilitate the transition toward a growth-promoting state.

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.

Experimental and computational workflow for the analysis of tRNA pools from eukaryotic cells by mim-tRNAseq

Quantifying tRNAs is crucial for understanding how they regulate mRNA translation but is hampered by their extensive sequence similarity and premature termination of reverse transcription at multiple modified nucleotides. Here, we describe the use of modification-induced misincorporation tRNA sequencing (mim-tRNAseq), which overcomes these limitations with optimized library construction and a comprehensive toolkit for data analysis and visualization. We outline algorithm improvements that enhance the efficiency and accuracy of read alignment and provide details on data analysis outputs using example datasets.

For complete details on the use and execution of this protocol, please refer to Behrens et al. (2021).

•

mim-tRNAseq enables efficient quantitation and analysis of eukaryotic tRNA pools

•

Combined library generation and computational package ensures user-friendly workflow

•

Improved multiplexing enables more cost-effective sample analysis

•

Novel deconvolution algorithm enhances read alignment accuracy and sensitivity

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

A split prime editor with untethered reverse transcriptase and circular RNA template

Delivery and optimization of prime editors (PEs) have been hampered by their large size and complexity. Although split versions of genome-editing tools can reduce construct size, they require special engineering to tether the binding and catalytic domains. Here we report a split PE (sPE) in which the Cas9 nickase (nCas9) remains untethered from the reverse transcriptase (RT). The sPE showed similar efficiencies in installing precise edits as the parental unsplit PE3 and no increase in insertion-deletion (indel) byproducts. Delivery of sPE to the mouse liver with hydrodynamic injection to modify β-catenin drove tumor formation with similar efficiency as PE3. Delivery with two adeno-associated virus (AAV) vectors corrected the disease-causing mutation in a mouse model of type I tyrosinemia. Similarly, prime editing guide RNAs (pegRNAs) can be split into a single guide RNA (sgRNA) and a circular RNA RT template to increase flexibility and stability. Compared to previous sPEs, ours lacks inteins, protein-protein affinity modules and nuclease-sensitive pegRNA extensions, which increase construct complexity and might reduce efficiency. Our modular system will facilitate the delivery and optimization of PEs.

Human dyskerin binds to cytoplasmic H/ACA-box-containing transcripts affecting nuclear hormone receptor dependence

Background

Dyskerin is a nuclear protein involved in H/ACA box snoRNA-guided uridine modification of RNA. In humans, its defective function is associated with cancer development and induces specific post-transcriptional alterations of gene expression. In this study, we seek to unbiasedly identify mRNAs regulated by dyskerin in human breast cancer-derived cells.

Results

We find that dyskerin depletion affects the expression and the association with polysomes of selected mRNA isoforms characterized by the retention of H/ACA box snoRNA-containing introns. These snoRNA retaining transcripts (snoRTs) are bound by dyskerin in the cytoplasm in the form of shorter 3′ snoRT fragments. We then characterize the whole cytoplasmic dyskerin RNA interactome and find both H/ACA box snoRTs and protein-coding transcripts which may be targeted by the snoRTs’ guide properties. Since a fraction of these protein-coding transcripts is involved in the nuclear hormone receptor binding, we test to see if this specific activity is affected by dyskerin. Obtained results indicate that dyskerin dysregulation may alter the dependence on nuclear hormone receptor ligands in breast cancer cells. These results are paralleled by consistent observations on the outcome of primary breast cancer patients stratified according to their tumor hormonal status. Accordingly, experiments in nude mice show that the reduction of dyskerin levels in estrogen-dependent cells favors xenograft development in the absence of estrogen supplementation.

Conclusions

Our work suggests a cytoplasmic function for dyskerin which could affect mRNA post-transcriptional networks relevant for nuclear hormone receptor functions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02746-3.

Direct tracking of reverse-transcriptase speed and template sensitivity: implications for sequencing and analysis of long RNA molecules

Although reverse-transcriptase (RT) enzymes are critical reagents for research and biotechnology, their mechanical properties are not well understood. In particular, we know little about their relative speed and response to structural obstacles in the template. Commercial retroviral RTs stop at many positions along mixed sequence templates, resulting in truncated cDNA products that complicate downstream analysis. By contrast, group II intron-encoded RTs appear to copy long RNAs with high processivity and minimal stops. However, their speed, consistency and pausing behavior have not been explored. Here, we analyze RT velocity as the enzyme moves through heterogeneous sequences and structures that are embedded within a long noncoding RNA transcript. We observe that heterogeneities in the template are highly disruptive to primer extension by retroviral RTs. However, sequence composition and template structure have negligible effects on behavior of group II intron RTs, such as MarathonRT (MRT). Indeed, MRT copies long RNAs in a single pass, and displays synchronized primer extension at a constant speed of 25 nt/sec. In addition, it passes through stable RNA structural motifs without perturbation of velocity. Taken together, the results demonstrate that consistent, robust translocative behavior is a hallmark of group II intron-encoded RTs, some of which operate at high velocity.

Integrated Transcriptome and Metabolomic Analysis Reveal Anti-Angiogenic Properties of Disarib, a Novel Bcl2-Specific Inhibitor

Transcriptomic profiling of several drugs in cancer cell lines has been utilised to obtain drug-specific signatures and guided combination therapy to combat drug resistance and toxicity. Global metabolomics reflects changes due to altered activity of enzymes, environmental factors, etc. Integrating transcriptomics and metabolomics can provide genotype-phenotype correlation, providing meaningful insights into alterations in gene expression and its outcome to understand differential metabolism and guide therapy. This study uses a multi-omics approach to understand the global gene expression and metabolite changes induced by Disarib, a novel Bcl2-specific inhibitor in the Ehrlich adenocarcinoma (EAC) breast cancer mouse model. RNAseq analysis was performed on EAC mouse tumours treated with Disarib and compared to the controls. The expression of 6 oncogenes and 101 tumour suppressor genes interacting with Bcl2 and Bak were modulated upon Disarib treatment. Cancer hallmark pathways like DNA repair, Cell cycle, angiogenesis, and mitochondrial metabolism were downregulated, and programmed cell death platelet-related pathways were upregulated. Global metabolomic profiling using LC-MS revealed that Oncometabolites like carnitine, oleic acid, glycine, and arginine were elevated in tumour mice compared to normal and were downregulated upon Disarib treatment. Integrated transcriptomic and metabolomic profiles identified arginine metabolism, histidine, and purine metabolism to be altered upon Disarib treatment. Pro-angiogenic metabolites, arginine, palmitic acid, oleic acid, and myristoleic acid were downregulated in Disarib-treated mice. We further validated the effect of Disarib on angiogenesis by qRT-PCR analysis of genes in the VEGF pathway. Disarib treatment led to the downregulation of pro-angiogenic markers. Furthermore, the chorioallantoic membrane assay displayed a reduction in the formation of the number of secondary blood vessels upon Disarib treatment. Disarib reduces tumours by reducing oncometabolite and activating apoptosis and downregulating angiogenesis.

Meiotic cDNA libraries reveal gene truncations and mitochondrial proteins important for competitive fitness in Saccharomyces cerevisiae

Gametogenesis is an evolutionarily conserved developmental program whereby a diploid progenitor cell undergoes meiosis and cellular remodeling to differentiate into haploid gametes, the precursors for sexual reproduction. Even in the simple eukaryotic organism Saccharomyces cerevisiae, the meiotic transcriptome is very rich and complex, thereby necessitating new tools for functional studies. Here, we report the construction of 5 stage-specific, inducible complementary DNA libraries from meiotic cells that represent over 84% of the genes found in the budding yeast genome. We employed computational strategies to detect endogenous meiotic transcript isoforms as well as library-specific gene truncations. Furthermore, we developed a robust screening pipeline to test the effect of each complementary DNA on competitive fitness. Our multiday proof-of-principle time course revealed 877 complementary DNAs that were detrimental for competitive fitness when overexpressed. The list included mitochondrial proteins that cause dose-dependent disruption of cellular respiration as well as library-specific gene truncations that expose a dominant negative effect on competitive growth. Together, these high-quality complementary DNA libraries provide an important tool for systematically identifying meiotic genes, transcript isoforms, and protein domains that are important for a specific biological function.

Development of Mild Chemical Catalysis Conditions for m1A-to-m6A Rearrangement on RNA

The conversion of N¹-methyladenosine (m¹A) to N⁶-methyladenosine (m⁶A) on RNA is an important step for both allowing efficient reverse transcription read-though for sequencing analysis and mapping modifications in the transcriptome. Enzymatic transformation is often used, but the efficiency of the removal can depend on local sequence context. Chemical conversion through the application of the Dimroth rearrangement, in which m¹A rearranges into m⁶A under heat and alkaline conditions, is an alternative, but the required alkaline conditions result in significant RNA degradation by hydrolysis of the phosphodiester backbone. Here, we report novel, mild pH conditions that catalyze m¹A-to-m⁶A arrangement using 4-nitrothiophenol as a catalyst. We demonstrate the efficient rearrangement in mononucleosides, synthetic RNA oligonucleotides, and RNAs isolated from human cell lines, thereby validating a new approach for converting m¹A-to-m⁶A in RNA samples for sequencing analyses.