Plastid: nucleotide-resolution analysis of next-generation sequencing and genomics data

Bigtools: a high-performance BigWig and BigBed library in Rust

MOTIVATION

The BigWig and BigBed file formats were originally designed for the visualization of next-generation sequencing data through a genome browser. Due to their versatility, these formats have long since become ubiquitous for the storage of processed sequencing data and regularly serve as the basis for downstream data analysis. As the number and size of sequencing experiments continues to accelerate, there is an increasing demand to efficiently generate and query BigWig and BigBed files in a scalable and robust manner, and to efficiently integrate these functionalities into data analysis environments and third-party applications.

RESULTS

Here, we present Bigtools, a feature-complete, high-performance, and integrable software library for generating and querying both BigWig and BigBed files. Bigtools is written in the Rust programming language and includes a flexible suite of command line tools as well as bindings to Python.

AVAILABILITY AND IMPLEMENTATION

Bigtools is cross-platform and released under the MIT license. It is distributed on Crates.io, Bioconda, and the Python Package Index, and the source code is available at https://github.com/jackh726/bigtools.

The ABCF ATPase New1 resolves translation termination defects associated with specific tRNAArg and tRNALys isoacceptors in the P site

The efficiency of translation termination is determined by the nature of the stop codon as well as its context. In eukaryotes, recognition of the A-site stop codon and release of the polypeptide are mediated by release factors eRF1 and eRF3, respectively. Translation termination is modulated by other factors which either directly interact with release factors or bind to the E-site and modulate the activity of the peptidyl transferase center. Previous studies suggested that the Saccharomyces cerevisiae ABCF ATPase New1 is involved in translation termination and/or ribosome recycling, however, the exact function remained unclear. Here, we have applied 5PSeq, single-particle cryo-EM and readthrough reporter assays to provide insight into the biological function of New1. We show that the lack of New1 results in ribosomal stalling at stop codons preceded by a lysine or arginine codon and that the stalling is not defined by the nature of the C-terminal amino acid but rather by the identity of the tRNA isoacceptor in the P-site. Collectively, our results suggest that translation termination is inefficient when ribosomes have specific tRNA isoacceptors in the P-site and that the recruitment of New1 rescues ribosomes at these problematic termination contexts.

Annotation of nuclear lncRNAs based on chromatin interactions

The human genome is pervasively transcribed and produces a wide variety of long non-coding RNAs (lncRNAs), constituting the majority of transcripts across human cell types. Some specific nuclear lncRNAs have been shown to be important regulatory components acting locally. As RNA-chromatin interaction and Hi-C chromatin conformation data showed that chromatin interactions of nuclear lncRNAs are determined by the local chromatin 3D conformation, we used Hi-C data to identify potential target genes of lncRNAs. RNA-protein interaction data suggested that nuclear lncRNAs act as scaffolds to recruit regulatory proteins to target promoters and enhancers. Nuclear lncRNAs may therefore play a role in directing regulatory factors to locations spatially close to the lncRNA gene. We provide the analysis results through an interactive visualization web portal at https://fantom.gsc.riken.jp/zenbu/reports/#F6_3D_lncRNA.

Proteins à la carte: riboproteogenomic exploration of bacterial N-terminal proteoform expression

In recent years, it has become evident that the true complexity of bacterial proteomes remains underestimated. Gene annotation tools are known to propagate biases and overlook certain classes of truly expressed proteins, particularly proteoforms-protein isoforms arising from a single gene. Recent (re-)annotation efforts heavily rely on ribosome profiling by providing a direct readout of translation to fully describe bacterial proteomes. In this study, we employ a robust riboproteogenomic pipeline to conduct a systematic census of expressed N-terminal proteoform pairs, representing two isoforms encoded by a single gene raised by annotated and alternative translation initiation, in Salmonella. Intriguingly, conditional-dependent changes in relative utilization of annotated and alternative translation initiation sites (TIS) were observed in several cases. This suggests that TIS selection is subject to regulatory control, adding yet another layer of complexity to our understanding of bacterial proteomes.

IMPORTANCE

With the emerging theme of genes within genes comprising the existence of alternative open reading frames (ORFs) generated by translation initiation at in-frame start codons, mechanisms that control the relative utilization of annotated and alternative TIS need to be unraveled and our molecular understanding of resulting proteoforms broadened. Utilizing complementary ribosome profiling strategies to map ORF boundaries, we uncovered dual-encoding ORFs generated by in-frame TIS usage in Salmonella. Besides demonstrating that alternative TIS usage may generate proteoforms with different characteristics, such as differential localization and specialized function, quantitative aspects of conditional retapamulin-assisted ribosome profiling (Ribo-RET) translation initiation maps offer unprecedented insights into the relative utilization of annotated and alternative TIS, enabling the exploration of gene regulatory mechanisms that control TIS usage and, consequently, the translation of N-terminal proteoform pairs.

Deep learning to decode sites of RNA translation in normal and cancerous tissues

The biological process of RNA translation is fundamental to cellular life and has wide-ranging implications for human disease. Yet, accurately delineating the variation in RNA translation represents a significant challenge. Here, we develop RiboTIE, a transformer model-based approach to map global RNA translation. We find that RiboTIE offers unparalleled precision and sensitivity for ribosome profiling data. Application of RiboTIE to normal brain and medulloblastoma cancer samples enables high-resolution insights into disease regulation of RNA translation.

No country for old methods: New tools for studying microproteins

Microproteins encoded by small open reading frames (sORFs) have emerged as a fascinating frontier in genomics. Traditionally overlooked due to their small size, recent technological advancements such as ribosome profiling, mass spectrometry-based strategies and advanced computational approaches have led to the annotation of more than 7000 sORFs in the human genome. Despite the vast progress, only a tiny portion of these microproteins have been characterized and an important challenge in the field lies in identifying functionally relevant microproteins and understanding their role in different cellular contexts. In this review, we explore the recent advancements in sORF research, focusing on the new methodologies and computational approaches that have facilitated their identification and functional characterization. Leveraging these new tools hold great promise for dissecting the diverse cellular roles of microproteins and will ultimately pave the way for understanding their role in the pathogenesis of diseases and identifying new therapeutic targets.

Bigtools: a high-performance BigWig and BigBed library in Rust

The BigWig and BigBed file formats were originally designed for the visualization of next-generation sequencing data through a genome browser. Due to their versatility, these formats have long since become ubiquitous for the storage of processed sequencing data and regularly serve as the basis for downstream data analysis. As the number and size of sequencing experiments continues to accelerate, there is an increasing demand to efficiently generate and query BigWig and BigBed files in a scalable and robust manner, and to efficiently integrate these functionalities into data analysis environments and third-party applications. Here, we present Bigtools, a feature-complete, high-performance, and integrable software library for generating and querying both BigWig and BigBed files. Bigtools is written in the Rust programming language and includes a flexible suite of command line tools as well as bindings to Python. Bigtools is cross-platform and released under the MIT license. It is distributed on Crates.io and the Python Package Index, and the source code is available at https://github.com/jackh726/bigtools.

Y-Box-Binding Proteins Have a Dual Impact on Cellular Translation

Y-box-binding proteins (YB proteins) are multifunctional DNA- and RNA-binding proteins that play an important role in the regulation of gene expression. The high homology of their cold shock domains and the similarity between their long, unstructured C-terminal domains suggest that Y-box-binding proteins may have similar functions in a cell. Here, we consider the functional interchangeability of the somatic YB proteins YB-1 and YB-3. RNA-seq and Ribo-seq are used to track changes in the mRNA abundance or mRNA translation in HEK293T cells solely expressing YB-1, YB-3, or neither of them. We show that YB proteins have a dual effect on translation. Although the expression of YB proteins stimulates global translation, YB-1 and YB-3 inhibit the translation of their direct CLIP-identified mRNA targets. The impact of YB-1 and YB-3 on the translation of their mRNA targets is similar, which suggests that they can substitute each other in inhibiting the translation of their mRNA targets in HEK293T cells.

Deciphering the Functional Long Non-Coding RNAs Derived from MicroRNA Loci

Albeit the majority of eukaryotic genomes can be pervasively transcribed to a diverse population of lncRNAs and various subtypes of lncRNA are discovered. However, the genome-wide study of miRNA-derived lncRNAs is still lacking. Here, it is reported that over 800 miRNA gene-originated lncRNAs (molncRNAs) are generated from miRNA loci. One of them, molnc-301b from miR-301b and miR-130b, functions as an "RNA decoy" to facilitate dissociation of the chromatin remodeling protein SMARCA5 from chromatin and thereby sequester transcription and mRNA translation. Specifically, molnc-301b attenuates erythropoiesis by mitigating the transcription of erythropoietic and translation-associated genes, such as GATA1 and FOS. In addition, a useful and powerful CRISPR screen platform to characterize the biological functions of molncRNAs at large-scale and single-cell levels is established and 29 functional molncRNAs in hematopoietic cells are identified. Collectively, the focus is on miRNA-derived lncRNAs, deciphering their landscape during normal hematopoiesis, and comprehensively evaluating their potential roles.

A PARP14/TARG1-Regulated RACK1 MARylation Cycle Drives Stress Granule Dynamics in Ovarian Cancer Cells

Mono(ADP-ribosyl)ation (MARylation), a post-translational modification (PTM) of proteins, is emerging as a critical regulator of ribosome function and translation. Herein, we demonstrate that RACK1, a member of the tryptophan-aspartate repeat (WD-repeat) family of proteins and an integral component of the ribosome, is MARylated by the mono(ADP-ribosyl) transferase (MART) PARP14 in ovarian cancer cells. We mapped and confirmed the sites of MARylation, which occur on three acidic residues within blades 4 and 5 of β-propeller domain of RACK1, a chaperone that shuttles and anchors proteins where needed. Site-specific MARylation of RACK1 is required for stress granule formation and promotes the colocalization of RACK1 to stress granules with key components, such as G3BP1, eIF3η, and 40S ribosomal proteins. In parallel, we observed reduced translation of a subset of mRNAs, including those encoding key cancer regulators (e.g., AKT). Treatment with a PARP14 inhibitor or mutation of the sites of MARylation on RACK1 blocks these outcomes. To re-set the system after prolonged stress and recovery, the ADP-ribosyl hydrolase TARG1 deMARylates RACK1 to dissociate the stress granules and return RACK1 and the 40S ribosomal subunit to the cytoplasm, allowing for a restoration of translation. Collectively, our results highlight the discovery of a PARP14/TARG1-regulated RACK1 MARylation cycle that controls stress granule assembly and disassembly in ovarian cancer cells.

Pervasive downstream RNA hairpins dynamically dictate start-codon selection

Translational reprogramming allows organisms to adapt to changing conditions. Upstream start codons (uAUGs), which are prevalently present in mRNAs, have crucial roles in regulating translation by providing alternative translation start sites1-4. However, what determines this selective initiation of translation between conditions remains unclear. Here, by integrating transcriptome-wide translational and structural analyses during pattern-triggered immunity in Arabidopsis, we found that transcripts with immune-induced translation are enriched with upstream open reading frames (uORFs). Without infection, these uORFs are selectively translated owing to hairpins immediately downstream of uAUGs, presumably by slowing and engaging the scanning preinitiation complex. Modelling using deep learning provides unbiased support for these recognizable double-stranded RNA structures downstream of uAUGs (which we term uAUG-ds) being responsible for the selective translation of uAUGs, and allows the prediction and rational design of translating uAUG-ds. We found that uAUG-ds-mediated regulation can be generalized to human cells. Moreover, uAUG-ds-mediated start-codon selection is dynamically regulated. After immune challenge in plants, induced RNA helicases that are homologous to Ded1p in yeast and DDX3X in humans resolve these structures, allowing ribosomes to bypass uAUGs to translate downstream defence proteins. This study shows that mRNA structures dynamically regulate start-codon selection. The prevalence of this RNA structural feature and the conservation of RNA helicases across kingdoms suggest that mRNA structural remodelling is a general feature of translational reprogramming.

Extension of mRNA poly(A) tails and 3'UTRs during neuronal differentiation exhibits variable association with post-transcriptional dynamics

Differentiation of neural progenitor cells into mature neuronal phenotypes relies on extensive temporospatial coordination of mRNA expression to support the development of functional brain circuitry. Cleavage and polyadenylation of mRNA has tremendous regulatory capacity through the alteration of mRNA stability and modulation of microRNA (miRNA) function, however the extent of utilization in neuronal development is currently unclear. Here, we employed poly(A) tail sequencing, mRNA sequencing, ribosome profiling and small RNA sequencing to explore the functional relationship between mRNA abundance, translation, poly(A) tail length, alternative polyadenylation (APA) and miRNA expression in an in vitro model of neuronal differentiation. Differential analysis revealed a strong bias towards poly(A) tail and 3'UTR lengthening during differentiation, both of which were positively correlated with changes in mRNA abundance, but not translation. Globally, changes in miRNA expression were predominantly associated with mRNA abundance and translation, however several miRNA-mRNA pairings with potential to regulate poly(A) tail length were identified. Furthermore, 3'UTR lengthening was observed to significantly increase the inclusion of non-conserved miRNA binding sites, potentially enhancing the regulatory capacity of these molecules in mature neuronal cells. Together, our findings suggest poly(A) tail length and APA function as part of a rich post-transcriptional regulatory matrix during neuronal differentiation.

Essential roles of RNA cap-proximal ribose methylation in mammalian embryonic development and fertility

Eukaryotic RNA pol II transcripts are capped at the 5' end by the methylated guanosine (m⁷G) moiety. In higher eukaryotes, CMTR1 and CMTR2 catalyze cap-proximal ribose methylations on the first (cap1) and second (cap2) nucleotides, respectively. These modifications mark RNAs as "self," blocking the activation of the innate immune response pathway. Here, we show that loss of mouse Cmtr1 or Cmtr2 leads to embryonic lethality, with non-overlapping sets of transcripts being misregulated, but without activation of the interferon pathway. In contrast, Cmtr1 mutant adult mouse livers exhibit chronic activation of the interferon pathway, with multiple interferon-stimulated genes being expressed. Conditional deletion of Cmtr1 in the germline leads to infertility, while global translation is unaffected in the Cmtr1 mutant mouse liver and human cells. Thus, mammalian cap1 and cap2 modifications have essential roles in gene regulation beyond their role in helping cellular transcripts to evade the innate immune system.

Ribosome profiling reveals the translational landscape and allele-specific translational efficiency in rice

Translational regulation is a critical step in the process of gene expression and governs the synthesis of proteins from mRNAs. Many studies have revealed translational regulation in plants in response to various environmental stimuli. However, there have been no studies documenting the comprehensive landscape of translational regulation and allele-specific translational efficiency in multiple plant tissues, especially those of rice, a main staple crop that feeds nearly half of the world's population. Here we used RNA sequencing and ribosome profiling data to analyze the transcriptome and translatome of an elite hybrid rice, Shanyou 63 (SY63), and its parental varieties Zhenshan 97 and Minghui 63. The results revealed that gene expression patterns varied more among tissues than among varieties at the transcriptional and translational levels. We identified 3392 upstream open reading frames (uORFs), and the uORF-containing genes were enriched in transcription factors. Only 668 of 13 492 long non-coding RNAs could be translated into peptides. Finally, we discovered numerous genes with allele-specific translational efficiency in SY63 and demonstrated that some cis-regulatory elements may contribute to allelic divergence in translational efficiency. Overall, these findings may improve our understanding of translational regulation in rice and provide information for molecular breeding research.

Human mitochondria require mtRF1 for translation termination at non-canonical stop codons

The mitochondrial translation machinery highly diverged from its bacterial counterpart. This includes deviation from the universal genetic code, with AGA and AGG codons lacking cognate tRNAs in human mitochondria. The locations of these codons at the end of COX1 and ND6 open reading frames, respectively, suggest they might function as stop codons. However, while the canonical stop codons UAA and UAG are known to be recognized by mtRF1a, the release mechanism at AGA and AGG codons remains a debated issue. Here, we show that upon the loss of another member of the mitochondrial release factor family, mtRF1, mitoribosomes accumulate specifically at AGA and AGG codons. Stalling of mitoribosomes alters COX1 transcript and protein levels, but not ND6 synthesis. In addition, using an in vitro reconstituted mitochondrial translation system, we demonstrate the specific peptide release activity of mtRF1 at the AGA and AGG codons. Together, our results reveal the role of mtRF1 in translation termination at non-canonical stop codons in mitochondria.

Monitoring the 5'UTR landscape reveals isoform switches to drive translational efficiencies in cancer

Transcriptional and translational control are key determinants of gene expression, however, to what extent these two processes can be collectively coordinated is still poorly understood. Here, we use Nanopore long-read sequencing and cap analysis of gene expression (CAGE-seq) to document the landscape of 5' and 3' untranslated region (UTR) isoforms and transcription start sites of epidermal stem cells, wild-type keratinocytes and squamous cell carcinomas. Focusing on squamous cell carcinomas, we show that a small cohort of genes with alternative 5'UTR isoforms exhibit overall increased translational efficiencies and are enriched in ribosomal proteins and splicing factors. By combining polysome fractionations and CAGE-seq, we further characterize two of these UTR isoform genes with identical coding sequences and demonstrate that the underlying transcription start site heterogeneity frequently results in 5' terminal oligopyrimidine (TOP) and pyrimidine-rich translational element (PRTE) motif switches to drive mTORC1-dependent translation of the mRNA. Genome-wide, we show that highly translated squamous cell carcinoma transcripts switch towards increased use of 5'TOP and PRTE motifs, have generally shorter 5'UTRs and expose decreased RNA secondary structures. Notably, we found that the two 5'TOP motif-containing, but not the TOP-less, RPL21 transcript isoforms strongly correlated with overall survival in human head and neck squamous cell carcinoma patients. Our findings warrant isoform-specific analyses in human cancer datasets and suggest that switching between 5'UTR isoforms is an elegant and simple way to alter protein synthesis rates, set their sensitivity to the mTORC1-dependent nutrient-sensing pathway and direct the translational potential of an mRNA by the precise 5'UTR sequence.

A circadian clock translational control mechanism targets specific mRNAs to cytoplasmic messenger ribonucleoprotein granules

Phosphorylation of Neurospora crassa eukaryotic initiation factor 2 α (eIF2α), a conserved translation initiation factor, is clock controlled. To determine the impact of rhythmic eIF2α phosphorylation on translation, we performed temporal ribosome profiling and RNA sequencing (RNA-seq) in wild-type (WT), clock mutant Δfrq, eIF2α kinase mutant Δcpc-3, and constitutively active cpc-3c cells. About 14% of mRNAs are rhythmically translated in WT cells, and translation rhythms for ∼30% of these mRNAs, which we named circadian translation-initiation-controlled genes (cTICs), are dependent on the clock and CPC-3. Most cTICs are expressed from arrhythmic mRNAs and contain a P-body (PB) localization motif in their 5' leader sequence. Deletion of SNR-1, a component of cytoplasmic messenger ribonucleoprotein granules (cmRNPgs) that include PBs and stress granules (SGs), and the PB motif on one of the cTIC mRNAs, zip-1, significantly alters zip-1 rhythmic translation. These results reveal that the clock regulates rhythmic translation of specific mRNAs through rhythmic eIF2α activity and cmRNPg metabolism.

Global and gene-specific translational regulation in Escherichia coli across different conditions

How well mRNA transcript levels represent protein abundances has been a controversial issue. Particularly across different environments, correlations between mRNA and protein exhibit remarkable variability from gene to gene. Translational regulation is likely to be one of the key factors contributing to mismatches between mRNA level and protein abundance in bacteria. Here, we quantified genome-wide transcriptome and relative translation efficiency (RTE) under 12 different conditions in Escherichia coli. By quantifying the mRNA-RTE correlation both across genes and across conditions, we uncovered a diversity of gene-specific translational regulations, cooperating with transcriptional regulations, in response to carbon (C), nitrogen (N), and phosphate (P) limitations. Intriguingly, we found that many genes regulating translation are themselves subject to translational regulation, suggesting possible feedbacks. Furthermore, a random forest model suggests that codon usage partially predicts a gene's cross-condition variability in translation efficiency; such cross-condition variability tends to be an inherent quality of a gene, independent of the specific nutrient limitations. These findings broaden the understanding of translational regulation under different environments and provide novel strategies for the control of translation in synthetic biology. In addition, our data offers a resource for future multi-omics studies.

Tracing Translational Footprint by Ribo-Seq: Principle, Workflow, and Applications to Understand the Mechanism of Human Diseases

RNA-seq has been widely used as a high-throughput method to characterize transcript dynamic changes in a broad context, such as development and diseases. However, whether RNA-seq-estimated transcriptional dynamics can be translated into protein level changes is largely unknown. Ribo-seq (Ribosome profiling) is an emerging technology that allows for the investigation of the translational footprint via profiling ribosome-bounded mRNA fragments. Ribo-seq coupled with RNA-seq will allow us to understand the transcriptional and translational control of the fundamental biological process and human diseases. This review focuses on discussing the principle, workflow, and applications of Ribo-seq to study human diseases.

End-to-end pipeline for differential analysis of pausing in ribosome profiling data

Ribosome profiling is a powerful technique which maps the distribution of ribosomes along mRNAs to analyze translation genome-wide. Ribosome density can be affected by multiple factors, such as changes to translation initiation or elongation rates. We describe the application of a metric for identifying genes rate-limited by these rates by analyzing the relative distribution of ribosome footprints along transcripts. This protocol also details two sample analyses comparing gene translation efficiencies and the distribution of ribosome densities on downloadable datasets.

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

•

Pipeline for measuring translation efficiency using ribosome profiling data

•

Quantification of differential ribosome distributions along mRNAs

•

Determine whether changes are consistent with altered initiation or elongation

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

Assessing and assuring interoperability of a genomics file format

Motivation

Bioinformatics software tools operate largely through the use of specialized genomics file formats. Often these formats lack formal specification, making it difficult or impossible for the creators of these tools to robustly test them for correct handling of input and output. This causes problems in interoperability between different tools that, at best, wastes time and frustrates users. At worst, interoperability issues could lead to undetected errors in scientific results.

Results

We developed a new verification system, Acidbio, which tests for correct behavior in bioinformatics software packages. We crafted tests to unify correct behavior when tools encounter various edge cases—potentially unexpected inputs that exemplify the limits of the format. To analyze the performance of existing software, we tested the input validation of 80 Bioconda packages that parsed the Browser Extensible Data (BED) format. We also used a fuzzing approach to automatically perform additional testing. Of 80 software packages examined, 75 achieved less than 70% correctness on our test suite. We categorized multiple root causes for the poor performance of different types of software. Fuzzing detected other errors that the manually designed test suite could not. We also created a badge system that developers can use to indicate more precisely which BED variants their software accepts and to advertise the software’s performance on the test suite.

Availability and implementation

Acidbio is available at https://github.com/hoffmangroup/acidbio.

Supplementary information

Supplementary data are available at Bioinformatics online.

Critical contribution of 3' non-seed base pairing to the in vivo function of the evolutionarily conserved let-7a microRNA

Base pairing of the seed region (g2–g8) is essential for microRNA targeting; however, the in vivo function of the 3′ non-seed region (g9–g22) is less well understood. Here, we report a systematic investigation of the in vivo roles of 3′ non-seed nucleotides in microRNA let-7a, whose entire g9–g22 region is conserved among bilaterians. We find that the 3′ non-seed sequence functionally distinguishes let-7a from its family paralogs. The complete pairing of g11–g16 is essential for let-7a to fully repress multiple key targets, including evolutionarily conserved lin-41, daf-12, and hbl-1. Nucleotides at g17–g22 are less critical but may compensate for mismatches in the g11–g16 region. Interestingly, a certain minimal complementarity to let-7a 3′ non-seed sequence can be required even for sites with perfect seed pairing. These results provide evidence that the specific configurations of both seed and 3′ non-seed base pairing can critically influence microRNA-mediated gene regulation in vivo.

Duan et al. find that microRNA-target pairing at g11–g16 is critical for the function of evolutionarily conserved microRNA let-7a; 3′ pairing is required for both perfect and imperfect seed in regulating multiple targets. These findings provide evidence that base pairing of specific microRNA non-seed nucleotides can critically contribute to target regulation.

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

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

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

Functional single-cell genomics of human cytomegalovirus infection

Understanding how viral and host factors interact and how perturbations impact infection is the basis for designing antiviral interventions. Here we define the functional contribution of each viral and host factor involved in human cytomegalovirus infection in primary human fibroblasts through pooled CRISPR interference and nuclease screening. To determine how genetic perturbation of critical host and viral factors alters the timing, course and progression of infection, we applied Perturb-seq to record the transcriptomes of tens of thousands of CRISPR-modified single cells and found that, normally, most cells follow a stereotypical transcriptional trajectory. Perturbing critical host factors does not change the stereotypical transcriptional trajectory per se but can stall, delay or accelerate progression along the trajectory, allowing one to pinpoint the stage of infection at which host factors act. Conversely, perturbation of viral factors can create distinct, abortive trajectories. Our results reveal the roles of host and viral factors and provide a roadmap for the dissection of host-pathogen interactions.

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.

Small Open Reading Frames, How to Find Them and Determine Their Function

Advances in genomics and molecular biology have revealed an abundance of small open reading frames (sORFs) across all types of transcripts. While these sORFs are often assumed to be non-functional, many have been implicated in physiological functions and a significant number of sORFs have been described in human diseases. Thus, sORFs may represent a hidden repository of functional elements that could serve as therapeutic targets. Unlike protein-coding genes, it is not necessarily the encoded peptide of an sORF that enacts its function, sometimes simply the act of translating an sORF might have a regulatory role. Indeed, the most studied sORFs are located in the 5′UTRs of coding transcripts and can have a regulatory impact on the translation of the downstream protein-coding sequence. However, sORFs have also been abundantly identified in non-coding RNAs including lncRNAs, circular RNAs and ribosomal RNAs suggesting that sORFs may be diverse in function. Of the many different experimental methods used to discover sORFs, the most commonly used are ribosome profiling and mass spectrometry. These can confirm interactions between transcripts and ribosomes and the production of a peptide, respectively. Extensions to ribosome profiling, which also capture scanning ribosomes, have further made it possible to see how sORFs impact the translation initiation of mRNAs. While high-throughput techniques have made the identification of sORFs less difficult, defining their function, if any, is typically more challenging. Together, the abundance and potential function of many of these sORFs argues for the necessity of including sORFs in gene annotations and systematically characterizing these to understand their potential functional roles. In this review, we will focus on the high-throughput methods used in the detection and characterization of sORFs and discuss techniques for validation and functional characterization.

Ribo-Seq and RNA-Seq of TMA46 ( DFRP1) and GIR2 ( DFRP2) knockout yeast strains

In eukaryotes, stalled and collided ribosomes are recognized by several conserved multicomponent systems, which either block protein synthesis in situ and resolve the collision locally, or trigger a general stress response. Yeast ribosome-binding GTPases RBG1 (DRG1 in mammals) and RBG2 (DRG2) form two distinct heterodimers with TMA46 (DFRP1) and GIR2 (DFRP2), respectively, both involved in mRNA translation. Accumulated evidence suggests that the dimers play partially redundant roles in elongation processivity and resolution of ribosome stalling and collision events, as well as in the regulation of GCN1-mediated signaling involved in ribosome-associated quality control (RQC). They also genetically interact with SLH1 (ASCC3) helicase, a key component of RQC trigger (RQT) complex disassembling collided ribosomes. Here, we present RNA-Seq and ribosome profiling (Ribo-Seq) data from S. cerevisiae strains with individual deletions of the TMA46 and GIR2 genes. Raw RNA-Seq and Ribo-Seq data as well as gene-level read counts are available in NCBI Gene Expression Omnibus (GEO) repository under GEO accession GSE185458 and GSE185286.

Identification of Non-Canonical Translation Products in C. elegans Using Tandem Mass Spectrometry

Transcriptome and ribosome sequencing have revealed the existence of many non-canonical transcripts, mainly containing splice variants, ncRNA, sORFs and altORFs. However, identification and characterization of products that may be translated out of these remains a challenge. Addressing this, we here report on 552 non-canonical proteins and splice variants in the model organism C. elegans using tandem mass spectrometry. Aided by sequencing-based prediction, we generated a custom proteome database tailored to search for non-canonical translation products of C. elegans. Using this database, we mined available mass spectrometric resources of C. elegans, from which 51 novel, non-canonical proteins could be identified. Furthermore, we utilized diverse proteomic and peptidomic strategies to detect 40 novel non-canonical proteins in C. elegans by LC-TIMS-MS/MS, of which 6 were common with our meta-analysis of existing resources. Together, this permits us to provide a resource with detailed annotation of 467 splice variants and 85 novel proteins mapped onto UTRs, non-coding regions and alternative open reading frames of the C. elegans genome.

Testing the Genomic Shock Hypothesis Using Transposable Element Expression in Yeast Hybrids

Transposable element (TE) insertions are a source of structural variation and can cause genetic instability and gene expression changes. A host can limit the spread of TEs with various repression mechanisms. Many examples of plant and animal interspecific hybrids show disrupted TE repression leading to TE propagation. Recent studies in yeast did not find any increase in transposition rate in hybrids. However, this does not rule out the possibility that the transcriptional or translational activity of TEs increases following hybridization because of a disruption of the host TE control mechanisms. Thus, whether total expression of a TE family is higher in hybrids than in their parental species remains to be examined. We leveraged publically available RNA-seq and ribosomal profiling data on yeast artificial hybrids of the Saccharomyces genus and performed differential expression analysis of their LTR retrotransposons (Ty elements). Our analyses of total mRNA levels show that Ty elements are generally not differentially expressed in hybrids, even when the hybrids are exposed to a low temperature stress condition. Overall, only 2/26 Ty families show significantly higher expression in the S. cerevisiae × S. uvarum hybrids while there are 3/26 showing significantly lower expression in the S. cerevisiae x S. paradoxus hybrids. Our analysis of ribosome profiling data of S. cerevisiae × S. paradoxus hybrids shows similar translation efficiency of Ty in both parents and hybrids, except for Ty1_cer showing higher translation efficiency. Overall, our results do not support the hypothesis that hybridization could act as a systematic trigger of TE expression in yeast and suggest that the impact of hybridization on TE activity is strain and TE specific.

Activation of the antiviral factor RNase L triggers translation of non-coding mRNA sequences

Ribonuclease L (RNase L) is activated as part of the innate immune response and plays an important role in the clearance of viral infections. When activated, it endonucleolytically cleaves both viral and host RNAs, leading to a global reduction in protein synthesis. However, it remains unknown how widespread RNA decay, and consequent changes in the translatome, promote the elimination of viruses. To study how this altered transcriptome is translated, we assayed the global distribution of ribosomes in RNase L activated human cells with ribosome profiling. We found that RNase L activation leads to a substantial increase in the fraction of translating ribosomes in ORFs internal to coding sequences (iORFs) and ORFs within 5' and 3' UTRs (uORFs and dORFs). Translation of these alternative ORFs was dependent on RNase L's cleavage activity, suggesting that mRNA decay fragments are translated to produce short peptides that may be important for antiviral activity.

Trips-Viz: an environment for the analysis of public and user-generated ribosome profiling data

Trips-Viz (https://trips.ucc.ie/) is an interactive platform for the analysis and visualization of ribosome profiling (Ribo-Seq) and shotgun RNA sequencing (RNA-seq) data. This includes publicly available and user generated data, hence Trips-Viz can be classified as a database and as a server. As a database it provides access to many processed Ribo-Seq and RNA-seq data aligned to reference transcriptomes which has been expanded considerably since its inception. Here, we focus on the server functionality of Trips-viz which also has been greatly improved. Trips-viz now enables visualisation of proteomics data from a large number of processed mass spectrometry datasets. It can be used to support translation inferred from Ribo-Seq data. Users are now able to upload a custom reference transcriptome as well as data types other than Ribo-Seq/RNA-Seq. Incorporating custom data has been streamlined with RiboGalaxy (https://ribogalaxy.ucc.ie/) integration. The other new functionality is the rapid detection of translated open reading frames (ORFs) through a simple easy to use interface. The analysis of differential expression has been also improved via integration of DESeq2 and Anota2seq in addition to a number of other improvements of existing Trips-viz features.

miCLIP-MaPseq Identifies Substrates of Radical SAM RNA-Methylating Enzyme Using Mechanistic Cross-Linking and Mismatch Profiling

The family of radical SAM RNA-methylating enzymes comprises a large group of proteins that contains only a few functionally characterized members. Several enzymes in this family have been implicated in the regulation of translation and antibiotic susceptibility, emphasizing their significance in bacterial physiology and their relevance to human health. While few characterized enzymes have been shown to modify diverse RNA substrates, highlighting potentially broad substrate scope within the family, many enzymes in this class have no known substrates. The precise knowledge of RNA substrates and modification sites for uncharacterized family members is important for unraveling their biological function. Here, we describe a strategy for substrate identification that takes advantage of mechanism-based cross-linking between the enzyme and its RNA substrates, which we named individual-nucleotide-resolution cross-linking and immunoprecipitation combined with mutational profiling with sequencing (miCLIP-MaPseq). Identification of the position of the modification site is achieved using thermostable group II intron reverse transcriptase (TGIRT), which introduces a mismatch at the site of the cross-link.

ORFik: a comprehensive R toolkit for the analysis of translation

BACKGROUND

With the rapid growth in the use of high-throughput methods for characterizing translation and the continued expansion of multi-omics, there is a need for back-end functions and streamlined tools for processing, analyzing, and characterizing data produced by these assays.

RESULTS

Here, we introduce ORFik, a user-friendly R/Bioconductor API and toolbox for studying translation and its regulation. It extends GenomicRanges from the genome to the transcriptome and implements a framework that integrates data from several sources. ORFik streamlines the steps to process, analyze, and visualize the different steps of translation with a particular focus on initiation and elongation. It accepts high-throughput sequencing data from ribosome profiling to quantify ribosome elongation or RCP-seq/TCP-seq to also quantify ribosome scanning. In addition, ORFik can use CAGE data to accurately determine 5'UTRs and RNA-seq for determining translation relative to RNA abundance. ORFik supports and calculates over 30 different translation-related features and metrics from the literature and can annotate translated regions such as proteins or upstream open reading frames (uORFs). As a use-case, we demonstrate using ORFik to rapidly annotate the dynamics of 5' UTRs across different tissues, detect their uORFs, and characterize their scanning and translation in the downstream protein-coding regions.

CONCLUSION

In summary, ORFik introduces hundreds of tested, documented and optimized methods. ORFik is designed to be easily customizable, enabling users to create complete workflows from raw data to publication-ready figures for several types of sequencing data. Finally, by improving speed and scope of many core Bioconductor functions, ORFik offers enhancement benefiting the entire Bioconductor environment.

AVAILABILITY

http://bioconductor.org/packages/ORFik .

smORFer: a modular algorithm to detect small ORFs in prokaryotes

Emerging evidence places small proteins (≤50 amino acids) more centrally in physiological processes. Yet, their functional identification and the systematic genome annotation of their cognate small open-reading frames (smORFs) remains challenging both experimentally and computationally. Ribosome profiling or Ribo-Seq (that is a deep sequencing of ribosome-protected fragments) enables detecting of actively translated open-reading frames (ORFs) and empirical annotation of coding sequences (CDSs) using the in-register translation pattern that is characteristic for genuinely translating ribosomes. Multiple identifiers of ORFs that use the 3-nt periodicity in Ribo-Seq data sets have been successful in eukaryotic smORF annotation. They have difficulties evaluating prokaryotic genomes due to the unique architecture (e.g. polycistronic messages, overlapping ORFs, leaderless translation, non-canonical initiation etc.). Here, we present a new algorithm, smORFer, which performs with high accuracy in prokaryotic organisms in detecting putative smORFs. The unique feature of smORFer is that it uses an integrated approach and considers structural features of the genetic sequence along with in-frame translation and uses Fourier transform to convert these parameters into a measurable score to faithfully select smORFs. The algorithm is executed in a modular way, and dependent on the data available for a particular organism, different modules can be selected for smORF search.

Persistence of Ambigrammatic Narnaviruses Requires Translation of the Reverse Open Reading Frame

Narnaviruses are RNA viruses detected in diverse fungi, plants, protists, arthropods, and nematodes. Though initially described as simple single-gene nonsegmented viruses encoding RNA-dependent RNA polymerase (RdRp), a subset of narnaviruses referred to as "ambigrammatic" harbor a unique genomic configuration consisting of overlapping open reading frames (ORFs) encoded on opposite strands. Phylogenetic analysis supports selection to maintain this unusual genome organization, but functional investigations are lacking. Here, we establish the mosquito-infecting Culex narnavirus 1 (CxNV1) as a model to investigate the functional role of overlapping ORFs in narnavirus replication. In CxNV1, a reverse ORF without homology to known proteins covers nearly the entire 3.2-kb segment encoding the RdRp. Additionally, two opposing and nearly completely overlapping novel ORFs are found on the second putative CxNV1 segment, the 0.8-kb "Robin" RNA. We developed a system to launch CxNV1 in a naive mosquito cell line and then showed that functional RdRp is required for persistence of both segments, and an intact reverse ORF is required on the RdRp segment for persistence. Mass spectrometry of persistently CxNV1-infected cells provided evidence for translation of this reverse ORF. Finally, ribosome profiling yielded a striking pattern of footprints for all four CxNV1 RNA strands that was distinct from actively translating ribosomes on host mRNA or coinfecting RNA viruses. Taken together, these data raise the possibility that the process of translation itself is important for persistence of ambigrammatic narnaviruses, potentially by protecting viral RNA with ribosomes, thus suggesting a heretofore undescribed viral tactic for replication and transmission. IMPORTANCE Fundamental to our understanding of RNA viruses is a description of which strand(s) of RNA are transmitted as the viral genome relative to which encode the viral proteins. Ambigrammatic narnaviruses break the mold. These viruses, found broadly in fungi, plants, and insects, have the unique feature of two overlapping genes encoded on opposite strands, comprising nearly the full length of the viral genome. Such extensive overlap is not seen in other RNA viruses and comes at the cost of reduced evolutionary flexibility in the sequence. The present study is motivated by investigating the benefits which balance that cost. We show for the first time a functional requirement for the ambigrammatic genome configuration in Culex narnavirus 1, which suggests a model for how translation of both strands might benefit this virus. Our work highlights a new blueprint for viral persistence, distinct from strategies defined by canonical definitions of the coding strand.

Genetic removal of p70 S6K1 corrects coding sequence length-dependent alterations in mRNA translation in fragile X syndrome mice

Loss of the fragile X mental retardation protein (FMRP) causes fragile X syndrome (FXS). FMRP is widely thought to repress protein synthesis, but its translational targets and modes of control remain in dispute. We previously showed that genetic removal of p70 S6 kinase 1 (S6K1) corrects altered protein synthesis as well as synaptic and behavioral phenotypes in FXS mice. In this study, we examined the gene specificity of altered messenger RNA (mRNA) translation in FXS and the mechanism of rescue with genetic reduction of S6K1 by carrying out ribosome profiling and RNA sequencing on cortical lysates from wild-type, FXS, S6K1 knockout, and double knockout mice. We observed reduced ribosome footprint (RF) abundance in the majority of differentially translated genes in the cortices of FXS mice. We used molecular assays to discover evidence that the reduction in RF abundance reflects an increased rate of ribosome translocation, which is captured as a decrease in the number of translating ribosomes at steady state and is normalized by inhibition of S6K1. We also found that genetic removal of S6K1 prevented a positive-to-negative gradation of alterations in translation efficiencies (RF/mRNA) with coding sequence length across mRNAs in FXS mouse cortices. Our findings reveal the identities of dysregulated mRNAs and a molecular mechanism by which reduction of S6K1 prevents altered translation in FXS.

CRISPR-based functional genomics in human dendritic cells

Dendritic cells (DCs) regulate processes ranging from antitumor and antiviral immunity to host-microbe communication at mucosal surfaces. It remains difficult, however, to genetically manipulate human DCs, limiting our ability to probe how DCs elicit specific immune responses. Here, we develop a CRISPR-Cas9 genome editing method for human monocyte-derived DCs (moDCs) that mediates knockouts with a median efficiency of >94% across >300 genes. Using this method, we perform genetic screens in moDCs, identifying mechanisms by which DCs tune responses to lipopolysaccharides from the human microbiome. In addition, we reveal donor-specific responses to lipopolysaccharides, underscoring the importance of assessing immune phenotypes in donor-derived cells, and identify candidate genes that control this specificity, highlighting the potential of our method to pinpoint determinants of inter-individual variation in immunity. Our work sets the stage for a systematic dissection of the immune signaling at the host-microbiome interface and for targeted engineering of DCs for neoantigen vaccination.

Synonymous mutations that regulate translation speed might play a non-negligible role in liver cancer development

Background

Synonymous mutations do not change the protein sequences. Automatically, they have been regarded as neutral events and are ignored in the mutation-based cancer studies. However, synonymous mutations will change the codon optimality, resulting in altered translational velocity.

Methods

We fully utilized the transcriptome and translatome of liver cancer and normal tissue from ten patients. We profiled the mutation spectrum and examined the effect of synonymous mutations on translational velocity.

Results

Synonymous mutations that increase the codon optimality significantly enhanced the translational velocity, and were enriched in oncogenes. Meanwhile, synonymous mutations decreasing codon optimality slowed down translation, and were enriched in tumor suppressor genes. These synonymous mutations significantly contributed to the translational changes in tumor samples compared to normal samples.

Conclusions

Synonymous mutations might play a role in liver cancer development by altering codon optimality and translational velocity. Synonymous mutations should no longer be ignored in the genome-wide studies.

RiboA: a web application to identify ribosome A-site locations in ribosome profiling data

Background

Translation is a fundamental process in gene expression. Ribosome profiling is a method that enables the study of transcriptome-wide translation. A fundamental, technical challenge in analyzing Ribo-Seq data is identifying the A-site location on ribosome-protected mRNA fragments. Identification of the A-site is essential as it is at this location on the ribosome where a codon is translated into an amino acid. Incorrect assignment of a read to the A-site can lead to lower signal-to-noise ratio and loss of correlations necessary to understand the molecular factors influencing translation. Therefore, an easy-to-use and accurate analysis tool is needed to accurately identify the A-site locations.

Results

We present RiboA, a web application that identifies the most accurate A-site location on a ribosome-protected mRNA fragment and generates the A-site read density profiles. It uses an Integer Programming method that reflects the biological fact that the A-site of actively translating ribosomes is generally located between the second codon and stop codon of a transcript, and utilizes a wide range of mRNA fragment sizes in and around the coding sequence (CDS). The web application is containerized with Docker, and it can be easily ported across platforms.

Conclusions

The Integer Programming method that RiboA utilizes is the most accurate in identifying the A-site on Ribo-Seq mRNA fragments compared to other methods. RiboA makes it easier for the community to use this method via a user-friendly and portable web application. In addition, RiboA supports reproducible analyses by tracking all the input datasets and parameters, and it provides enhanced visualization to facilitate scientific exploration. RiboA is available as a web service at https://a-site.vmhost.psu.edu/. The code is publicly available at https://github.com/obrien-lab/aip_web_docker under the MIT license.

Regulation of translation by methylation multiplicity of 18S rRNA

N6-methyladenosine (m6A) is a conserved ribonucleoside modification that regulates many facets of RNA metabolism. Using quantitative mass spectrometry, we find that the universally conserved tandem adenosines at the 3' end of 18S rRNA, thought to be constitutively di-methylated (m62A), are also mono-methylated (m6A). Although present at substoichiometric amounts, m6A at these positions increases significantly in response to sulfur starvation in yeast cells and mammalian cell lines. Combining yeast genetics and ribosome profiling, we provide evidence to suggest that m6A-bearing ribosomes carry out translation distinctly from m62A-bearing ribosomes, featuring a striking specificity for sulfur metabolism genes. Our work thus reveals methylation multiplicity as a mechanism to regulate translation.

The ribosome collision sensor Hel2 functions as preventive quality control in the secretory pathway

Ribosome collision because of translational stalling is recognized as a problematic event in translation by the E3 ubiquitin ligase Hel2, leading to non-canonical subunit dissociation followed by targeting of the faulty nascent peptides for degradation. Although Hel2-mediated quality control greatly contributes to maintenance of cellular protein homeostasis, its physiological role in dealing with endogenous substrates remains unclear. This study utilizes genome-wide analysis, based on selective ribosome profiling, to survey the endogenous substrates for Hel2. This survey reveals that Hel2 binds preferentially to the pre-engaged secretory ribosome-nascent chain complexes (RNCs), which translate upstream of targeting signals. Notably, Hel2 recruitment into secretory RNCs is elevated under signal recognition particle (SRP)-deficient conditions. Moreover, the mitochondrial defects caused by insufficient SRP are enhanced by hel2 deletion, along with mistargeting of secretory proteins into mitochondria. These findings provide insights into risk management in the secretory pathway that maintains cellular protein homeostasis.

Determinants of genome-wide distribution and evolution of uORFs in eukaryotes

Upstream open reading frames (uORFs) play widespread regulatory functions in modulating mRNA translation in eukaryotes, but the principles underlying the genomic distribution and evolution of uORFs remain poorly understood. Here, we analyze ~17 million putative canonical uORFs in 478 eukaryotic species that span most of the extant taxa of eukaryotes. We demonstrate how positive and purifying selection, coupled with differences in effective population size (Ne), has shaped the contents of uORFs in eukaryotes. Besides, gene expression level is important in influencing uORF occurrences across genes in a species. Our analyses suggest that most uORFs might play regulatory roles rather than encode functional peptides. We also show that the Kozak sequence context of uORFs has evolved across eukaryotic clades, and that noncanonical uORFs tend to have weaker suppressive effects than canonical uORFs in translation regulation. This study provides insights into the driving forces underlying uORF evolution in eukaryotes.

A stress-induced tyrosine-tRNA depletion response mediates codon-based translational repression and growth suppression

Eukaryotic transfer RNAs can become selectively fragmented upon various stresses, generating tRNA-derived small RNA fragments. Such fragmentation has been reported to impact a small fraction of the tRNA pool and thus presumed to not directly impact translation. We report that oxidative stress can rapidly generate tyrosine-tRNAGUA fragments in human cells-causing significant depletion of the precursor tRNA. Tyrosine-tRNAGUA depletion impaired translation of growth and metabolic genes enriched in cognate tyrosine codons. Depletion of tyrosine tRNAGUA or its translationally regulated targets USP3 and SCD repressed proliferation-revealing a dedicated tRNA-regulated growth-suppressive pathway for oxidative stress response. Tyrosine fragments are generated in a DIS3L2 exoribonuclease-dependent manner and inhibit hnRNPA1-mediated transcript destabilization. Moreover, tyrosine fragmentation is conserved in C. elegans. Thus, tRNA fragmentation can coordinately generate trans-acting small RNAs and functionally deplete a tRNA. Our findings reveal the existence of an underlying adaptive codon-based regulatory response inherent to the genetic code.

LncRNA-Encoded Peptide: Functions and Predicting Methods

Long non-coding RNA (lncRNA) was originally defined as the representative of the non-coding RNAs and unable to encode. However, recent reports suggest that some lncRNAs actually contain open reading frames that encode peptides. These coding products play important roles in the pathogenesis of many diseases. Here, we summarize the regulatory pathways of mammalian lncRNA-encoded peptides in influencing muscle function, mRNA stability, gene expression, and so on. We also address the promoting and inhibiting functions of the peptides in different cancers and other diseases. Then we introduce the computational predicting methods and data resources to predict the coding ability of lncRNA. The intention of this review is to provide references for further coding research and contribute to reveal the potential prospects for targeted tumor therapy.

Codon Resolution Analysis of Ribosome Profiling Data

Translation is a central biological process in living cells. Ribosome profiling approach enables assessing translation on a global, cell-wide level. Extracting versatile information from the ribosome profiling data usually requires specialized expertise for handling the sequencing data that is not available to the broad community of experimentalists. Here, we provide an easy-to-use and modifiable workflow that uses a small set of commands and enables full data analysis in a standardized way, including precise positioning of the ribosome-protected fragments, for determining codon-specific translation features. The workflow is complemented with simple step-by-step explanations and is accessible to scientists with no computational background.

Protein Synthesis in the Developing Neocortex at Near-Atomic Resolution Reveals Ebp1-Mediated Neuronal Proteostasis at the 60S Tunnel Exit

Protein synthesis must be finely tuned in the developing nervous system as the final essential step of gene expression. This study investigates the architecture of ribosomes from the neocortex during neurogenesis, revealing Ebp1 as a high-occupancy 60S peptide tunnel exit (TE) factor during protein synthesis at near-atomic resolution by cryoelectron microscopy (cryo-EM). Ribosome profiling demonstrated Ebp1-60S binding is highest during start codon initiation and N-terminal peptide elongation, regulating ribosome occupancy of these codons. Membrane-targeting domains emerging from the 60S tunnel, which recruit SRP/Sec61 to the shared binding site, displace Ebp1. Ebp1 is particularly abundant in the early-born neural stem cell (NSC) lineage and regulates neuronal morphology. Ebp1 especially impacts the synthesis of membrane-targeted cell adhesion molecules (CAMs), measured by pulsed stable isotope labeling by amino acids in cell culture (pSILAC)/bioorthogonal noncanonical amino acid tagging (BONCAT) mass spectrometry (MS). Therefore, Ebp1 is a central component of protein synthesis, and the ribosome TE is a focal point of gene expression control in the molecular specification of neuronal morphology during development.

RiboDiPA: a novel tool for differential pattern analysis in Ribo-seq data

Ribosome profiling, also known as Ribo-seq, has become a popular approach to investigate regulatory mechanisms of translation in a wide variety of biological contexts. Ribo-seq not only provides a measurement of translation efficiency based on the relative abundance of ribosomes bound to transcripts, but also has the capacity to reveal dynamic and local regulation at different stages of translation based on positional information of footprints across individual transcripts. While many computational tools exist for the analysis of Ribo-seq data, no method is currently available for rigorous testing of the pattern differences in ribosome footprints. In this work, we develop a novel approach together with an R package, RiboDiPA, for Differential Pattern Analysis of Ribo-seq data. RiboDiPA allows for quick identification of genes with statistically significant differences in ribosome occupancy patterns for model organisms ranging from yeast to mammals. We show that differential pattern analysis reveals information that is distinct and complimentary to existing methods that focus on translational efficiency analysis. Using both simulated Ribo-seq footprint data and three benchmark data sets, we illustrate that RiboDiPA can uncover meaningful pattern differences across multiple biological conditions on a global scale, and pinpoint characteristic ribosome occupancy patterns at single codon resolution.

What determines eukaryotic translation elongation: recent molecular and quantitative analyses of protein synthesis

Protein synthesis from mRNA is an energy-intensive and tightly controlled cellular process. Translation elongation is a well-coordinated, multifactorial step in translation that undergoes dynamic regulation owing to cellular state and environmental determinants. Recent studies involving genome-wide approaches have uncovered some crucial aspects of translation elongation including the mRNA itself and the nascent polypeptide chain. Additionally, these studies have fuelled quantitative and mathematical modelling of translation elongation. In this review, we provide a comprehensive overview of the key determinants of translation elongation. We discuss consequences of ribosome stalling or collision, and how the cells regulate translation in case of such events. Next, we review theoretical approaches and widely used mathematical models that have become an essential ingredient to interpret complex molecular datasets and study translation dynamics quantitatively. Finally, we review recent advances in live-cell reporter and related analysis techniques, to monitor the translation dynamics of single cells and single-mRNA molecules in real time.

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.

Improved computational analysis of ribosome dynamics from 5'P degradome data using fivepseq

In eukaryotes, 5′–3′ co-translation degradation machinery follows the last translating ribosome providing an in vivo footprint of its position. Thus, 5′ monophosphorylated (5′P) degradome sequencing, in addition to informing about RNA decay, also provides information regarding ribosome dynamics. Multiple experimental methods have been developed to investigate the mRNA degradome; however, computational tools for their reproducible analysis are lacking. Here, we present fivepseq: an easy-to-use application for analysis and interactive visualization of 5′P degradome data. This tool performs both metagene- and gene-specific analysis, and enables easy investigation of codon-specific ribosome pauses. To demonstrate its ability to provide new biological information, we investigate gene-specific ribosome pauses in Saccharomyces cerevisiae after eIF5A depletion. In addition to identifying pauses at expected codon motifs, we identify multiple genes with strain-specific degradation frameshifts. To show its wide applicability, we investigate 5′P degradome from Arabidopsis thaliana and discover both motif-specific ribosome protection associated with particular developmental stages and generally increased ribosome protection at termination level associated with age. Our work shows how the use of improved analysis tools for the study of 5′P degradome can significantly increase the biological information that can be derived from such datasets and facilitate its reproducible analysis.