Ribosome Profiling as a Tool to Decipher Viral Complexity

Ribosome profiling reveals downregulation of UMP biosynthesis as the major early response to phage infection

Bacteria have evolved diverse defense mechanisms to counter bacteriophage attacks. Genetic programs activated upon infection characterize phage-host molecular interactions and ultimately determine the outcome of the infection. In this study, we applied ribosome profiling to monitor protein synthesis during the early stages of sk1 bacteriophage infection in Lactococcus cremoris. Our analysis revealed major changes in gene expression within 5 minutes of sk1 infection. Notably, we observed a specific and severe downregulation of several pyr operons which encode enzymes required for uridine monophosphate biosynthesis. Consistent with previous findings, this is likely an attempt of the host to starve the phage of nucleotides it requires for propagation. We also observed a gene expression response that we expect to benefit the phage. This included the upregulation of 40 ribosome proteins that likely increased the host's translational capacity, concurrent with a downregulation of genes that promote translational fidelity (lepA and raiA). In addition to the characterization of host-phage gene expression responses, the obtained ribosome profiling data enabled us to identify two putative recoding events as well as dozens of loci currently annotated as pseudogenes that are actively translated. Furthermore, our study elucidated alterations in the dynamics of the translation process, as indicated by time-dependent changes in the metagene profile, suggesting global shifts in translation rates upon infection. Additionally, we observed consistent modifications in the ribosome profiles of individual genes, which were apparent as early as 2 minutes post-infection. The study emphasizes our ability to capture rapid alterations of gene expression during phage infection through ribosome profiling.

IMPORTANCE

The ribosome profiling technology has provided invaluable insights for understanding cellular translation and eukaryotic viral infections. However, its potential for investigating host-phage interactions remains largely untapped. Here, we applied ribosome profiling to Lactococcus cremoris cultures infected with sk1, a major infectious agent in dairy fermentation processes. This revealed a profound downregulation of genes involved in pyrimidine nucleotide synthesis at an early stage of phage infection, suggesting an anti-phage program aimed at restricting nucleotide availability and, consequently, phage propagation. This is consistent with recent findings and contributes to our growing appreciation for the role of nucleotide limitation as an anti-viral strategy. In addition to capturing rapid alterations in gene expression levels, we identified translation occurring outside annotated regions, as well as signatures of non-standard translation mechanisms. The gene profiles revealed specific changes in ribosomal densities upon infection, reflecting alterations in the dynamics of the translation process.

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.

A hidden translatome in tumors-the coding lncRNAs

Long noncoding RNAs (lncRNAs) have been extensively identified in eukaryotic genomes and have been shown to play critical roles in the development of multiple cancers. Through the application and development of ribosome analysis and sequencing technologies, advanced studies have discovered the translation of lncRNAs. Although lncRNAs were originally defined as noncoding RNAs, many lncRNAs actually contain small open reading frames that are translated into peptides. This opens a broad area for the functional investigation of lncRNAs. Here, we introduce prospective methods and databases for screening lncRNAs with functional polypeptides. We also summarize the specific lncRNA-encoded proteins and their molecular mechanisms that promote or inhibit cancerous. Importantly, the role of lncRNA-encoded peptides/proteins holds promise in cancer research, but some potential challenges remain unresolved. This review includes reports on lncRNA-encoded peptides or proteins in cancer, aiming to provide theoretical basis and related references to facilitate the discovery of more functional peptides encoded by lncRNA, and to further develop new anti-cancer therapeutic targets as well as clinical biomarkers of diagnosis and prognosis.

Translation initiation landscape profiling reveals hidden open-reading frames required for the pathogenesis of tomato yellow leaf curl Thailand virus

Plant viruses with densely packed genomes employ noncanonical translational strategies to increase the coding capacity for viral function. However, the diverse translational strategies used make it challenging to define the full set of viral genes. Here, using tomato yellow leaf curl Thailand virus (TYLCTHV, genus Begomovirus) as a model system, we identified genes beyond the annotated gene sets by experimentally profiling in vivo translation initiation sites (TISs). We found that unanticipated AUG TISs were prevalent and determined that their usage involves alternative transcriptional and/or translational start sites and is associated with flanking mRNA sequences. Specifically, two downstream in-frame TISs were identified in the viral gene AV2. These TISs were conserved in the begomovirus lineage and led to the translation of different protein isoforms localized to cytoplasmic puncta and at the cell periphery, respectively. In addition, we found translational evidence of an unexplored gene, BV2. BV2 is conserved among TYLCTHV isolates and localizes to the endoplasmic reticulum and plasmodesmata. Mutations of AV2 isoforms and BV2 significantly attenuated disease symptoms in tomato (Solanum lycopersicum). In conclusion, our study pinpointing in vivo TISs untangles the coding complexity of a plant viral genome and, more importantly, illustrates the biological significance of the hidden open-reading frames encoding viral factors for pathogenicity.

The complex biology of human cytomegalovirus latency

While many viral infections are limited and eventually resolved by the host immune response or by death of the host, other viruses establish long-term relationships with the host by way of a persistent infection, that range from chronic viruses that may be eventually cleared to those that establish life-long persistent or latent infection. Viruses infecting hosts from bacteria to humans establish quiescent infections that must be reactivated to produce progeny. For mammalian viruses, most notably herpesviruses, this quiescent maintenance of viral genomes in the absence of virus replication is referred to as latency. The latent strategy allows the virus to persist quiescently within a single host until conditions indicate a need to reactivate to reach a new host or, to re-seed a reservoir within the host. Here, I review common themes in viral strategies to regulate the latent cycle and reactivate from it ranging from bacteriophage to herpesviruses with a focus on human cytomegalovirus (HCMV). Themes central to herpesvirus latency include, epigenetic repression of viral gene expression and mechanisms to regulate host signaling and survival. Critical to the success of a latent program are mechanisms by which the virus can "sense" fluctuations in host biology (within the host) or environment (outside the host) and make appropriate "decisions" to maintain latency or re-initiate the replicative program. The signals or environments that indicate the establishment of a latent state, the very nature of the latent state, as well as the signals driving reactivation have been topics of intense study from bacteriophage to human viruses, as these questions encompass the height of complexity in virus-host interactions-where the host and the virus coexist.

RiboChat: a chat-style web interface for analysis and annotation of ribosome profiling data

The increasing volume of ribosome profiling (Ribo-seq) data, computational complexity of its data processing and operational handicap of related analytical procedures present a daunting set of informatics challenges. These impose a substantial barrier to researchers particularly with no or limited bioinformatics expertise in analyzing and decoding translation information from Ribo-seq data, thus driving the need for a new research paradigm for data computation and information extraction. In this knowledge base, we herein present a novel interactive web platform, RiboChat (https://db.cngb.org/ribobench/chat.html), for direct analyzing and annotating Ribo-seq data in the form of a chat conversation. It consists of a user-friendly web interface and a backend cloud-computing service. When typing a data analysis question into the chat window, the object-text detection module will be run to recognize relevant keywords from the input text. Based on the features identified in the input, individual analytics modules are then scored to find the perfect-matching candidate. The corresponding analytics module will be further executed after checking the completion status of the uploading of datasets and configured parameters. Overall, RiboChat represents an important step forward in the emerging direction of next-generation data analytics and will enable the broad research community to conveniently decipher translation information embedded within Ribo-seq data.

Estimation of peptide elongation times from ribosome profiling spectra

Ribosome profiling spectra bear rich information on translation control and dynamics. Yet, due to technical biases in library generation, extracting quantitative measures of discrete translation events has remained elusive. Using maximum likelihood statistics and data set from Escherichia coli we develop a robust method for neutralizing technical biases (e.g. base specific RNase preferences in ribosome-protected mRNA fragments (RPF) generation), which allows for correct estimation of translation times at single codon resolution. Furthermore, we validated the method with available datasets from E. coli treated with antibiotic to inhibit isoleucyl-tRNA synthetase, and two datasets from Saccharomyces cerevisiae treated with two RNases with distinct cleavage signatures. We demonstrate that our approach accounts for RNase cleavage preferences and provides bias-corrected translation times estimates. Our approach provides a solution to the long-standing problem of extracting reliable information about peptide elongation times from highly noisy and technically biased ribosome profiling spectra.

The effects of codon bias and optimality on mRNA and protein regulation

The central dogma of molecular biology entails that genetic information is transferred from nucleic acid to proteins. Notwithstanding retro-transcribing genetic elements, DNA is transcribed to RNA which in turn is translated into proteins. Recent advancements have shown that each stage is regulated to control protein abundances for a variety of essential physiological processes. In this regard, mRNA regulation is essential in fine-tuning or calibrating protein abundances. In this review, we would like to discuss one of several mRNA-intrinsic features of mRNA regulation that has been gaining traction of recent-codon bias and optimality. Specifically, we address the effects of codon bias with regard to codon optimality in several biological processes centred on translation, such as mRNA stability and protein folding among others. Finally, we examine how different organisms or cell types, through this system, are able to coordinate physiological pathways to respond to a variety of stress or growth conditions.

Identification of a herpes simplex virus 1 gene encoding neurovirulence factor by chemical proteomics

Identification of the complete set of translated genes of viruses is important to understand viral replication and pathogenesis as well as for therapeutic approaches to control viral infection. Here, we use chemical proteomics, integrating bio-orthogonal non-canonical amino acid tagging and high-resolution mass spectrometry, to characterize the newly synthesized herpes simplex virus 1 (HSV-1) proteome in infected cells. In these infected cells, host cellular protein synthesis is shut-off, increasing the chance to preferentially detect viral proteomes. We identify nine previously cryptic orphan protein coding sequences whose translated products are expressed in HSV-1-infected cells. Functional characterization of one identified protein, designated piUL49, shows that it is critical for HSV-1 neurovirulence in vivo by regulating the activity of virally encoded dUTPase, a key enzyme that maintains accurate DNA replication. Our results demonstrate that cryptic orphan protein coding genes of HSV-1, and probably other large DNA viruses, remain to be identified.

Here the authors use chemical proteomics to identify the herpes simplex virus 1 encoded proteome in infected cells. Functional characterization of one of the nine identified proteins, designated piUL49, shows that it acts as neurovirulence factor in mice by regulating a virally encoded dUTPase.

Translation of small downstream ORFs enhances translation of canonical main open reading frames

In addition to canonical open reading frames (ORFs), thousands of translated small ORFs (containing less than 100 codons) have been identified in untranslated mRNA regions (UTRs) across eukaryotes. Small ORFs in 5′ UTRs (upstream (u)ORFs) often repress translation of the canonical ORF within the same mRNA. However, the function of translated small ORFs in the 3′ UTRs (downstream (d)ORFs) is unknown. Contrary to uORFs, we find that translation of dORFs enhances translation of their corresponding canonical ORFs. This translation stimulatory effect of dORFs depends on the number of dORFs, but not the length or peptide they encode. We propose that dORFs represent a new, strong, and universal translation regulatory mechanism in vertebrates.

Trips-Viz: a transcriptome browser for exploring Ribo-Seq data

Ribosome profiling (Ribo-Seq) is a technique that allows for the isolation and sequencing of mRNA fragments protected from nuclease digestion by actively translating ribosomes. Mapping these ribosome footprints to a genome or transcriptome generates quantitative information on translated regions. To provide access to publicly available ribosome profiling data in the context of transcriptomes we developed Trips-Viz (transcriptome-wide information on protein synthesis-visualized). Trips-Viz provides a large range of graphical tools for exploring global properties of translatomes and of individual transcripts. It enables analysis of aligned footprints to evaluate datasets quality, differential gene expression detection, visual identification of upstream ORFs and alternative proteoforms. Trips-Viz is available at https://trips.ucc.ie.

A simplified qPCR method revealing tRNAome remodeling upon infection by genotype 3 hepatitis E virus

The landscape of tRNA-viral codons regulates viral adaption at the translational level, presumably through adapting to host codon usage or modulating the host tRNAome. We found that the major zoonotic genotype of hepatitis E virus (HEV) has not adapted to host codon usage, prompting exploration of the effects of HEV infection on the host tRNAome. However, tRNAome quantification is largely impeded by the extremely short sequences of tRNAs and redundancy of tRNA genes. Here, we present a length-extension and stepwise simplified qPCR method that utilizes a universal DNA/RNA hybrid tRNA adaptor and degenerate primers. Using this novel methodology, we observe that HEV infection dramatically reprograms the hepatic tRNAome, which is likely to facilitate translation of viral RNAs. This tRNAome quantification method bears broad implications for future tRNA research and possibly tRNA-based diagnostics.

GWIPS-viz: 2018 update

The GWIPS-viz browser (http://gwips.ucc.ie/) is an on-line genome browser which is tailored for exploring ribosome profiling (Ribo-seq) data. Since its publication in 2014, GWIPS-viz provides Ribo-seq data for an additional 14 genomes bringing the current total to 23. The integration of new Ribo-seq data has been automated thereby increasing the number of available tracks to 1792, a 10-fold increase in the last three years. The increase is particularly substantial for data derived from human sources. Following user requests, we added the functionality to download these tracks in bigWig format. We also incorporated new types of data (e.g. TCP-seq) as well as auxiliary tracks from other sources that help with the interpretation of Ribo-seq data. Improvements in the visualization of the data have been carried out particularly for bacterial genomes where the Ribo-seq data are now shown in a strand specific manner. For higher eukaryotic datasets, we provide characteristics of individual datasets using the RUST program which includes the triplet periodicity, sequencing biases and relative inferred A-site dwell times. This information can be used for assessing the quality of Ribo-seq datasets. To improve the power of the signal, we aggregate Ribo-seq data from several studies into Global aggregate tracks for each genome.

Overlapping genes and the proteins they encode differ significantly in their sequence composition from non-overlapping genes

Overlapping genes represent a fascinating evolutionary puzzle, since they encode two functionally unrelated proteins from the same DNA sequence. They originate by a mechanism of overprinting, in which point mutations in an existing frame allow the expression (the "birth") of a completely new protein from a second frame. In viruses, in which overlapping genes are abundant, these new proteins often play a critical role in infection, yet they are frequently overlooked during genome annotation. This results in erroneous interpretation of mutational studies and in a significant waste of resources. Therefore, overlapping genes need to be correctly detected, especially since they are now thought to be abundant also in eukaryotes. Developing better detection methods and conducting systematic evolutionary studies require a large, reliable benchmark dataset of known cases. We thus assembled a high-quality dataset of 80 viral overlapping genes whose expression is experimentally proven. Many of them were not present in databases. We found that overall, overlapping genes differ significantly from non-overlapping genes in their nucleotide and amino acid composition. In particular, the proteins they encode are enriched in high-degeneracy amino acids and depleted in low-degeneracy ones, which may alleviate the evolutionary constraints acting on overlapping genes. Principal component analysis revealed that the vast majority of overlapping genes follow a similar composition bias, despite their heterogeneity in length and function. Six proven mammalian overlapping genes also followed this bias. We propose that this apparently near-universal composition bias may either favour the birth of overlapping genes, or/and result from selection pressure acting on them.

Following Ribosome Footprints to Understand Translation at a Genome Wide Level

Protein translation is a key step in gene expression. The development of Ribosome Profiling has allowed the global analysis of this process at sub-codon resolution. In the last years the method has been applied to several models ranging from bacteria to mammalian cells yielding a surprising amount of insight on the mechanism and the regulation of translation. In this review we describe the key aspects of the experimental protocol and comment on the main conclusions raised in different models.

Evidence of translation efficiency adaptation of the coding regions of the bacteriophage lambda

Deciphering the way gene expression regulatory aspects are encoded in viral genomes is a challenging mission with ramifications related to all biomedical disciplines. Here, we aimed to understand how the evolution shapes the bacteriophage lambda genes by performing a high resolution analysis of ribosomal profiling data and gene expression related synonymous/silent information encoded in bacteriophage coding regions.

We demonstrated evidence of selection for distinct compositions of synonymous codons in early and late viral genes related to the adaptation of translation efficiency to different bacteriophage developmental stages. Specifically, we showed that evolution of viral coding regions is driven, among others, by selection for codons with higher decoding rates; during the initial/progressive stages of infection the decoding rates in early/late genes were found to be superior to those in late/early genes, respectively. Moreover, we argued that selection for translation efficiency could be partially explained by adaptation to Escherichia coli tRNA pool and the fact that it can change during the bacteriophage life cycle.

An analysis of additional aspects related to the expression of viral genes, such as mRNA folding and more complex/longer regulatory signals in the coding regions, is also reported. The reported conclusions are likely to be relevant also to additional viruses.

Proteogenomics in Aid of Host-Pathogen Interaction Studies: A Bacterial Perspective

By providing useful tools to study host-pathogen interactions, next-generation omics has recently enabled the study of gene expression changes in both pathogen and infected host simultaneously. However, since great discriminative power is required to study pathogen and host simultaneously throughout the infection process, the depth of quantitative gene expression profiling has proven to be unsatisfactory when focusing on bacterial pathogens, thus preferentially requiring specific strategies or the development of novel methodologies based on complementary omics approaches. In this review, we focus on the difficulties encountered when making use of proteogenomics approaches to study bacterial pathogenesis. In addition, we review different omics strategies (i.e., transcriptomics, proteomics and secretomics) and their applications for studying interactions of pathogens with their host.

Optimal Down Regulation of mRNA Translation

Down regulation of mRNA translation is an important problem in various bio-medical domains ranging from developing effective medicines for tumors and for viral diseases to developing attenuated virus strains that can be used for vaccination. Here, we study the problem of down regulation of mRNA translation using a mathematical model called the ribosome flow model (RFM). In the RFM, the mRNA molecule is modeled as a chain of n sites. The flow of ribosomes between consecutive sites is regulated by n + 1 transition rates. Given a set of feasible transition rates, that models the outcome of all possible mutations, we consider the problem of maximally down regulating protein production by altering the rates within this set of feasible rates. Under certain conditions on the feasible set, we show that an optimal solution can be determined efficiently. We also rigorously analyze two special cases of the down regulation optimization problem. Our results suggest that one must focus on the position along the mRNA molecule where the transition rate has the strongest effect on the protein production rate. However, this rate is not necessarily the slowest transition rate along the mRNA molecule. We discuss some of the biological implications of these results.

Long and Short Isoforms of the Human Cytomegalovirus UL138 Protein Silence IE Transcription and Promote Latency

The UL133-138 locus present in clinical strains of human cytomegalovirus (HCMV) encodes proteins required for latency and reactivation in CD34(+) hematopoietic progenitor cells and virion maturation in endothelial cells. The encoded proteins form multiple homo- and hetero-interactions and localize within secretory membranes. One of these genes, UL136 gene, is expressed as at least five different protein isoforms with overlapping and unique functions. Here we show that another gene from this locus, the UL138 gene, also generates more than one protein isoform. A long form of UL138 (pUL138-L) initiates translation from codon 1, possesses an amino-terminal signal sequence, and is a type one integral membrane protein. Here we identify a short protein isoform (pUL138-S) initiating from codon 16 that displays a subcellular localization similar to that of pUL138-L. Reporter, short-term transcription, and long-term virus production assays revealed that both pUL138-L and pUL138-S are able to suppress major immediate early (IE) gene transcription and the generation of infectious virions in cells in which HCMV latency is studied. The long form appears to be more potent at silencing IE transcription shortly after infection, while the short form seems more potent at restricting progeny virion production at later times, indicating that both isoforms of UL138 likely cooperate to promote HCMV latency.

IMPORTANCE

Latency allows herpesviruses to persist for the lives of their hosts in the face of effective immune control measures for productively infected cells. Controlling latent reservoirs is an attractive antiviral approach complicated by knowledge deficits for how latently infected cells are established, maintained, and reactivated. This is especially true for betaherpesviruses. The functional consequences of HCMV UL138 protein expression during latency include repression of viral IE1 transcription and suppression of virus replication. Here we show that short and long isoforms of UL138 exist and can themselves support latency but may do so in temporally distinct manners. Understanding the complexity of gene expression and its impact on latency is important for considering potential antivirals targeting latent reservoirs.