An optimized kit-free method for making strand-specific deep sequencing libraries from RNA fragments

Non-AUG HIV-1 uORF translation elicits specific T cell immune response and regulates viral transcript expression

Human immunodeficiency virus type-1 (HIV-1) is a complex retrovirus that relies on alternative splicing, translational, and post-translational mechanisms to produce over 15 functional proteins from its single ~10 kb transcriptional unit. Using ribosome profiling, nascent protein labeling, RNA sequencing, and whole-proteomics of infected CD4 + T lymphocytes, we characterized the transcriptional, translational, and post-translational landscape during infection. While viral infection exerts a significant impact on host transcript abundance, global translation rates are only modestly affected. Proteomics data reveal extensive transcriptional and post-translational regulation, with many genes showing opposing trends between transcript/ribosome profiling and protein abundance. These findings highlight a complex regulatory network orchestrating gene expression at multiple levels. Viral ribosome profiling further uncovered extensive non-AUG translation of small peptides from upstream open reading frames (uORFs) within the 5' long terminal repeat, which elicit specific T cell responses in people living with HIV. Conservation of uORF translation among retroviruses, along with TAR sequences, shapes DDX3 dependency for efficient translation of the main viral open reading frames.

Temperature-Dependent tRNA Modifications in Bacillales

Transfer RNA (tRNA) modifications are essential for the temperature adaptation of thermophilic and psychrophilic organisms as they control the rigidity and flexibility of transcripts. To further understand how specific tRNA modifications are adjusted to maintain functionality in response to temperature fluctuations, we investigated whether tRNA modifications represent an adaptation of bacteria to different growth temperatures (minimal, optimal, and maximal), focusing on closely related psychrophilic (P. halocryophilus and E. sibiricum), mesophilic (B. subtilis), and thermophilic (G. stearothermophilus) Bacillales. Utilizing an RNA sequencing approach combined with chemical pre-treatment of tRNA samples, we systematically profiled dihydrouridine (D), 4-thiouridine (s4U), 7-methyl-guanosine (m7G), and pseudouridine (Ψ) modifications at single-nucleotide resolution. Despite their close relationship, each bacterium exhibited a unique tRNA modification profile. Our findings revealed increased tRNA modifications in the thermophilic bacterium at its optimal growth temperature, particularly showing elevated levels of s4U8 and Ψ55 modifications compared to non-thermophilic bacteria, indicating a temperature-dependent regulation that may contribute to thermotolerance. Furthermore, we observed higher levels of D modifications in psychrophilic and mesophilic bacteria, indicating an adaptive strategy for cold environments by enhancing local flexibility in tRNAs. Our method demonstrated high effectiveness in identifying tRNA modifications compared to an established tool, highlighting its potential for precise tRNA profiling studies.

Metabolism-dependent secondary effect of anti-MAPK cancer therapy on DNA repair

Amino acid bioavailability impacts mRNA translation in a codon-dependent manner. Here, we report that the anti-cancer MAPK inhibitors (MAPKi) decrease the intracellular concentration of aspartate and glutamate in melanoma cells. This coincides with the accumulation of ribosomes on codons corresponding to these amino acids and triggers the translation-dependent degradation of mRNAs encoding aspartate- and glutamate-rich proteins, involved in DNA metabolism such as DNA replication and repair. Consequently, cells that survive MAPKi degrade aspartate and glutamate likely to generate energy, which simultaneously decreases their requirement for amino acids due to the downregulation of aspartate- and glutamate-rich proteins involved in cell proliferation. Concomitantly, the downregulation of aspartate- and glutamate-rich proteins involved in DNA repair increases DNA damage loads. Thus, DNA repair defects, and therefore mutations, are at least in part a secondary effect of the metabolic adaptation of cells exposed to MAPKi.

Stable structures or PABP1 loading protects cellular and viral RNAs against ISG20-mediated decay

ISG20 is an IFN-induced 3'-5' RNA exonuclease that acts as a broad antiviral factor. At present, the features that expose RNA to ISG20 remain unclear, although recent studies have pointed to the modulatory role of epitranscriptomic modifications in the susceptibility of target RNAs to ISG20. These findings raise the question as to how cellular RNAs, on which these modifications are abundant, cope with ISG20. To obtain an unbiased perspective on this topic, we used RNA-seq and biochemical assays to identify elements that regulate the behavior of RNAs against ISG20. RNA-seq analyses not only indicate a general preservation of the cell transcriptome, but they also highlight a small, but detectable, decrease in the levels of histone mRNAs. Contrarily to all other cellular ones, histone mRNAs are non-polyadenylated and possess a short stem-loop at their 3' end, prompting us to examine the relationship between these features and ISG20 degradation. The results we have obtained indicate that poly(A)-binding protein loading on the RNA 3' tail provides a primal protection against ISG20, easily explaining the overall protection of cellular mRNAs observed by RNA-seq. Terminal stem-loop RNA structures have been associated with ISG20 protection before. Here, we re-examined this question and found that the balance between resistance and susceptibility to ISG20 depends on their thermodynamic stability. These results shed new light on the complex interplay that regulates the susceptibility of different classes of viruses against ISG20.

Translation-dependent and -independent mRNA decay occur through mutually exclusive pathways defined by ribosome density during T cell activation

mRNA translation and decay are tightly interconnected processes both in the context of mRNA quality-control pathways and for the degradation of functional mRNAs. Cotranslational mRNA degradation through codon usage, ribosome collisions, and the recruitment of specific proteins to ribosomes is an important determinant of mRNA turnover. However, the extent to which translation-dependent mRNA decay (TDD) and translation-independent mRNA decay (TID) pathways participate in the degradation of mRNAs has not been studied yet. Here we describe a comprehensive analysis of basal and signal-induced TDD and TID in mouse primary CD4+ T cells. Our results indicate that most cellular transcripts are decayed to some extent in a translation-dependent manner. Our analysis further identifies the length of untranslated regions, the density of ribosomes, and GC3 content as important determinants of TDD magnitude. Consistently, all transcripts that undergo changes in ribosome density within their coding sequence upon T cell activation display a corresponding change in their TDD level. Moreover, we reveal a dynamic modulation in the relationship between GC3 content and TDD upon T cell activation, with a reversal in the impact of GC3- and AU3-rich codons. Altogether, our data show a strong and dynamic interconnection between mRNA translation and decay in mammalian primary cells.

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.

Advances in methods for tRNA sequencing and quantification

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

Host-Directed Targeting of LincRNA-MIR99AHG Suppresses Intracellular Growth of Mycobacterium tuberculosis

Tuberculosis (TB) caused by Mycobacterium tuberculosis (Mtb) kills 1.6 million people worldwide every year, and there is an urgent need for targeting host-pathogen interactions as a strategy to reduce mycobacterial resistance to current antimicrobials. Noncoding RNAs are emerging as important regulators of numerous biological processes and avenues for exploitation in host-directed therapeutics. Although long noncoding RNAs (lncRNAs) are abundantly expressed in immune cells, their functional role in gene regulation and bacterial infections remains understudied. In this study, we identify an immunoregulatory long intergenic noncoding RNA, lincRNA-MIR99AHG, which is upregulated in mouse and human macrophages upon IL-4/IL-13 stimulation and downregulated after clinical Mtb HN878 strain infection and in peripheral blood mononuclear cells from active TB patients. To evaluate the functional role of lincRNA-MIR99AHG, we used antisense locked nucleic acid (LNA) GapmeR-mediated antisense oligonucleotide (ASO) lncRNA knockdown experiments. Knockdown of lincRNA-MIR99AHG with ASOs significantly reduced intracellular Mtb growth in mouse and human macrophages and reduced pro-inflammatory cytokine production. In addition, in vivo treatment of mice with MIR99AHG ASOs reduced the mycobacterial burden in the lung and spleen. Furthermore, in macrophages, lincRNA-MIR99AHG is translocated to the nucleus and interacts with high affinity to hnRNPA2/B1 following IL-4/IL-13 stimulation and Mtb HN878 infection. Together, these findings identify lincRNA-MIR99AHG as a positive regulator of inflammation and macrophage polarization to promote Mtb growth and a possible target for adjunctive host-directed therapy against TB.

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.

Ribosome profiling reveals novel regulation of C9ORF72 GGGGCC repeat-containing RNA translation

GGGGCC (G₄C₂) repeat expansion in the first intron of C9ORF72 causes amyotrophic lateral sclerosis and frontotemporal dementia. Repeat-containing RNA is translated into dipeptide repeat (DPR) proteins, some of which are neurotoxic. Using dynamic ribosome profiling, we identified three translation initiation sites in the intron upstream of (G₄C₂) repeats; these sites are detected irrespective of the presence or absence of the repeats. During translocation, ribosomes appear to be stalled on the repeats. An AUG in the preceding C9ORF72 exon initiates a uORF that inhibits downstream translation. Polysome isolation indicates that unspliced (G₄C₂) repeat-containing RNA is a substrate for DPR protein synthesis. (G₄C₂) repeat-containing RNA translation is 5' cap-independent but inhibited by the initiation factor DAP5, suggesting an interplay with uORF function. These results define novel translational mechanisms of expanded (G₄C₂) repeat-containing RNA in disease.

RIPiT-Seq: A tandem immunoprecipitation approach to reveal global binding landscape of multisubunit ribonucleoproteins

RNA-binding proteins (RBPs) regulate all aspects of RNA metabolism. The ability to identify RNA targets bound by RBPs is critical for understanding RBP function. While powerful techniques are available to identify binding sites of individual RBPs at high resolution, it remains challenging to unravel binding sites of multicomponent ribonucleoproteins (RNPs) where multiple RBPs or proteins function cooperatively to bind to target RNAs. To fill this gap, we have previously developed RNA Immunoprecipitation in Tandem followed by high-throughput sequencing (RIPiT-seq) to characterize RNA targets of compositionally distinct RNP complexes by sequentially immunoprecipitating two proteins from the same RNP and sequencing the co-purifying RNA footprints. Here, we provide an updated and improved protocol for RIPiT-seq. In this protocol, we have used CRISPR-Cas9 to introduce affinity tag to endogenous protein of interest to capture a more representative state of an RNP complex. We present a modified protocol for library preparation for high-throughput sequencing so that it exclusively uses equipment and reagents available in a standard molecular biology lab. This updated custom library preparation protocol is compatible with commercial PCR multiplexing systems for Illumina sequencing platform for simultaneous and cost-effective analysis of large number of samples.

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

FMRP links optimal codons to mRNA stability in neurons

Fragile X syndrome (FXS) is caused by inactivation of the FMR1 gene and loss of encoded FMRP, an RNA binding protein that represses translation of some of its target transcripts. Here we use ribosome profiling and RNA sequencing to investigate the dysregulation of translation in the mouse brain cortex. We find that most changes in ribosome occupancy on hundreds of mRNAs are largely driven by dysregulation in transcript abundance. Many down-regulated mRNAs, which are mostly responsible for neuronal and synaptic functions, are highly enriched for FMRP binding targets. RNA metabolic labeling demonstrates that, in FMRP-deficient cortical neurons, mRNA down-regulation is caused by elevated degradation and is correlated with codon optimality. Moreover, FMRP preferentially binds mRNAs with optimal codons, suggesting that it stabilizes such transcripts through direct interactions via the translational machinery. Finally, we show that the paradigm of genetic rescue of FXS-like phenotypes in FMRP-deficient mice by deletion of the Cpeb1 gene is mediated by restoration of steady-state RNA levels and consequent rebalancing of translational homeostasis. Our data establish an essential role of FMRP in codon optimality-dependent mRNA stability as an important factor in FXS.

Retrograde signaling by a mtDNA-encoded non-coding RNA preserves mitochondrial bioenergetics

Alveolar epithelial type II (AETII) cells are important for lung epithelium maintenance and function. We demonstrate that AETII cells from mouse lungs exposed to cigarette smoke (CS) increase the levels of the mitochondria-encoded non-coding RNA, mito-RNA-805, generated by the control region of the mitochondrial genome. The protective effects of mito-ncR-805 are associated with positive regulation of mitochondrial energy metabolism, and respiration. Levels of mito-ncR-805 do not relate to steady-state transcription or replication of the mitochondrial genome. Instead, CS-exposure causes the redistribution of mito-ncR-805 from mitochondria to the nucleus, which correlated with the increased expression of nuclear-encoded genes involved in mitochondrial function. These studies reveal an unrecognized mitochondria stress associated retrograde signaling, and put forward the idea that mito-ncRNA-805 represents a subtype of small non coding RNAs that are regulated in a tissue- or cell-type specific manner to protect cells under physiological stress.

Ribosome profiling in mouse hippocampus: plasticity-induced regulation and bidirectional control by TSC2 and FMRP

BACKGROUND

Mutations in TSC2 are the most common cause of tuberous sclerosis (TSC), a disorder with a high incidence of autism and intellectual disability. TSC2 regulates mRNA translation required for group 1 metabotropic glutamate receptor-dependent synaptic long-term depression (mGluR-LTD) and behavior, but the identity of mRNAs responsive to mGluR-LTD signaling is largely unknown.

METHODS

We utilized Tsc2+/- mice as a mouse model of TSC and prepared hippocampal slices from these animals. We induced mGluR-LTD synaptic plasticity in slices and processed the samples for RNA-seq and ribosome profiling to identify differentially expressed genes in Tsc2+/- and following mGluR-LTD synaptic plasticity.

RESULTS

Ribosome profiling reveals that in Tsc2+/- mouse hippocampal slices, the expression of several mRNAs was dysregulated: terminal oligopyrimidine (TOP)-containing mRNAs decreased, while FMRP-binding targets increased. Remarkably, we observed the opposite changes of FMRP binding targets in Fmr1-/y hippocampi. In wild-type hippocampus, induction of mGluR-LTD caused rapid changes in the steady-state levels of hundreds of mRNAs, many of which are FMRP targets. Moreover, mGluR-LTD failed to promote phosphorylation of eukaryotic elongation factor 2 (eEF2) in TSC mice, and chemically mimicking phospho-eEF2 with low cycloheximide enhances mGluR-LTD in TSC mice.

CONCLUSION

These results suggest a molecular basis for bidirectional regulation of synaptic plasticity and behavior by TSC2 and FMRP. Our study also suggests that altered mGluR-regulated translation elongation contributes to impaired synaptic plasticity in Tsc2+/- mice.

Lnc-FAM84B-4 acts as an oncogenic lncRNA by interacting with protein hnRNPK to restrain MAPK phosphatases-DUSP1 expression

The mitogen activated protein kinase (MAPK) pathway has been reported to be involved in many cancer developments. Normally, MAPK activity is self-limited between rapid phosphorylation and dephosphorylation. In abnormal conditions, however, this dynamic equilibrium is broken, trigging tumor-suppressing or -promoting roles. While dual-specificity MAPK phosphatases (MKP/DUSPs) are important for cascade control in MAPK pathway, their role in colorectal cancer (CRC) remains largely unknown. Here, we investigated lnc-FAM84B-4 and DUSP1 to systematically elucidate their underlying roles in MAPK singling pathway and functions in CRC. Upregulated lnc-FAM84B-4 was identified by re-mining CRC microarray. Functional assays were performed in vitro and in vivo. RNA-Seq, RNA pull-down, and RIP assays were used to investigate the mechanisms of Lnc-FAM84B-4 in regulating expression of DUSP1. The results indicated that Lnc-FAM84B-4 regulates MAPK pathway by restraining DUSP1 expression. Mechanistically, RNA pull-down followed by mass spectrum determined hnRNPK functions as a binding partner of lnc-FAM84B-4 in mediating DUSP1 expression. Our findings demonstrate the important role of lnc-FAM84B-4-hnRNPK-DUSP1 axis in CRC development, and suggest a therapeutic target for CRC treatment.

METTL6 is a tRNA m3C methyltransferase that regulates pluripotency and tumor cell growth

Recently, covalent modifications of RNA, such as methylation, have emerged as key regulators of all aspects of RNA biology and have been implicated in numerous diseases, for instance, cancer. Here, we undertook a combination of in vitro and in vivo screens to test 78 potential methyltransferases for their roles in hepatocellular carcinoma (HCC) cell proliferation. We identified methyltransferase-like protein 6 (METTL6) as a crucial regulator of tumor cell growth. We show that METTL6 is a bona fide transfer RNA (tRNA) methyltransferase, catalyzing the formation of 3-methylcytidine at C32 of specific serine tRNA isoacceptors. Deletion of Mettl6 in mouse stem cells results in changes in ribosome occupancy and RNA levels, as well as impaired pluripotency. In mice, Mettl6 knockout results in reduced energy expenditure. We reveal a previously unknown pathway in the maintenance of translation efficiency with a role in maintaining stem cell self-renewal, as well as impacting tumor cell growth profoundly.

FMRP Control of Ribosome Translocation Promotes Chromatin Modifications and Alternative Splicing of Neuronal Genes Linked to Autism

Silencing of FMR1 and loss of its gene product, FMRP, results in fragile X syndrome (FXS). FMRP binds brain mRNAs and inhibits polypeptide elongation. Using ribosome profiling of the hippocampus, we find that ribosome footprint levels in Fmr1-deficient tissue mostly reflect changes in RNA abundance. Profiling over a time course of ribosome runoff in wild-type tissue reveals a wide range of ribosome translocation rates; on many mRNAs, the ribosomes are stalled. Sucrose gradient ultracentrifugation of hippocampal slices after ribosome runoff reveals that FMRP co-sediments with stalled ribosomes, and its loss results in decline of ribosome stalling on specific mRNAs. One such mRNA encodes SETD2, a lysine methyltransferase that catalyzes H3K36me3. Chromatin immunoprecipitation sequencing (ChIP-seq) demonstrates that loss of FMRP alters the deployment of this histone mark. H3K36me3 is associated with alternative pre-RNA processing, which we find occurs in an FMRP-dependent manner on transcripts linked to neural function and autism spectrum disorders.

Control of noncoding RNA production and histone levels by a 5' tRNA fragment

In this study Boskovic et al. set out to elucidate the functions of a small RNA derived from the 5’ end of mature tRNA-Gly-GCC. Using several genomic, biochemical, and molecular methods, the authors reveal a conserved mechanism for 5’ tRNA fragment control of noncoding RNA biogenesis and global chromatin organization.

Small RNAs derived from mature tRNAs, referred to as tRNA fragments or “tRFs,” are an emerging class of regulatory RNAs with poorly understood functions. We recently identified a role for one specific tRF—5′ tRF-Gly-GCC, or tRF-GG—as a repressor of genes associated with the endogenous retroelement MERVL, but the mechanistic basis for this regulation was unknown. Here, we show that tRF-GG plays a role in production of a wide variety of noncoding RNAs—snoRNAs, scaRNAs, and snRNAs—that are dependent on Cajal bodies for stability and activity. Among these noncoding RNAs, regulation of the U7 snRNA by tRF-GG modulates heterochromatin-mediated transcriptional repression of MERVL elements by supporting an adequate supply of histone proteins. Importantly, the effects of inhibiting tRF-GG on histone mRNA levels, on activity of a histone 3′ UTR reporter, and ultimately on MERVL regulation could all be suppressed by manipulating U7 RNA levels. We additionally show that the related RNA-binding proteins hnRNPF and hnRNPH bind directly to tRF-GG, and are required for Cajal body biogenesis, positioning these proteins as strong candidates for effectors of tRF-GG function in vivo. Together, our data reveal a conserved mechanism for 5′ tRNA fragment control of noncoding RNA biogenesis and, consequently, global chromatin organization.

Antibody cross-reactivity accounts for widespread appearance of m1A in 5'UTRs

N1-methyladenosine (m1A) was proposed to be a highly prevalent modification in mRNA 5'UTRs based on mapping studies using an m1A-binding antibody. We developed a bioinformatic approach to discover m1A and other modifications in mRNA throughout the transcriptome by analyzing preexisting ultra-deep RNA-Seq data for modification-induced misincorporations. Using this approach, we detected appreciable levels of m1A only in one mRNA: the mitochondrial MT-ND5 transcript. As an alternative approach, we also developed an antibody-based m1A-mapping approach to detect m1A at single-nucleotide resolution, and confirmed that the commonly used m1A antibody maps sites to the transcription-start site in mRNA 5'UTRs. However, further analysis revealed that these were false-positives caused by binding of the antibody to the m7G-cap. A different m1A antibody that lacks cap-binding cross-reactivity does not show enriched binding in 5'UTRs. These results demonstrate that high-stoichiometry m1A sites are exceedingly rare in mRNAs and that previous mappings of m1A to 5'UTRs were the result of antibody cross-reactivity to the 5' cap.

Free circular introns with an unusual branchpoint in neuronal projections

The polarized structure of axons and dendrites in neuronal cells depends in part on RNA localization. Previous studies have looked at which polyadenylated RNAs are enriched in neuronal projections or at synapses, but less is known about the distribution of non-adenylated RNAs. By physically dissecting projections from cell bodies of primary rat hippocampal neurons and sequencing total RNA, we found an unexpected set of free circular introns with a non-canonical branchpoint enriched in neuronal projections. These introns appear to be tailless lariats that escape debranching. They lack ribosome occupancy, sequence conservation, and known localization signals, and their function, if any, is not known. Nonetheless, their enrichment in projections has important implications for our understanding of the mechanisms by which RNAs reach distal compartments of asymmetric cells.

Staphylococcus aureus Small Colony Variants (SCVs): News From a Chronic Prosthetic Joint Infection

Small colony variants (SCV) of Staphylococcus aureus have been reported as implicated in chronic infections. Here, we investigated the genomic and transcriptomic changes involved in the evolution from a wild-type to a SCV from in a patient with prosthetic joint infection relapse. The SCV presented a stable phenotype with no classical auxotrophy and the emergence of rifampicin resistance. Whole Genome Sequencing (WGS) analysis showed only the loss of a 42.5 kb phage and 3 deletions, among which one targeting the rpoB gene, known to be the target of rifampicin and to be associated to SCV formation in the context of a constitutively active stringent response. Transcriptomic analysis highlighted a specific signature in the SCV strain including a complex, multi-level strategy of survival and adaptation to chronicity within the host including a protection from the inflammatory response, an evasion of the immune response, a constitutively activated stringent response and a scavenging of iron sources.

Senataxin homologue Sen1 is required for efficient termination of RNA polymerase III transcription

R-loop disassembly by the human helicase Senataxin contributes to genome integrity and to proper transcription termination at a subset of RNA polymerase II genes. Whether Senataxin also contributes to transcription termination at other classes of genes has remained unclear. Here, we show that Sen1, one of two fission yeast homologues of Senataxin, promotes efficient termination of RNA polymerase III (RNAP3) transcription in vivo. In the absence of Sen1, RNAP3 accumulates downstream of RNAP3-transcribed genes and produces long exosome-sensitive 3'-extended transcripts. Importantly, neither of these defects was affected by the removal of R-loops. The finding that Sen1 acts as an ancillary factor for RNAP3 transcription termination in vivo challenges the pre-existing view that RNAP3 terminates transcription autonomously. We propose that Sen1 is a cofactor for transcription termination that has been co-opted by different RNA polymerases in the course of evolution.

Identification of Footprints of RNA:Protein Complexes via RNA Immunoprecipitation in Tandem Followed by Sequencing (RIPiT-Seq)

RNA immunoprecipitation in tandem (RIPiT) is a method for enriching RNA footprints of a pair of proteins within an RNA:protein (RNP) complex. RIPiT employs two purification steps. First, immunoprecipitation of a tagged RNP subunit is followed by mild RNase digestion and subsequent non-denaturing affinity elution. A second immunoprecipitation of another RNP subunit allows for enrichment of a defined complex. Following a denaturing elution of RNAs and proteins, the RNA footprints are converted into high-throughput DNA sequencing libraries. Unlike the more popular ultraviolet (UV) crosslinking followed by immunoprecipitation (CLIP) approach to enrich RBP binding sites, RIPiT is UV-crosslinking independent. Hence RIPiT can be applied to numerous proteins present in the RNA interactome and beyond that are essential to RNA regulation but do not directly contact the RNA or UV-crosslink poorly to RNA. The two purification steps in RIPiT provide an additional advantage of identifying binding sites where a protein of interest acts in partnership with another cofactor. The double purification strategy also serves to enhance signal by limiting background. Here, we provide a step-wise procedure to perform RIPiT and to generate high-throughput sequencing libraries from isolated RNA footprints. We also outline RIPiT's advantages and applications and discuss some of its limitations.

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

Thermostable group II intron reverse transcriptases (TGIRTs) with high fidelity and processivity have been used for a variety of RNA sequencing (RNA-seq) applications, including comprehensive profiling of whole-cell, exosomal, and human plasma RNAs; quantitative tRNA-seq based on the ability of TGIRT enzymes to give full-length reads of tRNAs and other structured small ncRNAs; high-throughput mapping of post-transcriptional modifications; and RNA structure mapping. Here, we improved TGIRT-seq methods for comprehensive transcriptome profiling by rationally designing RNA-seq adapters that minimize adapter dimer formation. Additionally, we developed biochemical and computational methods for remediating 5′- and 3′-end biases, the latter based on a random forest regression model that provides insight into the contribution of different factors to these biases. These improvements, some of which may be applicable to other RNA-seq methods, increase the efficiency of TGIRT-seq library construction and improve coverage of very small RNAs, such as miRNAs. Our findings provide insight into the biochemical basis of 5′- and 3′-end biases in RNA-seq and suggest general approaches for remediating biases and decreasing adapter dimer formation.

Optimization of ribosome profiling using low-input brain tissue from fragile X syndrome model mice

Dysregulated protein synthesis is a major underlying cause of many neurodevelopmental diseases including fragile X syndrome. In order to capture subtle but biologically significant differences in translation in these disorders, a robust technique is required. One powerful tool to study translational control is ribosome profiling, which is based on deep sequencing of mRNA fragments protected from ribonuclease (RNase) digestion by ribosomes. However, this approach has been mainly applied to rapidly dividing cells where translation is active and large amounts of starting material are readily available. The application of ribosome profiling to low-input brain tissue where translation is modest and gene expression changes between genotypes are expected to be small has not been carefully evaluated. Using hippocampal tissue from wide type and fragile X mental retardation 1 (Fmr1) knockout mice, we show that variable RNase digestion can lead to significant sample batch effects. We also establish GC content and ribosome footprint length as quality control metrics for RNase digestion. We performed RNase titration experiments for low-input samples to identify optimal conditions for this critical step that is often improperly conducted. Our data reveal that optimal RNase digestion is essential to ensure high quality and reproducibility of ribosome profiling for low-input brain tissue.

Genome editing in primary cells and in vivo using viral-derived Nanoblades loaded with Cas9-sgRNA ribonucleoproteins

Programmable nucleases have enabled rapid and accessible genome engineering in eukaryotic cells and living organisms. However, their delivery into target cells can be technically challenging when working with primary cells or in vivo. Here, we use engineered murine leukemia virus-like particles loaded with Cas9-sgRNA ribonucleoproteins (Nanoblades) to induce efficient genome-editing in cell lines and primary cells including human induced pluripotent stem cells, human hematopoietic stem cells and mouse bone-marrow cells. Transgene-free Nanoblades are also capable of in vivo genome-editing in mouse embryos and in the liver of injected mice. Nanoblades can be complexed with donor DNA for "all-in-one" homology-directed repair or programmed with modified Cas9 variants to mediate transcriptional up-regulation of target genes. Nanoblades preparation process is simple, relatively inexpensive and can be easily implemented in any laboratory equipped for cellular biology.

Regulatory discrimination of mRNAs by FMRP controls mouse adult neural stem cell differentiation

Fragile X syndrome (FXS) is caused by the loss of fragile X mental retardation protein (FMRP), an RNA binding protein whose deficiency impacts many brain functions, including differentiation of adult neural stem cells (aNSCs). However, the mechanism by which FMRP influences these processes remains unclear. Here, we performed ribosome profiling and transcriptomic analysis of aNSCs in parallel from wild-type and Fmr1 knockout mice. Our data revealed diverse gene expression changes at both mRNA and translation levels. Many mitosis and neurogenesis genes were dysregulated primarily at the mRNA level, while numerous synaptic genes were mostly dysregulated at the translation level. Translational "buffering", whereby changes in ribosome association with mRNA are compensated by alterations in RNA abundance, was also evident. Knockdown of NECDIN, an FMRP-repressed transcriptional factor, rescued neuronal differentiation. In addition, we discovered that FMRP regulates mitochondrial mRNA expression and energy homeostasis. Thus, FMRP controls diverse transcriptional and posttranscriptional gene expression programs critical for neural differentiation.

Decreasing miRNA sequencing bias using a single adapter and circularization approach

The ability to accurately quantify all the microRNAs (miRNAs) in a sample is important for understanding miRNA biology and for development of new biomarkers and therapeutic targets. We develop a new method for preparing miRNA sequencing libraries, RealSeq®-AC, that involves ligating the miRNAs with a single adapter and circularizing the ligation products. When compared to other methods, RealSeq®-AC provides greatly reduced miRNA sequencing bias and allows the identification of the largest variety of miRNAs in biological samples. This reduced bias also allows robust quantification of miRNAs present in samples across a wide range of RNA input levels.

Transcriptome-wide Interrogation of the Functional Intronome by Spliceosome Profiling

Full understanding of eukaryotic transcriptomes and how they respond to different conditions requires deep knowledge of all sites of intron excision. Although RNA sequencing (RNA-seq) provides much of this information, the low abundance of many spliced transcripts (often due to their rapid cytoplasmic decay) limits the ability of RNA-seq alone to reveal the full repertoire of spliced species. Here, we present "spliceosome profiling," a strategy based on deep sequencing of RNAs co-purifying with late-stage spliceosomes. Spliceosome profiling allows for unambiguous mapping of intron ends to single-nucleotide resolution and branchpoint identification at unprecedented depths. Our data reveal hundreds of new introns in S. pombe and numerous others that were previously misannotated. By providing a means to directly interrogate sites of spliceosome assembly and catalysis genome-wide, spliceosome profiling promises to transform our understanding of RNA processing in the nucleus, much as ribosome profiling has transformed our understanding mRNA translation in the cytoplasm.

Cloning and Identification of Recombinant Argonaute-Bound Small RNAs Using Next-Generation Sequencing

Argonaute proteins (AGOs) are loaded with small RNAs as guides to recognize target mRNAs. Since the target specificity heavily depends on the base complementarity between two strands, it is important to identify small guide and long target RNAs bound to AGOs. For this purpose, next-generation sequencing (NGS) technologies have extended our appreciation truly to the nucleotide level. However, the identification of RNAs via NGS from scarce RNA samples remains a challenge. Further, most commercial and published methods are compatible with either small RNAs or long RNAs, but are not equally applicable to both. Therefore, a single method that yields quantitative, bias-free NGS libraries to identify small and long RNAs from low levels of input will be of wide interest. Here, we introduce such a procedure that is based on several modifications of two published protocols and allows robust, sensitive, and reproducible cloning and sequencing of small amounts of RNAs of variable lengths. The method was applied to the identification of small RNAs bound to a purified eukaryotic AGO. Following ligation of a DNA adapter to RNA 3'-end, the key feature of this method is to use the adapter for priming reverse transcription (RT) wherein biotinylated deoxyribonucleotides specifically incorporated into the extended complementary DNA. Such RT products are enriched on streptavidin beads, circularized while immobilized on beads and directly used for PCR amplification. We provide a stepwise guide to generate RNA-Seq libraries, their purification, quantification, validation, and preparation for next-generation sequencing. We also provide basic steps in post-NGS data analyses using Galaxy, an open-source, web-based platform.

Large Scale Mutational and Kinetic Analysis of a Self-Hydrolyzing Deoxyribozyme

Deoxyribozymes are catalytic DNA sequences whose atomic structures are generally difficult to elucidate. Mutational analysis remains a principal approach for understanding and engineering deoxyribozymes with diverse catalytic activities. However, laborious preparation and biochemical characterization of individual sequences severely limit the number of mutants that can be studied biochemically. Here, we applied deep sequencing to directly measure the activities of self-hydrolyzing deoxyribozyme sequences in high throughput. First, all single and double mutants within the 15-base catalytic core of the deoxyribozyme I-R3 were assayed to unambiguously determine the tolerated and untolerated mutations at each position. Subsequently, 4096 deoxyribozyme variants with tolerated base substitutions at seven positions were kinetically assayed in parallel. We identified 533 active mutants whose first-order rate constants and activation energies were determined. The results indicate an isolated and narrow peak in the deoxyribozyme sequence space and provide a quantitative view of the effects of multiple mutations on the deoxyribozyme activity for the first time.

Transcriptome-wide Analysis of Roles for tRNA Modifications in Translational Regulation

Covalent nucleotide modifications in noncoding RNAs affect a plethora of biological processes, and new functions continue to be discovered even for well-known modifying enzymes. To systematically compare the functions of a large set of noncoding RNA modifications in gene regulation, we carried out ribosome profiling in budding yeast to characterize 57 nonessential genes involved in tRNA modification. Deletion mutants exhibited a range of translational phenotypes, with enzymes known to modify anticodons, or non-tRNA substrates such as rRNA, exhibiting the most dramatic translational perturbations. Our data build on prior reports documenting translational upregulation of the nutrient-responsive transcription factor Gcn4 in response to numerous tRNA perturbations, and identify many additional translationally regulated mRNAs throughout the yeast genome. Our data also uncover unexpected roles for tRNA-modifying enzymes in regulation of TY retroelements, and in rRNA 2'-O-methylation. This dataset should provide a rich resource for discovery of additional links between tRNA modifications and gene regulation.

FLT1 and transcriptome-wide polyadenylation site (PAS) analysis in preeclampsia

Maternal symptoms of preeclampsia (PE) are primarily driven by excess anti-angiogenic factors originating from the placenta. Chief among these are soluble Flt1 proteins (sFlt1s) produced from alternatively polyadenylated mRNA isoforms. Here we used polyadenylation site sequencing (PAS-Seq) of RNA from normal and PE human placentae to interrogate transcriptome-wide gene expression and alternative polyadenylation signatures associated with early-onset PE (EO-PE; symptom onset < 34 weeks) and late-onset PE (LO-PE; symptom onset > 34 weeks) cohorts. While we observed no general shift in alternative polyadenylation associated with PE, the EO-PE and LO-PE cohorts do exhibit gene expression profiles distinct from both each other and from normal placentae. The only two genes upregulated across all transcriptome-wide PE analyses to date (microarray, RNA-Seq and PAS-Seq) are NRIP1 (RIP140), a transcriptional co-regulator linked to metabolic syndromes associated with obesity, and Flt1. Consistent with sFlt1 overproduction being a significant driver of clinical symptoms, placental Flt1 mRNA levels strongly correlate with maternal blood pressure. For Flt1, just three mRNA isoforms account for > 94% of all transcripts, with increased transcription of the entire locus driving Flt1 upregulation in both EO-PE and LO-PE. These three isoforms thus represent potential targets for therapeutic RNA interference (RNAi) in both early and late presentations.

Transcriptome-wide measurement of translation by ribosome profiling

Translation is one of the fundamental processes of life. It comprises the assembly of polypeptides whose amino acid sequence corresponds to the codon sequence of an mRNA's ORF. Translation is performed by the ribosome; therefore, in order to understand translation and its regulation we must be able to determine the numbers and locations of ribosomes on mRNAs in vivo. Furthermore, we must be able to examine their redistribution in different physiological contexts and in response to experimental manipulations. The ribosome profiling method provides us with an opportunity to learn these locations, by sequencing a cDNA library derived from the short fragments of mRNA covered by the ribosome. Since its original description, the ribosome profiling method has undergone continuing development; in this article we describe the method's current state. Important improvements include: the incorporation of sample barcodes to enable library multiplexing, the incorporation of unique molecular identifiers to enable to removal of duplicated sequences, and the replacement of a gel-purification step with the enzymatic degradation of unligated linker.

Reversible methylation of m6Am in the 5' cap controls mRNA stability

Internal bases in mRNA can be subjected to modifications that influence the fate of mRNA in cells. One of the most prevalent modified bases is found at the 5' end of mRNA, at the first encoded nucleotide adjacent to the 7-methylguanosine cap. Here we show that this nucleotide, N6,2'-O-dimethyladenosine (m6Am), is a reversible modification that influences cellular mRNA fate. Using a transcriptome-wide map of m6Am we find that m6Am-initiated transcripts are markedly more stable than mRNAs that begin with other nucleotides. We show that the enhanced stability of m6Am-initiated transcripts is due to resistance to the mRNA-decapping enzyme DCP2. Moreover, we find that m6Am is selectively demethylated by fat mass and obesity-associated protein (FTO). FTO preferentially demethylates m6Am rather than N6-methyladenosine (m6A), and reduces the stability of m6Am mRNAs. Together, these findings show that the methylation status of m6Am in the 5' cap is a dynamic and reversible epitranscriptomic modification that determines mRNA stability.

Capture and sequencing of NAD-capped RNA sequences with NAD captureSeq

Here we describe a protocol for NAD captureSeq that allows for the identification of nicotinamide-adenine dinucleotide (NAD)-capped RNA sequences in total RNA samples from different organisms. NAD-capped RNA is first chemo-enzymatically biotinylated with high efficiency, permitting selective capture on streptavidin beads. Then, a highly efficient library preparation protocol tailored to immobilized, 5'-modified RNA is applied, with adaptor ligation to the RNA's 3' terminus and reverse transcription (RT) performed on-bead. Then, cDNA is released into solution, tailed, ligated to a second adaptor and PCR-amplified. After next-generation sequencing (NGS) of the DNA library, enriched sequences are identified by comparison with a control sample in which the first step of chemo-enzymatic biotinylation is omitted. Because the downstream protocol does not necessarily rely on NAD-modified but on 'clickable' or biotin-modified RNA, it can be applied to other RNA modifications or RNA-biomolecule interactions. The central part of this protocol can be completed in ∼7 d, excluding preparatory steps, sequencing and bioinformatic analysis.

A Long Noncoding RNA lincRNA-EPS Acts as a Transcriptional Brake to Restrain Inflammation

Long intergenic noncoding RNAs (lincRNAs) are important regulators of gene expression. Although lincRNAs are expressed in immune cells, their functions in immunity are largely unexplored. Here, we identify an immunoregulatory lincRNA, lincRNA-EPS, that is precisely regulated in macrophages to control the expression of immune response genes (IRGs). Transcriptome analysis of macrophages from lincRNA-EPS-deficient mice, combined with gain-of-function and rescue experiments, revealed a specific role for this lincRNA in restraining IRG expression. Consistently, lincRNA-EPS-deficient mice manifest enhanced inflammation and lethality following endotoxin challenge in vivo. lincRNA-EPS localizes at regulatory regions of IRGs to control nucleosome positioning and repress transcription. Further, lincRNA-EPS mediates these effects by interacting with heterogeneous nuclear ribonucleoprotein L via a CANACA motif located in its 3' end. Together, these findings identify lincRNA-EPS as a repressor of inflammatory responses, highlighting the importance of lincRNAs in the immune system.

Dual randomization of oligonucleotides to reduce the bias in ribosome-profiling libraries

Protein translation is at the heart of cellular metabolism and its in-depth characterization is key for many lines of research. Recently, ribosome profiling became the state-of-the-art method to quantitatively characterize translation dynamics at a transcriptome-wide level. However, the strategy of library generation affects its outcomes. Here, we present a modified ribosome-profiling protocol starting from yeast, human cells and vertebrate brain tissue. We use a DNA linker carrying four randomized positions at its 5′ end and a reverse-transcription (RT) primer with three randomized positions to reduce artifacts during library preparation. The use of seven randomized nucleotides allows to efficiently detect library-generation artifacts. We find that the effect of polymerase chain reaction (PCR) artifacts is relatively small for global analyses when sufficient input material is used. However, when input material is limiting, our strategy improves the sensitivity of gene-specific analyses. Furthermore, randomized nucleotides alleviate the skewed frequency of specific sequences at the 3′ end of ribosome-protected fragments (RPFs) likely resulting from ligase specificity. Finally, strategies that rely on dual ligation show a high degree of gene-coverage variation. Taken together, our approach helps to remedy two of the main problems associated with ribosome-profiling data. This will facilitate the analysis of translational dynamics and increase our understanding of the influence of RNA modifications on translation.

Redefining the Translational Status of 80S Monosomes

Fully assembled ribosomes exist in two populations: polysomes and monosomes. While the former has been studied extensively, to what extent translation occurs on monosomes and its importance for overall translational output remain controversial. Here, we used ribosome profiling to examine the translational status of 80S monosomes in Saccharomyces cerevisiae. We found that the vast majority of 80S monosomes are elongating, not initiating. Further, most mRNAs exhibit some degree of monosome occupancy, with monosomes predominating on nonsense-mediated decay (NMD) targets, upstream open reading frames (uORFs), canonical ORFs shorter than ∼ 590 nt, and ORFs for which the total time required to complete elongation is substantially shorter than that required for initiation. Importantly, mRNAs encoding low-abundance regulatory proteins tend to be enriched in the monosome fraction. Our data highlight the importance of monosomes for the translation of highly regulated mRNAs.

RNA-Seq methods for transcriptome analysis

Deep sequencing has been revolutionizing biology and medicine in recent years, providing single base-level precision for our understanding of nucleic acid sequences in high throughput fashion. Sequencing of RNA, or RNA-Seq, is now a common method to analyze gene expression and to uncover novel RNA species. Aspects of RNA biogenesis and metabolism can be interrogated with specialized methods for cDNA library preparation. In this study, we review current RNA-Seq methods for general analysis of gene expression and several specific applications, including isoform and gene fusion detection, digital gene expression profiling, targeted sequencing and single-cell analysis. In addition, we discuss approaches to examine aspects of RNA in the cell, technical challenges of existing RNA-Seq methods, and future directions. WIREs RNA 2017, 8:e1364. doi: 10.1002/wrna.1364 For further resources related to this article, please visit the WIREs website.

Methods to Study Long Noncoding RNA Biology in Cancer

Thousands of long noncoding RNAs (lncRNAs) have been discovered in recent years. The functions of lncRNAs range broadly from regulating chromatin structure and gene expression in the nucleus to controlling messenger RNA (mRNA) processing, mRNA posttranscriptional regulation, cellular signaling, and protein activity in the cytoplasm. Experimental and computational techniques have been developed to characterize lncRNAs in high-throughput scale, to study the lncRNA function in vitro and in vivo, to map lncRNA binding sites on the genome, and to capture lncRNA-protein interactions with the identification of lncRNA-binding partners, binding sites, and interaction determinants. In this chapter, we will discuss these technologies and their applications in decoding the functions of lncRNAs. Understanding these techniques including their advantages and disadvantages and developing them in the future will be essential to elaborate the roles of lncRNAs in cancer and other diseases.