The use of duplex-specific nuclease in ribosome profiling and a user-friendly software package for Ribo-seq data analysis

Advances in ribosome profiling technologies

Ribosome profiling (or Ribo-seq) has emerged as a powerful approach for revealing the regulatory mechanisms of protein synthesis, on the basis of deep sequencing of ribosome footprints. Recent innovations in Ribo-seq technologies have significantly enhanced their sensitivity, specificity, and resolution. In this review, we outline emerging Ribo-seq derivatives that overcome barriers in low inputs, rRNA contamination, data calibration, and single-cell applications. These advances enable detailed insights into translational control across diverse biological contexts.

Ribo-seq guided design of enhanced protein secretion in Komagataellaphaffii

The production of recombinant proteins requires the precise coordination of various biological processes, including protein synthesis, folding, trafficking, and secretion. The overproduction of a heterologous protein can impose various bottlenecks on these networks. Identifying and alleviating these bottlenecks can guide strain engineering efforts to enhance protein production. The methylotrophic yeast Komagataella phaffii is used for its high capacity to produce recombinant proteins. Here, we use ribosome profiling to identify bottlenecks in protein secretion during heterologous expression of human serum albumin (HSA). Validation of this analysis showed that the knockout of non-essential genes whose gene products target the ER, through co- and post-translational mechanisms, and have high ribosome utilization can increase production of a heterologous protein, HSA. A triple knockout in co-translationally translocated carbohydrate and acetate transporter Gal2p, cell wall maintenance protein Ydr134cp, and the post-translationally translocated cell wall protein Aoa65896.1 increased HSA production by 35 %. This data-driven strain engineering approach uses cell-level information to identify gene targets for phenotype improvement. This specific case identifies hits and creates strains with improved HSA production, with Ribo-seq and bioinformatic analysis to identify non-essential ER targeted proteins that are high ribosome utilizers.

Zika viruses encode 5' upstream open reading frames affecting infection of human brain cells

Zika virus (ZIKV), an emerging mosquito-borne flavivirus, is associated with congenital neurological complications. Here, we investigate potential pathological correlates of virus gene expression in representative ZIKV strains through RNA sequencing and ribosome profiling. In addition to the single long polyprotein found in all flaviviruses, we identify the translation of unrecognised upstream open reading frames (uORFs) in the genomic 5' region. In Asian/American strains, ribosomes translate uORF1 and uORF2, whereas in African strains, the two uORFs are fused into one (African uORF). We use reverse genetics to examine the impact on ZIKV fitness of different uORFs mutant viruses. We find that expression of the African uORF and the Asian/American uORF1 modulates virus growth and tropism in human cortical neurons and cerebral organoids, suggesting a potential role in neurotropism. Although the uORFs are expressed in mosquito cells, we do not see a measurable effect on transmission by the mosquito vector in vivo. The discovery of ZIKV uORFs sheds new light on the infection of the human brain cells by this virus and raises the question of their existence in other neurotropic flaviviruses.

Optimization of ribosome profiling in plants including structural analysis of rRNA fragments

BACKGROUND

Ribosome profiling (or Ribo-seq) is a technique that provides genome-wide information on the translational landscape (translatome). Across different plant studies, variable methodological setups have been described which raises questions about the general comparability of data that were generated from diverging methodologies. Furthermore, a common problem when performing Ribo-seq are abundant rRNA fragments that are wastefully incorporated into the libraries and dramatically reduce sequencing depth. To remove these rRNA contaminants, it is common to perform preliminary trials to identify these fragments because they are thought to vary depending on nuclease treatment, tissue source, and plant species.

RESULTS

Here, we compile valuable insights gathered over years of generating Ribo-seq datasets from different species and experimental setups. We highlight which technical steps are important for maintaining cross experiment comparability and describe a highly efficient approach for rRNA removal. Furthermore, we provide evidence that many rRNA fragments are structurally preserved over diverse nuclease regimes, as well as across plant species. Using a recently published cryo-electron microscopy (cryo-EM) structure of the tobacco 80S ribosome, we show that the most abundant rRNA fragments are spatially derived from the solvent-exposed surface of the ribosome.

CONCLUSION

The guidelines presented here shall aid newcomers in establishing ribosome profiling in new plant species and provide insights that will help in customizing the methodology for individual research goals.

Calibrated ribosome profiling assesses the dynamics of ribosomal flux on transcripts

Ribosome profiling, which is based on deep sequencing of ribosome footprints, has served as a powerful tool for elucidating the regulatory mechanism of protein synthesis. However, the current method has substantial issues: contamination by rRNAs and the lack of appropriate methods to measure ribosome numbers in transcripts. Here, we overcome these hurdles through the development of "Ribo-FilterOut", which is based on the separation of footprints from ribosome subunits by ultrafiltration, and "Ribo-Calibration", which relies on external spike-ins of stoichiometrically defined mRNA-ribosome complexes. A combination of these approaches estimates the number of ribosomes on a transcript, the translation initiation rate, and the overall number of translation events before its decay, all in a genome-wide manner. Moreover, our method reveals the allocation of ribosomes under heat shock stress, during aging, and across cell types. Our strategy of modified ribosome profiling measures kinetic and stoichiometric parameters of cellular translation across the transcriptome.

A novel ilarvirus protein CP-RT is expressed via stop codon readthrough and suppresses RDR6-dependent RNA silencing

Ilarviruses are a relatively understudied but important group of plant RNA viruses that includes a number of crop pathogens. Their genomes comprise three RNA segments encoding two replicase subunits, movement protein, coat protein (CP), and (in some ilarvirus subgroups) a protein that suppresses RNA silencing. Here we report that, in many ilarviruses, RNA3 encodes an additional protein (termed CP-RT) as a result of ribosomal readthrough of the CP stop codon into a short downstream readthrough (RT) ORF. Using asparagus virus 2 as a model, we find that CP-RT is expressed in planta where it functions as a weak suppressor of RNA silencing. CP-RT expression is essential for persistent systemic infection in leaves and shoot apical meristem. CP-RT function is dependent on a putative zinc-finger motif within RT. Replacing the asparagus virus 2 RT with the RT of an ilarvirus from a different subgroup restored the ability to establish persistent infection. These findings open up a new avenue for research on ilarvirus silencing suppression, persistent meristem invasion and vertical transmission.

Utilizing high-resolution ribosome profiling for the global investigation of gene expression in Chlamydomonas

Ribosome profiling (Ribo-seq) is a powerful method for the deep analysis of translation mechanisms and regulatory circuits during gene expression. Extraction and sequencing of ribosome-protected fragments (RPFs) and parallel RNA-seq yields genome-wide insight into translational dynamics and post-transcriptional control of gene expression. Here, we provide details on the Ribo-seq method and the subsequent analysis with the unicellular model alga Chlamydomonas reinhardtii (Chlamydomonas) for generating high-resolution data covering more than 10 000 different transcripts. Detailed analysis of the ribosomal offsets on transcripts uncovers presumable transition states during translocation of elongating ribosomes within the 5' and 3' sections of transcripts and characteristics of eukaryotic translation termination, which are fundamentally distinct for chloroplast translation. In chloroplasts, a heterogeneous RPF size distribution along the coding sequence indicates specific regulatory phases during protein synthesis. For example, local accumulation of small RPFs correlates with local slowdown of psbA translation, possibly uncovering an uncharacterized regulatory step during PsbA/D1 synthesis. Further analyses of RPF distribution along specific cytosolic transcripts revealed characteristic patterns of translation elongation exemplified for the major light-harvesting complex proteins, LHCs. By providing high-quality datasets for all subcellular genomes and attaching our data to the Chlamydomonas reference genome, we aim to make ribosome profiles easily accessible for the broad research community. The data can be browsed without advanced bioinformatic background knowledge for translation output levels of specific genes and their splice variants and for monitoring genome annotation.

Improved super-resolution ribosome profiling reveals prevalent translation of upstream ORFs and small ORFs in Arabidopsis

A crucial step in functional genomics is identifying actively translated ORFs and linking them to biological functions. The challenge lies in identifying short ORFs, as their identification is greatly influenced by data quality and depth. Here, we improved the coverage of super-resolution Ribo-seq in Arabidopsis (Arabidopsis thaliana), revealing uncharacterized translation events for nuclear, chloroplastic, and mitochondrial genes. Assisted by a transcriptome assembly, we identified 7,751 unconventional translation events, comprising 6,996 upstream ORFs (uORFs) and 209 downstream ORFs on annotated protein-coding genes, as well as 546 ORFs in presumed noncoding RNAs. Proteomic data confirmed the production of stable proteins from some of these unannotated translation events. We present evidence of active translation from primary transcripts of trans-acting small interfering RNAs (TAS1-4) and microRNAs (pri-MIR163 and pri-MIR169) and periodic ribosome stalling supporting cotranslational decay. Additionally, we developed a method for identifying extremely short uORFs, including 370 minimum uORFs (AUG-stop), and 2,921 tiny uORFs (2 to 10 amino acids) and 681 uORFs that overlap with each other. Remarkably, these short uORFs exhibit strong translational repression as do longer uORFs. We also systematically discovered 594 uORFs regulated by alternative splicing, suggesting widespread isoform-specific translational control. Finally, these prevalent uORFs are associated with numerous important pathways. In summary, our improved Arabidopsis translational landscape provides valuable resources to study gene expression regulation.

Ribosome Profiling of Plants

Translation is a key step in control of gene expression, yet most analyses of global responses to a stimulus focus on transcription and the transcriptome. For RNA viruses in particular, which have no DNA-templated transcriptional control, control of viral and host translation is crucial. Here, we describe the method of ribosome profiling (ribo-seq) in plants, applied to virus infection. Ribo-seq is a deep sequencing technique that reveals the translatome by presenting a snapshot of the positions and relative amounts of translating ribosomes on all mRNAs in the cell. In contrast to RNA-seq, a crude cell extract is first digested with ribonuclease to degrade all mRNA not protected by a translating 80S ribosome. The resulting ribosome-protected fragments (RPFs) are deep sequenced. The number of reads mapping to a specific mRNA compared to the standard RNA-seq reads reveals the translational efficiency of that mRNA. Moreover, the precise positions of ribosome pause sites, previously unknown translatable open reading frames, and noncanonical translation events can be characterized quantitatively using ribo-seq. As this technique requires meticulous technique, here we present detailed step-by-step instructions for cell lysate preparation by flash freezing of samples, nuclease digestion of cell lysate, monosome collection by sucrose cushion ultracentrifugation, size-selective RNA extraction and rRNA depletion, library preparation for sequencing and finally quality control of sequenced data. These experimental methods apply to many plant systems, with minor nuclease digestion modifications depending on the plant tissue and species. This protocol should be valuable for studies of plant virus gene expression, and the global translational response to virus infection, or any other biotic or abiotic stress, by the host plant.

Genome-Wide Identification, Evolutionary Analysis, and Functional Studies of APX Genes in Melon (Cucuis melo L.)

The antioxidative enzyme ascorbate peroxidase (APX) exerts a critically important function through scavenging reactive oxygen species (ROS), alleviating oxidative damage in plants, and enhancing their tolerance to salinity. Here, we identified 28 CmAPX genes that display an uneven distribution pattern throughout the 12 chromosomes of the melon genome by carrying out a bioinformatics analysis. Phylogenetic analyses revealed that the CmAPX gene family comprised seven different clades, with each clade of genes exhibiting comparable motifs and structures. We cloned 28 CmAPX genes to infer their encoded protein sequences; we then compared these sequences with proteins encoded by rice APX proteins (OsAPX2), Puccinellia tenuiflora APX proteins (PutAPX) and with pea APX proteins. We found that the CmAPX17, CmAPX24, and CmAPX27 genes in Clade I were closely related, and their structures were highly conserved. CmAPX27 (MELO3C020719.2.1) was found to promote resistance to 150 mM NaCl salt stress, according to quantitative real-time fluorescence PCR. Transcriptome data revealed that CmAPX27 was differentially expressed among tissues, and the observed differences in expression were significant. Virus-induced gene silencing of CmAPX27 significantly decreased salinity tolerance, and CmAPX27 exhibited differential expression in the leaf, stem, and root tissues of melon plants. This finding demonstrates that CmAPX27 exerts a key function in melon's tolerance to salt stress. Generally, CmAPX27 could be a target in molecular breeding efforts aimed at improving the salt tolerance of melon; further studies of CmAPX27 could unveil novel physiological mechanisms through which antioxidant enzymes mitigate the deleterious effects of ROS stress.

The distinct translational landscapes of gram-negative Salmonella and gram-positive Listeria

Translational control in pathogenic bacteria is fundamental to gene expression and affects virulence and other infection phenotypes. We used an enhanced ribosome profiling protocol coupled with parallel transcriptomics to capture accurately the global translatome of two evolutionarily distant pathogenic bacteria-the Gram-negative bacterium Salmonella and the Gram-positive bacterium Listeria. We find that the two bacteria use different mechanisms to translationally regulate protein synthesis. In Salmonella, in addition to the expected correlation between translational efficiency and cis-regulatory features such as Shine-Dalgarno (SD) strength and RNA secondary structure around the initiation codon, our data reveal an effect of the 2nd and 3rd codons, where the presence of tandem lysine codons (AAA-AAA) enhances translation in both Salmonella and E. coli. Strikingly, none of these features are seen in efficiently translated Listeria transcripts. Instead, approximately 20% of efficiently translated Listeria genes exhibit 70 S footprints seven nt upstream of the authentic start codon, suggesting that these genes may be subject to a novel translational initiation mechanism. Our results show that SD strength is not a direct hallmark of translational efficiency in all bacteria. Instead, Listeria has evolved additional mechanisms to control gene expression level that are distinct from those utilised by Salmonella and E. coli.

What Can Ribo-Seq, Immunopeptidomics, and Proteomics Tell Us About the Noncanonical Proteome?

Ribosome profiling (Ribo-Seq) has proven transformative for our understanding of the human genome and proteome by illuminating thousands of noncanonical sites of ribosome translation outside the currently annotated coding sequences (CDSs). A conservative estimate suggests that at least 7000 noncanonical ORFs are translated, which, at first glance, has the potential to expand the number of human protein CDSs by 30%, from ∼19,500 annotated CDSs to over 26,000 annotated CDSs. Yet, additional scrutiny of these ORFs has raised numerous questions about what fraction of them truly produce a protein product and what fraction of those can be understood as proteins according to conventional understanding of the term. Adding further complication is the fact that published estimates of noncanonical ORFs vary widely by around 30-fold, from several thousand to several hundred thousand. The summation of this research has left the genomics and proteomics communities both excited by the prospect of new coding regions in the human genome but searching for guidance on how to proceed. Here, we discuss the current state of noncanonical ORF research, databases, and interpretation, focusing on how to assess whether a given ORF can be said to be "protein coding."

What can Ribo-seq and proteomics tell us about the non-canonical proteome?

Ribosome profiling (Ribo-seq) has proven transformative for our understanding of the human genome and proteome by illuminating thousands of non-canonical sites of ribosome translation outside of the currently annotated coding sequences (CDSs). A conservative estimate suggests that at least 7,000 non-canonical open reading frames (ORFs) are translated, which, at first glance, has the potential to expand the number of human protein-coding sequences by 30%, from ∼19,500 annotated CDSs to over 26,000. Yet, additional scrutiny of these ORFs has raised numerous questions about what fraction of them truly produce a protein product and what fraction of those can be understood as proteins according to conventional understanding of the term. Adding further complication is the fact that published estimates of non-canonical ORFs vary widely by around 30-fold, from several thousand to several hundred thousand. The summation of this research has left the genomics and proteomics communities both excited by the prospect of new coding regions in the human genome, but searching for guidance on how to proceed. Here, we discuss the current state of non-canonical ORF research, databases, and interpretation, focusing on how to assess whether a given ORF can be said to be "protein-coding".

In brief

The human genome encodes thousands of non-canonical open reading frames (ORFs) in addition to protein-coding genes. As a nascent field, many questions remain regarding non-canonical ORFs. How many exist? Do they encode proteins? What level of evidence is needed for their verification? Central to these debates has been the advent of ribosome profiling (Ribo-seq) as a method to discern genome-wide ribosome occupancy, and immunopeptidomics as a method to detect peptides that are processed and presented by MHC molecules and not observed in traditional proteomics experiments. This article provides a synthesis of the current state of non-canonical ORF research and proposes standards for their future investigation and reporting.

Highlights

Combined use of Ribo-seq and proteomics-based methods enables optimal confidence in detecting non-canonical ORFs and their protein products.Ribo-seq can provide more sensitive detection of non-canonical ORFs, but data quality and analytical pipelines will impact results.Non-canonical ORF catalogs are diverse and span both high-stringency and low-stringency ORF nominations.A framework for standardized non-canonical ORF evidence will advance the research field.

Ribosome footprint profiling enables elucidating the systemic regulation of fatty acid accumulation in Acer truncatum

BACKGROUND

The accumulation of fatty acids in plants covers a wide range of functions in plant physiology and thereby affects adaptations and characteristics of species. As the famous woody oilseed crop, Acer truncatum accumulates unsaturated fatty acids and could serve as the model to understand the regulation and trait formation in oil-accumulation crops. Here, we performed Ribosome footprint profiling combing with a multi-omics strategy towards vital time points during seed development, and finally constructed systematic profiling from transcription to proteomes. Additionally, we characterized the small open reading frames (ORFs) and revealed that the translational efficiencies of focused genes were highly influenced by their sequence features.

RESULTS

The comprehensive multi-omics analysis of lipid metabolism was conducted in A. truncatum. We applied the Ribo-seq and RNA-seq techniques, and the analyses of transcriptional and translational profiles of seeds collected at 85 and 115 DAF were compared. Key members of biosynthesis-related structural genes (LACS, FAD2, FAD3, and KCS) were characterized fully. More meaningfully, the regulators (MYB, ABI, bZIP, and Dof) were identified and revealed to affect lipid biosynthesis via post-translational regulations. The translational features results showed that translation efficiency tended to be lower for the genes with a translated uORF than for the genes with a non-translated uORF. They provide new insights into the global mechanisms underlying the developmental regulation of lipid metabolism.

CONCLUSIONS

We performed Ribosome footprint profiling combing with a multi-omics strategy in A. truncatum seed development, which provides an example of the use of Ribosome footprint profiling in deciphering the complex regulation network and will be useful for elucidating the metabolism of A. truncatum seed oil and the regulatory mechanisms.

Genome-wide identification, evolution, and expression analysis of MLO gene family in melon (Cucumis melo L.)

Powdery mildew (PM) is one of the main fungal diseases that appear during the cultivation of the melon fruit crop. Mildew Resistance Locus "O" (MLO) is known as a gene family and has seven conserved transmembrane domains. An induced functional loss of a specific MLO gene could mainly confer PM resistance to melons. However, the genomic structure of MLO genes and its main role in PM resistance still remain unclear in melon. In this study, bioinformatic analysis identified a total of 14 MLO gene family members in the melon genome sequence, and these genes were distributed in an uneven manner on eight chromosomes. The phylogenetic analysis divided the CmMLO genes into five different clades, and gene structural analysis showed that genes in the same clade had similar intron and exon distribution patterns. In addition, by cloning the CmMLO gene sequence in four melon lines, analyzing the CmMLO gene expression pattern after infection, and making microscopic observations of the infection pattern of PM, we concluded that the CmMLO5 (MELO3C012438) gene plays a negative role in regulating PM-resistance in the susceptible melon line (Topmark), and the critical time point for gene function was noticed at 24 and 72 hours after PM infection. The mutational analysis exhibited a single base mutation at 572 bp, which further results in loss of protein function, thus conferring PM resistance in melon. In summary, our research evidence provides a thorough understanding of the CmMLO gene family and demonstrates their potential role in disease resistance, as well as a theoretical foundation for melon disease resistance breeding.

Discovering microproteins: making the most of ribosome profiling data

Building a reference set of protein-coding open reading frames (ORFs) has revolutionized biological process discovery and understanding. Traditionally, gene models have been confirmed using cDNA sequencing and encoded translated regions inferred using sequence-based detection of start and stop combinations longer than 100 amino-acids to prevent false positives. This has led to small ORFs (smORFs) and their encoded proteins left un-annotated. Ribo-seq allows deciphering translated regions from untranslated irrespective of the length. In this review, we describe the power of Ribo-seq data in detection of smORFs while discussing the major challenge posed by data-quality, -depth and -sparseness in identifying the start and end of smORF translation. In particular, we outline smORF cataloguing efforts in humans and the large differences that have arisen due to variation in data, methods and assumptions. Although current versions of smORF reference sets can already be used as a powerful tool for hypothesis generation, we recommend that future editions should consider these data limitations and adopt unified processing for the community to establish a canonical catalogue of translated smORFs.

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.

Identification and analysis of smORFs in Chlamydomonas reinhardtii

Small open reading frames (smORFs) have been acknowledged as an important partner in organism functions ranging from bacteria to higher eukaryotes. However, lack of investigation of smORFs in green algae, despite their importance in ecology and evolution. We applied bioinformatic analysis, ribosome profiling, and small peptide proteomics to provide a genome-wide and high-confident smORF database in the model green alga Chlamydomonas reinhardtii. The whole genome was screened first to mine potential coding smORFs. Then conservative analysis, ribosome profiling, and proteomics data were processed to identify conserved smORFs and generate translation evidence. The combination of procedures resulted in 2014 smORFs that might exist in the C. reinhardtii genome. The expression of smORFs in Cd treatment suggested that two smORFs might participate in redox reaction, three in inorganic phosphate transport, and one in DNA repair under stress. Our study built a genome-widely database in C. reinhardtii, providing target smORFs for further research.

Ciliogenesis requires sphingolipid-dependent membrane and axoneme interaction

Cilium formation and regeneration requires new protein synthesis, but the underlying cytosolic translational reprogramming remains largely unknown. Using ribosome footprinting, we performed global translatome profiling during cilia regeneration in Chlamydomonas and uncovered that flagellar genes undergo an early transcriptional activation but late translational repression. This pattern guided our identification of sphingolipid metabolism enzymes, including serine palmitoyltransferase (SPT), as essential regulators for ciliogenesis. Cryo-electron tomography showed that ceramide loss abnormally increased the membrane-axoneme distance and generated bulged cilia. We found that ceramides interact with intraflagellar transport (IFT) particle proteins that IFT motors transport along axoneme microtubules (MTs), suggesting that ceramide-IFT particle-IFT motor-MT interactions connect the ciliary membrane with the axoneme to form rod-shaped cilia. SPT-deficient vertebrate cells were defective in ciliogenesis, and SPT mutations from patients with hereditary sensory neuropathy disrupted cilia, which could be restored by sphingolipid supplementation. These results reveal a conserved role of sphingolipid in cilium formation and link compromised sphingolipid production with ciliopathies.

Comparison of rRNA depletion methods for efficient bacterial mRNA sequencing

Current methods of high-throughput RNA sequencing of prokaryotes, including transcriptome analysis or ribosomal profiling, need deep sequencing to achieve sufficient numbers of effective reads (e.g., mapping to mRNA) in order to also find weakly expressed genetic elements. The fraction of high-quality reads mapping to coding RNAs (i.e., mRNA) is mainly influenced by the large content of rRNA and, to a lesser extent, tRNA in total RNA. Thus, depletion of rRNA increases coverage and thus sequencing costs. RiboZero, a depletion kit based on probe hybridisation and rRNA-removal was found to be most efficient in the past, but it was discontinued in 2018. To facilitate comparability with previous experiments and to help choose adequate replacements, we compare three commercially available rRNA depletion kits also based on hybridization and magnetic beads, i.e., riboPOOLs, RiboMinus and MICROBExpress, with the former RiboZero. Additionally, we constructed biotinylated probes for magnetic bead capture and rRNA depletion in this study. Based on E. coli, we found similar efficiencies in rRNA depletion for riboPOOLs and the self-made depletion method; both comparable to the former RiboZero, followed by RiboMinus, succeeded by MICROBExpress. Further, our in-house protocol allows customized species-specific rRNA or even tRNA depletion or depletion of other RNA targets. Both, the self-made biotinylated probes and riboPOOLs, were most successful in reducing the rRNA content and thereby increasing sequencing depth concerning mRNA reads. Additionally, the number of reads matching to weakly expressed genes are increased. In conclusion, the self-made specific biotinylated probes and riboPOOLs are an adequate replacement for the former RiboZero. Both are very efficient in depleting rRNAs, increasing mRNA reads and thus sequencing efficiency.

Short open reading frames (sORFs) and microproteins: an update on their identification and validation measures

A short open reading frame (sORFs) constitutes ≤ 300 bases, encoding a microprotein or sORF-encoded protein (SEP) which comprises ≤ 100 amino acids. Traditionally dismissed by genome annotation pipelines as meaningless noise, sORFs were found to possess coding potential with ribosome profiling (RIBO-Seq), which unveiled sORF-based transcripts at various genome locations. Nonetheless, the existence of corresponding microproteins that are stable and functional was little substantiated by experimental evidence initially. With recent advancements in multi-omics, the identification, validation, and functional characterisation of sORFs and microproteins have become feasible. In this review, we discuss the history and development of an emerging research field of sORFs and microproteins. In particular, we focus on an array of bioinformatics and OMICS approaches used for predicting, sequencing, validating, and characterizing these recently discovered entities. These strategies include RIBO-Seq which detects sORF transcripts via ribosome footprints, and mass spectrometry (MS)-based proteomics for sequencing the resultant microproteins. Subsequently, our discussion extends to the functional characterisation of microproteins by incorporating CRISPR/Cas9 screen and protein–protein interaction (PPI) studies. Our review discusses not only detection methodologies, but we also highlight on the challenges and potential solutions in identifying and validating sORFs and their microproteins. The novelty of this review lies within its validation for the functional role of microproteins, which could contribute towards the future landscape of microproteomics.

Ribosome profiling of porcine reproductive and respiratory syndrome virus reveals novel features of viral gene expression

The arterivirus porcine reproductive and respiratory syndrome virus (PRRSV) causes significant economic losses to the swine industry worldwide. Here we apply ribosome profiling (RiboSeq) and parallel RNA sequencing (RNASeq) to characterise the transcriptome and translatome of both species of PRRSV and to analyse the host response to infection. We calculated programmed ribosomal frameshift (PRF) efficiency at both sites on the viral genome. This revealed the nsp2 PRF site as the second known example where temporally regulated frameshifting occurs, with increasing -2 PRF efficiency likely facilitated by accumulation of the PRF-stimulatory viral protein, nsp1β. Surprisingly, we find that PRF efficiency at the canonical ORF1ab frameshift site also increases over time, in contradiction of the common assumption that RNA structure-directed frameshift sites operate at a fixed efficiency. This has potential implications for the numerous other viruses with canonical PRF sites. Furthermore, we discovered several highly translated additional viral ORFs, the translation of which may be facilitated by multiple novel viral transcripts. For example, we found a highly expressed 125-codon ORF overlapping nsp12, which is likely translated from novel subgenomic RNA transcripts that overlap the 3' end of ORF1b. Similar transcripts were discovered for both PRRSV-1 and PRRSV-2, suggesting a potential conserved mechanism for temporally regulating expression of the 3'-proximal region of ORF1b. We also identified a highly translated, short upstream ORF in the 5' UTR, the presence of which is highly conserved amongst PRRSV-2 isolates. These findings reveal hidden complexity in the gene expression programmes of these important nidoviruses.

[Ribosomal profiling as a tool for studying translation in plants: main results, problems and future prospects]

The expression of eukaryotic genes can be regulated at several stages, including the translation of mRNA. It is known that the structure of mRNA can affect both the efficiency of interaction with the translation apparatus in general and the choice of translation initiation sites. To study the translated fraction of the transcriptome, experimental methods of analysis were developed, the most informative of which is ribosomal profiling (RP, Ribo-seq). Originally developed for use in yeast systems, this method has been adapted for research in translation mechanisms in many plant species. This technology includes the isolation of the polysomal fraction and high-performance sequencing of a pool of mRNA fragments associated with ribosomes. Comparing the results of transcript coverage with reads obtained using the ribosome profiling with the transcriptional efficiency of genes allows the translation efficiency to be evaluated for each transcript. The exact positions of ribosomes determined on mRNA sequences allow determining the translation of open reading frames and switching between the translation of several reading frames - a phenomenon in which two or more overlapping frames are read from one mRNA and different proteins are synthesized. The advantage of this method is that it provides quantitative estimates of ribosome coverage of mRNA and can detect relatively rare translation events. Using this technology, it was possible to identify and classify plant genes by the type of regulation of their expression at the transcription, translation, or both levels. Features of the mRNA structure that affect translation levels have been revealed: the formation of G2 quadruplexes and the presence of specific motifs in the 5'-UTR region, GC content, the presence of alternative translation starts, and the influence of uORFs on the translation of downstream mORFs. In this review, we briefly reviewed the RP methodology and the prospects for its application to study the structural and functional organization and regulation of plant gene expression.

Translation complex stabilization on messenger RNA and footprint profiling to study the RNA responses and dynamics of protein biosynthesis in the cells

During protein biosynthesis, ribosomes bind to messenger (m)RNA, locate its protein-coding information, and translate the nucleotide triplets sequentially as codons into the corresponding sequence of amino acids, forming proteins. Non-coding mRNA features, such as 5' and 3' untranslated regions (UTRs), start sites or stop codons of different efficiency, stretches of slower or faster code and nascent polypeptide interactions can alter the translation rates transcript-wise. Most of the homeostatic and signal response pathways of the cells converge on individual mRNA control, as well as alter the global translation output. Among the multitude of approaches to study translational control, one of the most powerful is to infer the locations of translational complexes on mRNA based on the mRNA fragments protected by these complexes from endonucleolytic hydrolysis, or footprints. Translation complex profiling by high-throughput sequencing of the footprints allows to quantify the transcript-wise, as well as global, alterations of translation, and uncover the underlying control mechanisms by attributing footprint locations and sizes to different configurations of the translational complexes. The accuracy of all footprint profiling approaches critically depends on the fidelity of footprint generation and many methods have emerged to preserve certain or multiple configurations of the translational complexes, often in challenging biological material. In this review, a systematic summary of approaches to stabilize translational complexes on mRNA for footprinting is presented and major findings are discussed. Future directions of translation footprint profiling are outlined, focusing on the fidelity and accuracy of inference of the native in vivo translation complex distribution on mRNA.

Investigating molecular mechanisms of 2A-stimulated ribosomal pausing and frameshifting in Theilovirus

The 2A protein of Theiler's murine encephalomyelitis virus (TMEV) acts as a switch to stimulate programmed -1 ribosomal frameshifting (PRF) during infection. Here, we present the X-ray crystal structure of TMEV 2A and define how it recognises the stimulatory RNA element. We demonstrate a critical role for bases upstream of the originally predicted stem-loop, providing evidence for a pseudoknot-like conformation and suggesting that the recognition of this pseudoknot by beta-shell proteins is a conserved feature in cardioviruses. Through examination of PRF in TMEV-infected cells by ribosome profiling, we identify a series of ribosomal pauses around the site of PRF induced by the 2A-pseudoknot complex. Careful normalisation of ribosomal profiling data with a 2A knockout virus facilitated the identification, through disome analysis, of ribosome stacking at the TMEV frameshifting signal. These experiments provide unparalleled detail of the molecular mechanisms underpinning Theilovirus protein-stimulated frameshifting.

Simultaneous ribosome profiling of hundreds of microbes from the human microbiome

Ribosome profiling enables sequencing of ribosome-bound fragments of RNA, revealing which transcripts are being translated as well as the position of ribosomes along mRNAs. Although ribosome profiling has been applied to cultured bacterial isolates, its application to uncultured, mixed communities has been challenging. We present MetaRibo-Seq, a protocol that enables the application of ribosome profiling directly to the human fecal microbiome. MetaRibo-Seq is a benchmarked method that includes several modifications to existing ribosome profiling protocols, specifically addressing challenges involving fecal sample storage, purity and input requirements. We also provide a computational workflow to quality control and trim reads, de novo assemble a reference metagenome with metagenomic reads, align MetaRibo-Seq reads to the reference, and assess MetaRibo-Seq library quality ( https://github.com/bhattlab/bhattlab_workflows/tree/master/metariboseq ). This MetaRibo-Seq protocol enables researchers in standard molecular biology laboratories to study translation in the fecal microbiome in ~5 d.

Small-scale sequencing enables quality assessment of Ribo-Seq data: an example from Arabidopsis cell culture

BACKGROUND

Translation is a tightly regulated process, controlling the rate of protein synthesis in cells. Ribosome sequencing (Ribo-Seq) is a recently developed tool for studying actively translated mRNA and can thus directly address translational regulation. Ribo-Seq libraries need to be sequenced to a great depth due to high contamination by rRNA and other contaminating nucleic acid fragments. Deep sequencing is expensive, and it generates large volumes of data, making data analysis complicated and time consuming.

METHODS AND RESULTS

Here we developed a platform for Ribo-Seq library construction and data analysis to enable rapid quality assessment of Ribo-Seq libraries with the help of a small-scale sequencer. Our data show that several qualitative features of a Ribo-Seq library, such as read length distribution, P-site distribution, reading frame and triplet periodicity, can be effectively evaluated using only the data generated by a benchtop sequencer with a very limited number of reads.

CONCLUSION

Our pipeline enables rapid evaluation of Ribo-Seq libraries, opening up possibilities for optimization of Ribo-Seq library construction from difficult samples, and leading to better decision making prior to more costly deep sequencing.

HRIBO: high-throughput analysis of bacterial ribosome profiling data

MOTIVATION

Ribosome profiling (Ribo-seq) is a powerful approach based on deep sequencing of cDNA libraries generated from ribosome-protected RNA fragments to explore the translatome of a cell, and is especially useful for the detection of small proteins (50-100 amino acids) that are recalcitrant to many standard biochemical and in silico approaches. While pipelines are available to analyze Ribo-seq data, none are designed explicitly for the automatic processing and analysis of data from bacteria, nor are they focused on the discovery of unannotated open reading frames (ORFs).

RESULTS

We present HRIBO (High-throughput annotation by Ribo-seq), a workflow to enable reproducible and high-throughput analysis of bacterial Ribo-seq data. The workflow performs all required pre-processing and quality control steps. Importantly, HRIBO outputs annotation-independent ORF predictions based on two complementary bacteria-focused tools, and integrates them with additional feature information and expression values. This facilitates the rapid and high-confidence discovery of novel ORFs and their prioritization for functional characterization.

AVAILABILITY AND IMPLEMENTATION

HRIBO is a free and open source project available under the GPL-3 license at: https://github.com/RickGelhausen/HRIBO.

Chaperone-mediated coupling of subunit availability to activation of flagellar Type III secretion

Bacterial flagellar subunits are exported across the cell membrane by the flagellar Type III Secretion System (fT3SS), powered by the proton motive force (pmf) and a specialized ATPase that enables the flagellar export gate to utilize the pmf electric potential (ΔΨ). Export gate activation is mediated by the ATPase stalk, FliJ, but how this process is regulated to prevent wasteful dissipation of pmf in the absence of subunit cargo is not known. Here, we show that FliJ activation of the export gate is regulated by flagellar export chaperones. FliJ binds unladen chaperones and, by using novel chaperone variants specifically defective for FliJ binding, we show that disruption of this interaction attenuates motility and cognate subunit export. We demonstrate in vitro that chaperones and the FlhA export gate component compete for binding to FliJ, and show in vivo that unladen chaperones, which would be present in the cell when subunit levels are low, sequester FliJ to prevent activation of the export gate and attenuate subunit export. Our data indicate a mechanism whereby chaperones couple availability of subunit cargo to pmf-driven export by the fT3SS.

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 .

Manipulation of the unfolded protein response: A pharmacological strategy against coronavirus infection

Coronavirus infection induces the unfolded protein response (UPR), a cellular signalling pathway composed of three branches, triggered by unfolded proteins in the endoplasmic reticulum (ER) due to high ER load. We have used RNA sequencing and ribosome profiling to investigate holistically the transcriptional and translational response to cellular infection by murine hepatitis virus (MHV), often used as a model for the Betacoronavirus genus to which the recently emerged SARS-CoV-2 also belongs. We found the UPR to be amongst the most significantly up-regulated pathways in response to MHV infection. To confirm and extend these observations, we show experimentally the induction of all three branches of the UPR in both MHV- and SARS-CoV-2-infected cells. Over-expression of the SARS-CoV-2 ORF8 or S proteins alone is itself sufficient to induce the UPR. Remarkably, pharmacological inhibition of the UPR greatly reduced the replication of both MHV and SARS-CoV-2, revealing the importance of this pathway for successful coronavirus replication. This was particularly striking when both IRE1α and ATF6 branches of the UPR were inhibited, reducing SARS-CoV-2 virion release (~1,000-fold). Together, these data highlight the UPR as a promising antiviral target to combat coronavirus infection.

High-throughput 5'P sequencing enables the study of degradation-associated ribosome stalls

RNA degradation is critical for gene expression and mRNA quality control. mRNA degradation is connected to the translation process up to the degree that 5'-3' mRNA degradation follows the last translating ribosome. Here, we present an improved high-throughput 5'P degradome RNA-sequencing method (HT-5Pseq). HT-5Pseq is easy, scalable, and uses affordable duplex-specific nuclease-based rRNA depletion. We investigate in vivo ribosome stalls focusing on translation termination. By comparing ribosome stalls identified by ribosome profiling, disome-seq and HT-5Pseq, we find that degradation-associated ribosome stalls are often enriched in Arg preceding the stop codon. On the contrary, mRNAs depleted for those stalls use more frequently a TAA stop codon preceded by hydrophobic amino acids. Finally, we show that termination stalls found by HT-5Pseq, and not by other approaches, are associated with decreased mRNA stability. Our work suggests that ribosome stalls associated with mRNA decay can be easily captured by investigating the 5'P degradome.

Ribo-ODDR: Oligo design pipeline for experiment-specific rRNA depletion in ribo-seq

Motivation

Ribosome Profiling (Ribo-seq) has revolutionized the study of RNA translation by providing information on ribosome positions across all translated RNAs with nucleotide-resolution. Yet several technical limitations restrict the sequencing depth of such experiments, the most common of which is the overabundance of rRNA fragments. Various strategies can be employed to tackle this issue, including the use of commercial rRNA depletion kits. However, as they are designed for more standardized RNAseq experiments, they may perform suboptimally in Ribo-seq. In order to overcome this, it is possible to use custom biotinylated oligos complementary to the most abundant rRNA fragments, however currently no computational framework exists to aid the design of optimal oligos.

Results

Here, we first show that a major confounding issue is that the rRNA fragments generated via Ribo-seq vary significantly with differing experimental conditions, suggesting that a ‘one-size-fits-all’ approach may be inefficient. Therefore we developed Ribo-ODDR, an oligo design pipeline integrated with a user-friendly interface that assists in oligo selection for efficient experiment-specific rRNA depletion. Ribo-ODDR uses preliminary data to identify the most abundant rRNA fragments, and calculates the rRNA depletion efficiency of potential oligos. We experimentally show that Ribo-ODDR designed oligos outperform commercially available kits and lead to a significant increase in rRNA depletion in Ribo-seq.

Availability and implementation

Ribo-ODDR is freely accessible at https://github.com/fallerlab/Ribo-ODDR.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

RiboToolkit: an integrated platform for analysis and annotation of ribosome profiling data to decode mRNA translation at codon resolution

Ribosome profiling (Ribo-seq) is a powerful technology for globally monitoring RNA translation; ranging from codon occupancy profiling, identification of actively translated open reading frames (ORFs), to the quantification of translational efficiency under various physiological or experimental conditions. However, analyzing and decoding translation information from Ribo-seq data is not trivial. Although there are many existing tools to analyze Ribo-seq data, most of these tools are designed for specific or limited functionalities and an easy-to-use integrated tool to analyze Ribo-seq data is lacking. Fortunately, the small size (26–34 nt) of ribosome protected fragments (RPFs) in Ribo-seq and the relatively small amount of sequencing data greatly facilitates the development of such a web platform, which is easy to manipulate for users with or without bioinformatic expertise. Thus, we developed RiboToolkit (http://rnabioinfor.tch.harvard.edu/RiboToolkit), a convenient, freely available, web-based service to centralize Ribo-seq data analyses, including data cleaning and quality evaluation, expression analysis based on RPFs, codon occupancy, translation efficiency analysis, differential translation analysis, functional annotation, translation metagene analysis, and identification of actively translated ORFs. Besides, easy-to-use web interfaces were developed to facilitate data analysis and intuitively visualize results. Thus, RiboToolkit will greatly facilitate the study of mRNA translation based on ribosome profiling.

Nuclease-mediated depletion biases in ribosome footprint profiling libraries

Ribosome footprint profiling is a high-throughput sequencing-based technique that provides detailed and global views of translation in living cells. An essential part of this technology is removal of unwanted, normally very abundant, ribosomal RNA sequences that dominate libraries and increase sequencing costs. The most effective commercial solution (Ribo-Zero) has been discontinued as a standalone product and a number of new, experimentally distinct commercial applications have emerged on the market. Here we evaluated several commercially available alternatives designed for RNA-seq of human samples and find them generally unsuitable for ribosome footprint profiling. We instead recommend the use of custom-designed biotinylated oligos, which were widely used in early ribosome profiling studies. Importantly, we warn that depletion solutions based on targeted nuclease cleavage significantly perturb the high-resolution information that can be derived from the data, and thus do not recommend their use for any applications that require precise determination of the ends of RNA fragments.

Transcriptomic and Translatomic Analyses Reveal Insights into the Developmental Regulation of Secondary Metabolism in the Young Shoots of Tea Plants (Camellia sinensis L.)

Accumulation of secondary metabolites in the young shoots of tea plants is developmentally modulated, especially flavonoids. Here, we investigate the developmental regulation mechanism of secondary metabolism in the developing leaves of tea plants using an integrated multiomic approach. For the pair of Leaf2/Bud, the correlation coefficient of the fold change of mRNA and RPFs abundances involved in flavonoid biosynthesis was 0.9359, being higher than that of RPFs and protein (R² = 0.6941). These correlations were higher than the corresponding correlation coefficients for secondary metabolisms and genome-wide scale. Metabolomic analysis demonstrates that the developmental modulations of the structural genes for flavonoid biosynthesis-related pathways align with the concentration changes of catechin and flavonol glycoside groups. Relatively high translational efficiency (TE > 2) was observed in the four flavonoid structural genes (chalcone isomerase, dihydroflavonol 4-reductase, anthocyanidin synthase, and flavonol synthase). In addition, we originally provided the information on identified small open reading frames (small ORFs) and main ORFs in tea leaves and elaborated that the presence of upstream ORFs may have a repressive effect on the translation of downstream ORFs. Our data suggest that transcriptional regulation coordinates with translational regulation and may contribute to the elevation of translational efficiencies for the structural genes involved in the flavonoid biosynthesis pathways during tea leaf development.

Ribosome profiling in plants: what is not lost in translation?

Translation is a highly dynamic cellular process whereby genetic information residing in an mRNA molecule is converted into a protein that in turn executes specific functions. However, pre-synthesized mRNA levels do not always correlate with corresponding protein levels, suggesting that translational control plays an essential role in gene regulation. A better understanding of how gene expression is regulated during translation will enable the discovery of new genes and mechanisms that control important traits in plants. Therefore, in recent years, several methods have been developed to analyse the translatome; that is, all mRNAs being actively translated at a given time, tissue, and/or developmental stage. Ribosome profiling or ribo-seq is one such technology revolutionizing our ability to analyse the translatome and in turn understand translational control of gene expression. Ribo-seq involves isolating mRNA-ribosome complexes, treating them with a RNase, and then identifying ribosome-protected mRNA regions by deep sequencing. Here, we briefly review recent ribosome profiling studies that revealed new insights into plant biology. Manipulation of novel genes identified using ribosome profiling could prove useful for increasing yield through improved biotic and abiotic stress tolerance.

Selective Translation of Cell Fate Regulators Mediates Tolerance to Broad Oncogenic Stress

Human skin tolerates a surprisingly high burden of oncogenic lesions. Although adult epidermis can suppress the expansion of individual mutant clones, the mechanisms behind tolerance to oncogene activation across broader regions of tissue are unclear. Here, we uncover a dynamic translational mechanism that coordinates oncogenic HRAS-induced hyperproliferation with loss of progenitor self-renewal to restrain aberrant growth and tumorigenesis. We identify translation initiator eIF2B5 as a central co-regulator of HRAS proliferation and cell fate choice. By coupling in vivo ribosome profiling with genetic screening, we provide direct evidence that oncogene-induced loss of progenitor self-renewal is driven by eIF2B5-mediated translation of ubiquitination genes. Ubiquitin ligase FBXO32 specifically inhibits epidermal renewal without affecting overall proliferation, thus restraining HRAS-driven tumorigenesis while maintaining normal tissue growth. Thus, oncogene-driven translation is not necessarily inherently tumor promoting but instead can manage widespread oncogenic stress by steering progenitor fate to prolong normal tissue growth.

RiboMiner: a toolset for mining multi-dimensional features of the translatome with ribosome profiling data

Background

Ribosome profiling has been widely used for studies of translation under a large variety of cellular and physiological contexts. Many of these studies have greatly benefitted from a series of data-mining tools designed for dissection of the translatome from different aspects. However, as the studies of translation advance quickly, the current toolbox still falls in short, and more specialized tools are in urgent need for deeper and more efficient mining of the important and new features of the translation landscapes.

Results

Here, we present RiboMiner, a bioinformatics toolset for mining of multi-dimensional features of the translatome with ribosome profiling data. RiboMiner performs extensive quality assessment of the data and integrates a spectrum of tools for various metagene analyses of the ribosome footprints and for detailed analyses of multiple features related to translation regulation. Visualizations of all the results are available. Many of these analyses have not been provided by previous methods. RiboMiner is highly flexible, as the pipeline could be easily adapted and customized for different scopes and targets of the studies.

Conclusions

Applications of RiboMiner on two published datasets did not only reproduced the main results reported before, but also generated novel insights into the translation regulation processes. Therefore, being complementary to the current tools, RiboMiner could be a valuable resource for dissections of the translation landscapes and the translation regulations by mining the ribosome profiling data more comprehensively and with higher resolution. RiboMiner is freely available at https://github.com/xryanglab/RiboMiner and https://pypi.org/project/RiboMiner.

Control of translation during the unfolded protein response in maize seedlings: Life without PERKs

The accumulation of misfolded proteins in the endoplasmic reticulum (ER) defines a condition called ER stress that induces the unfolded protein response (UPR). The UPR in mammalian cells attenuates protein synthesis initiation, which prevents the piling up of misfolded proteins in the ER. Mammalian cells rely on Protein Kinase RNA-Like Endoplasmic Reticulum Kinase (PERK) phosphorylation of eIF2α to arrest protein synthesis, however, plants do not have a PERK homolog, so the question is whether plants control translation in response to ER stress. We compared changes in RNA levels in the transcriptome to the RNA levels protected by ribosomes and found a decline in translation efficiency, including many UPR genes, in response to ER stress. The decline in translation efficiency is due to the fact that many mRNAs are not loaded onto polyribosomes (polysomes) in proportion to their increase in total RNA, instead some of the transcripts accumulate in stress granules (SGs). The RNAs that populate SGs are not derived from the disassembly of polysomes because protein synthesis remains steady during stress. Thus, the surge in transcription of UPR genes in response to ER stress is accompanied by the formation of SGs, and the sequestration of mRNAs in SGs may serve to temporarily relieve the translation load during ER stress.

An RNA thermoswitch regulates daytime growth in Arabidopsis

Temperature is a major environmental cue affecting plant growth and development. Plants often experience higher temperatures in the context of a 24 h day-night cycle, with temperatures peaking in the middle of the day. Here, we find that the transcript encoding the bHLH transcription factor PIF7 undergoes a direct increase in translation in response to warmer temperature. Diurnal expression of PIF7 transcript gates this response, allowing PIF7 protein to quickly accumulate in response to warm daytime temperature. Enhanced PIF7 protein levels directly activate the thermomorphogenesis pathway by inducing the transcription of key genes such as the auxin biosynthetic gene YUCCA8, and are necessary for thermomorphogenesis to occur under warm cycling daytime temperatures. The temperature-dependent translational enhancement of PIF7 messenger RNA is mediated by the formation of an RNA hairpin within its 5' untranslated region, which adopts an alternative conformation at higher temperature, leading to increased protein synthesis. We identified similar hairpin sequences that control translation in additional transcripts including WRKY22 and the key heat shock regulator HSFA2, suggesting that this is a conserved mechanism enabling plants to respond and adapt rapidly to high temperatures.

Computational methods for ribosome profiling data analysis

Since the introduction of the ribosome profiling technique in 2009 its popularity has greatly increased. It is widely used for the comprehensive assessment of gene expression and for studying the mechanisms of regulation at the translational level. As the number of ribosome profiling datasets being produced continues to grow, so too does the need for reliable software that can provide answers to the biological questions it can address. This review describes the computational methods and tools that have been developed to analyze ribosome profiling data at the different stages of the process. It starts with initial routine processing of raw data and follows with more specific tasks such as the identification of translated open reading frames, differential gene expression analysis, or evaluation of local or global codon decoding rates. The review pinpoints challenges associated with each step and explains the ways in which they are currently addressed. In addition it provides a comprehensive, albeit incomplete, list of publicly available software applicable to each step, which may be a beneficial starting point to those unexposed to ribosome profiling analysis. The outline of current challenges in ribosome profiling data analysis may inspire computational biologists to search for novel, potentially superior, solutions that will improve and expand the bioinformatician's toolbox for ribosome profiling data analysis. This article is characterized under: Translation > Ribosome Structure/Function RNA Evolution and Genomics > Computational Analyses of RNA Translation > Translation Mechanisms Translation > Translation Regulation.

RiboFlow, RiboR and RiboPy: an ecosystem for analyzing ribosome profiling data at read length resolution

SUMMARY

Ribosome occupancy measurements enable protein abundance estimation and infer mechanisms of translation. Recent studies have revealed that sequence read lengths in ribosome profiling data are highly variable and carry critical information. Consequently, data analyses require the computation and storage of multiple metrics for a wide range of ribosome footprint lengths. We developed a software ecosystem including a new efficient binary file format named 'ribo'. Ribo files store all essential data grouped by ribosome footprint lengths. Users can assemble ribo files using our RiboFlow pipeline that processes raw ribosomal profiling sequencing data. RiboFlow is highly portable and customizable across a large number of computational environments with built-in capabilities for parallelization. We also developed interfaces for writing and reading ribo files in the R (RiboR) and Python (RiboPy) environments. Using RiboR and RiboPy, users can efficiently access ribosome profiling quality control metrics, generate essential plots and carry out analyses. Altogether, these components create a software ecosystem for researchers to study translation through ribosome profiling.

AVAILABILITY AND IMPLEMENTATION

For a quickstart, please see https://ribosomeprofiling.github.io. Source code, installation instructions and links to documentation are available on GitHub: https://github.com/ribosomeprofiling.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification, evolution and expression analyses of whole genome-wide TLP gene family in Brassica napus

BACKGROUND

Brassica is a very important genus of Brassicaceae, including many important oils, vegetables, forage crops, and ornamental horticultural plants. TLP family genes play important regulatory roles in the growth and development of plants. Therefore, this study used a bioinformatics approach to conduct the systematic comparative genomics analysis of TLP gene family in B. napus and other three important Brassicaceae crops.

RESULTS

Here, we identified a total of 29 TLP genes from B. napus genome, and they distributed on 16 chromosomes of B. napus. The evolutionary relationship showed that these genes could be divided into six groups from Group A to F. We found that the gene corresponding to Arabidopsis thaliana AT1G43640 was completely lost in B. rapa, B. oleracea and B. napus after whole genome triplication. The gene corresponding to AT1G25280 was retained in all the three species we analysed, belonging to 1:3:6 ratios. Our analyses suggested that there was a selective loss of some genes that might be redundant after genome duplication. This study proposed that the TLP genes in B. napus did not directly expansion compared with its diploid parents B. rapa, and B. oleracea. Instead, an indirect expansion of TLP gene family occurred in its two diploid parents. In addition, the study further utilized RNA-seq to detect the expression pattern of TLP genes between different tissues and two subgenomes.

CONCLUSIONS

This study systematically conducted the comparative analyses of TLP gene family in B. napus, discussed the loss and expansion of genes after genome duplication. It provided rich gene resources for exploring the molecular mechanism of TLP gene family. Meanwhile, it provided guidance and reference for the research of other gene families in B. napus.

RiboVIEW: a computational framework for visualization, quality control and statistical analysis of ribosome profiling data

Recently, newly developed ribosome profiling methods based on high-throughput sequencing of ribosome-protected mRNA footprints allow to study genome-wide translational changes in detail. However, computational analysis of the sequencing data still represents a bottleneck for many laboratories. Further, specific pipelines for quality control and statistical analysis of ribosome profiling data, providing high levels of both accuracy and confidence, are currently lacking. In this study, we describe automated bioinformatic and statistical diagnoses to perform robust quality control of ribosome profiling data (RiboQC), to efficiently visualize ribosome positions and to estimate ribosome speed (RiboMine) in an unbiased way. We present an R pipeline to setup and undertake the analyses that offers the user an HTML page to scan own data regarding the following aspects: periodicity, ligation and digestion of footprints; reproducibility and batch effects of replicates; drug-related artifacts; unbiased codon enrichment including variability between mRNAs, for A, P and E sites; mining of some causal or confounding factors. We expect our pipeline to allow an optimal use of the wealth of information provided by ribosome profiling experiments.

An ortholog of the Vasa intronic gene is required for small RNA-mediated translation repression in Chlamydomonas reinhardtii

Small RNAs (sRNAs) are a class of noncoding RNAs that regulate complementary mRNAs, by triggering translation repression and/or transcript decay, and influence multiple biological processes. In animals, land plants, and some protists like the alga Chlamydomonas, sRNAs can repress translation of polyribosome-associated mRNAs, without or with only minimal transcript destabilization. However, the precise silencing mechanism is poorly understood. We found that Chlamydomonas VIG1, a homolog of the Drosophila melanogaster Vasa intronic gene and a member of a widely conserved protein family in eukaryotes, is involved in this process. VIG1 appears to be an ancillary ribosomal constituent. Additionally, VIG1 copurifies with core components of sRNA effector complexes and plays a key role in the sRNA-mediated translation repression of polyribosomal transcripts.

Small RNAs (sRNAs) associate with Argonaute (AGO) proteins in effector complexes, termed RNA-induced silencing complexes (RISCs), which regulate complementary transcripts by translation inhibition and/or RNA degradation. In the unicellular alga Chlamydomonas, several metazoans, and land plants, emerging evidence indicates that polyribosome-associated transcripts can be translationally repressed by RISCs without substantial messenger RNA (mRNA) destabilization. However, the mechanism of translation inhibition in a polyribosomal context is not understood. Here we show that Chlamydomonas VIG1, an ortholog of the Drosophila melanogaster Vasa intronic gene (VIG), is required for this process. VIG1 localizes predominantly in the cytosol and comigrates with monoribosomes and polyribosomes by sucrose density gradient sedimentation. A VIG1-deleted mutant shows hypersensitivity to the translation elongation inhibitor cycloheximide, suggesting that VIG1 may have a nonessential role in ribosome function/structure. Additionally, FLAG-tagged VIG1 copurifies with AGO3 and Dicer-like 3 (DCL3), consistent with it also being a component of the RISC. Indeed, VIG1 is necessary for the repression of sRNA-targeted transcripts at the translational level but is dispensable for cleavage-mediated RNA interference and for the association of the AGO3 effector with polyribosomes or target transcripts. Our results suggest that VIG1 is an ancillary ribosomal component and plays a role in sRNA-mediated translation repression of polyribosomal transcripts.

Genome-wide translational reprogramming of genes important for myocyte functions in overload-induced heart failure

Genome-wide changes in gene translational efficiency during the development of heart failure are poorly understood. We tested the hypothesis that aberrant changes in translational efficiency of cardiac genes are associated with the development of myocyte decompensation in response to persistent stress stimuli. We demonstrated that chronic pressure overload in mice resulted in a genome-wide reprogramming of translational efficiency, with >50% of the translatome exhibiting decreased translational efficiencies during the transition from myocardial compensation to decompensation. Importantly, these translationally repressed genes included those involved in angiogenesis and energy metabolism. Moreover, we showed that the stress-induced translational reprogramming was accompanied by persistent activation of the eukaryotic initiation factor 2α (eIF2α)-mediated stress response pathway. Counteracting the endogenous eIF2α functions by cardiac-specific overexpression of an eIF2α-S51A mutant ameliorated the development of myocyte decompensation, with concomitant improvements in translation of cardiac functional genes and increases in angiogenic responses. These data suggest that the mismatch between transcription and translation of the cardiac genes with essential functions may represent a novel molecular mechanism underlying the development of myocyte decompensation in response to chronic stress stimuli, and the eIF2α pathway may be a viable therapeutic target for recovering the optimal translation of the repressed cardiac genes.