PausePred and Rfeet: webtools for inferring ribosome pauses and visualizing footprint density from ribosome profiling data

Dissecting infectious bronchitis virus-induced host shutoff at the translation level

Viruses have evolved a range of strategies to utilize or manipulate the host's cellular translational machinery for efficient infection, although the mechanisms by which infectious bronchitis virus (IBV) manipulates the host translation machinery remain unclear. In this study, we firstly demonstrate that IBV infection causes host shutoff, although viral protein synthesis is not affected. We then screened 23 viral proteins, and identified that more than one viral protein is responsible for IBV-induced host shutoff, the inhibitory effects of proteins Nsp15 were particularly pronounced. Ribosome profiling was used to draw the landscape of viral mRNA and cellular genes expression model, and the results showed that IBV mRNAs gradually dominated the cellular mRNA pool, the translation efficiency of the viral mRNAs was lower than the median efficiency (about 1) of cellular mRNAs. In the analysis of viral transcription and translation, higher densities of RNA sequencing (RNA-seq) and ribosome profiling (Ribo-seq) reads were observed for structural proteins and 5' untranslated regions, which conformed to the typical transcriptional characteristics of nested viruses. Translational halt events and the number of host genes increased significantly after viral infection. The translationally paused genes were enriched in translation, unfolded-protein-related response, and activation of immune response pathways. Immune- and inflammation-related mRNAs were inefficiently translated in infected cells, and IBV infection delayed the production of IFN-β and IFN-λ. Our results describe the translational landscape of IBV-infected cells and demonstrate new strategies by which IBV induces host gene shutoff to promote its replication.

IMPORTANCE

Infectious bronchitis virus (IBV) is a γ-coronavirus that causes huge economic losses to the poultry industry. Understanding how the virus manipulates cellular biological processes to facilitate its replication is critical for controlling viral infections. Here, we used Ribo-seq to determine how IBV infection remodels the host's biological processes and identified multiple viral proteins involved in host gene shutoff. Immune- and inflammation-related mRNAs were inefficiently translated, the translation halt of unfolded proteins and immune activation-related genes increased significantly, benefitting IBV replication. These data provide new insights into how IBV modulates its host's antiviral responses.

Polyglutamine-mediated ribotoxicity disrupts proteostasis and stress responses in Huntington's disease

Huntington's disease (HD) is a neurodegenerative disorder caused by expansion of a CAG trinucleotide repeat in the Huntingtin (HTT) gene, encoding a homopolymeric polyglutamine (polyQ) tract. Although mutant HTT (mHTT) protein is known to aggregate, the links between aggregation and neurotoxicity remain unclear. Here we show that both translation and aggregation of wild-type HTT and mHTT are regulated by a stress-responsive upstream open reading frame and that polyQ expansions cause abortive translation termination and release of truncated, aggregation-prone mHTT fragments. Notably, we find that mHTT depletes translation elongation factor eIF5A in brains of symptomatic HD mice and cultured HD cells, leading to pervasive ribosome pausing and collisions. Loss of eIF5A disrupts homeostatic controls and impairs recovery from acute stress. Importantly, drugs that inhibit translation initiation reduce premature termination and mitigate this escalating cascade of ribotoxic stress and dysfunction in HD.

The ribosome quality control factor Asc1 determines the fate of HSP70 mRNA on and off the ribosome

Cells survive harsh environmental conditions by potently upregulating molecular chaperones such as heat shock proteins (HSPs), particularly the inducible members of the HSP70 family. The life cycle of HSP70 mRNA in the cytoplasm is unique-it is translated during stress when most cellular mRNA translation is repressed and rapidly degraded upon recovery. Contrary to its 5' untranslated region's role in maximizing translation, we discovered that the HSP70 coding sequence (CDS) suppresses its translation via the ribosome quality control (RQC) mechanism. The CDS of the most inducible Saccharomyces cerevisiae HSP70 gene, SSA4, is uniquely enriched with low-frequency codons that promote ribosome stalling during heat stress. Stalled ribosomes are recognized by the RQC components Asc1p and Hel2p and two novel RQC components, the ribosomal proteins Rps28Ap and Rps19Bp. Surprisingly, RQC does not signal SSA4 mRNA degradation via No-Go-Decay. Instead, Asc1p destabilizes SSA4 mRNA during recovery from heat stress by a mechanism independent of ribosome binding and SSA4 codon optimality. Therefore, Asc1p operates in two pathways that converge to regulate the SSA4 mRNA life cycle during stress and recovery. Our research identifies Asc1p as a critical regulator of the stress response and RQC as the mechanism tuning HSP70 synthesis.

Comparative oxidation proteomics analyses suggest redox regulation of cytosolic translation in rice leaves upon Magnaporthe oryzae infection

Pathogen attack can increase plant levels of reactive oxygen species (ROS), which act as signaling molecules to activate plant defense mechanisms. Elucidating these processes is crucial for understanding redox signaling pathways in plant defense responses. Using an iodo-tandem mass tag (TMT)-based quantitative proteomics approach, we mapped 3362 oxidized cysteine sites in 2275 proteins in rice leaves. Oxidized proteins were involved in gene expression, peptide biosynthetic processes, stress responses, ROS metabolic processes, and translation pathways. Magnaporthe oryzae infection led to increased oxidative modification levels of 512 cysteine sites in 438 proteins, including many transcriptional regulators and ribosomal proteins. Ribosome profiling (Ribo-seq) analysis revealed that the oxidative modification of ribosomal proteins promoted the translational efficiency of many mRNAs involved in defense response pathways, thereby affecting rice immunity. Our results suggest that increased oxidative modification of ribosomal proteins in rice leaves promotes cytosolic translation, thus revealing a novel function of post-translational modifications. Furthermore, the oxidation-sensitive proteins identified here provide a valuable resource for research on protein redox regulation and can guide future mechanistic studies.

Velcrin-induced selective cleavage of tRNALeu(TAA) by SLFN12 causes cancer cell death

Velcrin compounds kill cancer cells expressing high levels of phosphodiesterase 3A (PDE3A) and Schlafen family member 12 (SLFN12) by inducing complex formation between these two proteins, but the mechanism of cancer cell killing by the PDE3A-SLFN12 complex is not fully understood. Here, we report that the physiological substrate of SLFN12 RNase is tRNALeu(TAA). SLFN12 selectively digests tRNALeu(TAA), and velcrin treatment promotes the cleavage of tRNALeu(TAA) by inducing PDE3A-SLFN12 complex formation in vitro. We found that distinct sequences in the variable loop and acceptor stem of tRNALeu(TAA) are required for substrate digestion. Velcrin treatment of sensitive cells results in downregulation of tRNALeu(TAA), ribosome pausing at Leu-TTA codons and global inhibition of protein synthesis. Velcrin-induced cleavage of tRNALeu(TAA) by SLFN12 and the concomitant global inhibition of protein synthesis thus define a new mechanism of apoptosis initiation.

METTL1-mediated m7G modification of Arg-TCT tRNA drives oncogenic transformation

The emerging "epitranscriptomics" field is providing insights into the biological and pathological roles of different RNA modifications. The RNA methyltransferase METTL1 catalyzes N7-methylguanosine (m7G) modification of tRNAs. Here we find METTL1 is frequently amplified and overexpressed in cancers and is associated with poor patient survival. METTL1 depletion causes decreased abundance of m7G-modified tRNAs and altered cell cycle and inhibits oncogenicity. Conversely, METTL1 overexpression induces oncogenic cell transformation and cancer. Mechanistically, we find increased abundance of m7G-modified tRNAs, in particular Arg-TCT-4-1, and increased translation of mRNAs, including cell cycle regulators that are enriched in the corresponding AGA codon. Accordingly, Arg-TCT expression is elevated in many tumor types and is associated with patient survival, and strikingly, overexpression of this individual tRNA induces oncogenic transformation. Thus, METTL1-mediated tRNA modification drives oncogenic transformation through a remodeling of the mRNA "translatome" to increase expression of growth-promoting proteins and represents a promising anti-cancer target.

Histone deacetylases control lysine acetylation of ribosomal proteins in rice

Lysine acetylation (Kac) is well known to occur in histones for chromatin function and epigenetic regulation. In addition to histones, Kac is also detected in a large number of proteins with diverse biological functions. However, Kac function and regulatory mechanism for most proteins are unclear. In this work, we studied mutation effects of rice genes encoding cytoplasm-localized histone deacetylases (HDAC) on protein acetylome and found that the HDAC protein HDA714 was a major deacetylase of the rice non-histone proteins including many ribosomal proteins (r-proteins) and translation factors that were extensively acetylated. HDA714 loss-of-function mutations increased Kac levels but reduced abundance of r-proteins. In vitro and in vivo experiments showed that HDA714 interacted with r-proteins and reduced their Kac. Substitutions of lysine by arginine (depleting Kac) in several r-proteins enhance, while mutations of lysine to glutamine (mimicking Kac) decrease their stability in transient expression system. Ribo-seq analysis revealed that the hda714 mutations resulted in increased ribosome stalling frequency. Collectively, the results uncover Kac as a functional posttranslational modification of r-proteins which is controlled by histone deacetylases, extending the role of Kac in gene expression to protein translational regulation.

Mutant Huntingtin stalls ribosomes and represses protein synthesis in a cellular model of Huntington disease

The polyglutamine expansion of huntingtin (mHTT) causes Huntington disease (HD) and neurodegeneration, but the mechanisms remain unclear. Here, we found that mHtt promotes ribosome stalling and suppresses protein synthesis in mouse HD striatal neuronal cells. Depletion of mHtt enhances protein synthesis and increases the speed of ribosomal translocation, while mHtt directly inhibits protein synthesis in vitro. Fmrp, a known regulator of ribosome stalling, is upregulated in HD, but its depletion has no discernible effect on protein synthesis or ribosome stalling in HD cells. We found interactions of ribosomal proteins and translating ribosomes with mHtt. High-resolution global ribosome footprint profiling (Ribo-Seq) and mRNA-Seq indicates a widespread shift in ribosome occupancy toward the 5' and 3' end and unique single-codon pauses on selected mRNA targets in HD cells, compared to controls. Thus, mHtt impedes ribosomal translocation during translation elongation, a mechanistic defect that can be exploited for HD therapeutics.

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.

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.

Cyclic GMP-AMP synthase promotes the inflammatory and autophagy responses in Huntington disease

Huntington disease (HD) is a genetic disorder caused by glutamine-expansion in the huntingtin (mHTT) protein, which affects motor, psychiatric, and cognitive function, but the mechanisms remain unclear. mHTT is known to induce DNA damage and affect autophagy, both associated with inflammatory responses, but what mediates all these were unknown. Here we report that cGAS, a DNA damage sensor, is highly upregulated in the striatum of a mouse model and HD human patient’s tissue. We found ribosomes, which make proteins, are robustly accumulated on the cGAS mRNA in HD cells. cGAS depletion decreases—and cGAS expression increases—both inflammatory and autophagy responses in HD striatal cells. Thus, cGAS is a therapeutic target for HD. Blocking cGAS will prevent/slow down HD symptoms.

Huntington disease (HD) is caused by an expansion mutation of the N-terminal polyglutamine of huntingtin (mHTT). mHTT is ubiquitously present, but it induces noticeable damage to the brain’s striatum, thereby affecting motor, psychiatric, and cognitive functions. The striatal damage and progression of HD are associated with the inflammatory response; however, the underlying molecular mechanisms remain unclear. Here, we report that cGMP-AMP synthase (cGAS), a DNA sensor, is a critical regulator of inflammatory and autophagy responses in HD. Ribosome profiling revealed that the cGAS mRNA has high ribosome occupancy at exon 1 and codon-specific pauses at positions 171 (CCG) and 172 (CGT) in HD striatal cells. Moreover, the protein levels and activity of cGAS (based on the phosphorylated STING and phosphorylated TBK1 levels), and the expression and ribosome occupancy of cGAS-dependent inflammatory genes (Ccl5 and Cxcl10) are increased in HD striatum. Depletion of cGAS diminishes cGAS activity and decreases the expression of inflammatory genes while suppressing the up-regulation of autophagy in HD cells. In contrast, reinstating cGAS in cGAS-depleted HD cells activates cGAS activity and promotes inflammatory and autophagy responses. Ribosome profiling also revealed that LC3A and LC3B, the two major autophagy initiators, show altered ribosome occupancy in HD cells. We also detected the presence of numerous micronuclei, which are known to induce cGAS, in the cytoplasm of neurons derived from human HD embryonic stem cells. Collectively, our results indicate that cGAS is up-regulated in HD and mediates inflammatory and autophagy responses. Thus, targeting the cGAS pathway may offer therapeutic benefits in HD.

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.

XPRESSyourself: Enhancing, standardizing, and automating ribosome profiling computational analyses yields improved insight into data

Ribosome profiling, an application of nucleic acid sequencing for monitoring ribosome activity, has revolutionized our understanding of protein translation dynamics. This technique has been available for a decade, yet the current state and standardization of publicly available computational tools for these data is bleak. We introduce XPRESSyourself, an analytical toolkit that eliminates barriers and bottlenecks associated with this specialized data type by filling gaps in the computational toolset for both experts and non-experts of ribosome profiling. XPRESSyourself automates and standardizes analysis procedures, decreasing time-to-discovery and increasing reproducibility. This toolkit acts as a reference implementation of current best practices in ribosome profiling analysis. We demonstrate this toolkit’s performance on publicly available ribosome profiling data by rapidly identifying hypothetical mechanisms related to neurodegenerative phenotypes and neuroprotective mechanisms of the small-molecule ISRIB during acute cellular stress. XPRESSyourself brings robust, rapid analysis of ribosome-profiling data to a broad and ever-expanding audience and will lead to more reproducible and accessible measurements of translation regulation. XPRESSyourself software is perpetually open-source under the GPL-3.0 license and is hosted at https://github.com/XPRESSyourself, where users can access additional documentation and report software issues.

Autism-Misregulated eIF4G Microexons Control Synaptic Translation and Higher Order Cognitive Functions

Microexons represent the most highly conserved class of alternative splicing, yet their functions are poorly understood. Here, we focus on closely related neuronal microexons overlapping prion-like domains in the translation initiation factors, eIF4G1 and eIF4G3, the splicing of which is activity dependent and frequently disrupted in autism. CRISPR-Cas9 deletion of these microexons selectively upregulates synaptic proteins that control neuronal activity and plasticity and further triggers a gene expression program mirroring that of activated neurons. Mice lacking the Eif4g1 microexon display social behavior, learning, and memory deficits, accompanied by altered hippocampal synaptic plasticity. We provide evidence that the eIF4G microexons function as a translational brake by causing ribosome stalling, through their propensity to promote the coalescence of cytoplasmic granule components associated with translation repression, including the fragile X mental retardation protein FMRP. The results thus reveal an autism-disrupted mechanism by which alternative splicing specializes neuronal translation to control higher order cognitive functioning.

Comparative Analysis of Gene Expression in Virulent and Attenuated Strains of Infectious Bronchitis Virus at Subcodon Resolution

Like all coronaviruses, avian infectious bronchitis virus (IBV) possesses a long, single-stranded, positive-sense RNA genome (∼27 kb) and has a complex replication strategy that includes the production of a nested set of subgenomic mRNAs (sgmRNAs). Here, we used whole-transcriptome sequencing (RNASeq) and ribosome profiling (RiboSeq) to delineate gene expression in the IBV M41-CK and Beau-R strains at subcodon resolution. RNASeq facilitated a comparative analysis of viral RNA synthesis and revealed two novel transcription junction sites in the attenuated Beau-R strain, one of which would generate a sgmRNA encoding a ribosomally occupied open reading frame (dORF) located downstream of the nucleocapsid coding region. RiboSeq permitted quantification of the translational efficiency of virus gene expression and identified, for the first time, sites of ribosomal pausing on the genome. Quantification of reads flanking the programmed ribosomal frameshifting (PRF) signal at the genomic RNA ORF1a/ORF1b junction revealed that PRF in IBV is highly efficient (33 to 40%). Triplet phasing of RiboSeq data allowed precise determination of reading frames and revealed the translation of two ORFs (ORF4b and ORF4c on sgmRNA IR), which are widely conserved across IBV isolates. Analysis of differential gene expression in infected primary chick kidney cells indicated that the host cell response to IBV occurs primarily at the level of transcription, with global upregulation of immune-related mRNA transcripts following infection and comparatively modest changes in the translation efficiencies of host genes. Cellular genes and gene networks differentially expressed during virus infection were also identified, giving insights into the host cell response to IBV infection.IMPORTANCE IBV is a major avian pathogen and presents a substantial economic burden to the poultry industry. Improved vaccination strategies are urgently needed to curb the global spread of this virus, and the development of suitable vaccine candidates will be aided by an improved understanding of IBV molecular biology. Our high-resolution data have enabled a precise study of transcription and translation in cells infected with both pathogenic and attenuated forms of IBV and expand our understanding of gammacoronaviral gene expression. We demonstrate that gene expression shows considerable intraspecies variation, with single nucleotide polymorphisms being associated with altered production of sgmRNA transcripts, and our RiboSeq data sets enabled us to uncover novel ribosomally occupied ORFs in both strains. The numerous cellular genes and gene networks found to be differentially expressed during virus infection provide insights into the host cell response to IBV infection.