{'fraction_non_alphanumeric': 0.03457577243746935, 'fraction_numerical': 0.012506130456105935, 'mean_word_length': 4.882211538461538, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 9, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
GWIPS-viz: 2018 update

Abstract 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.


INTRODUCTION
Ribosome profiling (Ribo-seq) is a biochemical technique that utilizes high throughput sequencing that captures the mRNA fragments that are protected by actively translating ribosomes (1) thereby providing Genome-Wide Information on Protein Synthesis (GWIPS) (2). Ribo-seq was first carried out in Saccharomyces cerevisiae (1) and has since been used in many organisms resulting in a substantial growth in the number of published datasets. The numerous applications of the ribosome profiling technique as well as its limitations are described in details elsewhere (3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14). While the majority of Ribo-seq datasets represent foot-prints of elongating ribosomes, a number of studies have used protocols for enriching footprints deriving from initiating ribosomes and more recently a modification of the ribosome profiling protocol allowed footprinting of scanning ribosomes (15).
To account for differences in mRNA abundance, most Ribo-seq studies also generate parallel datasets where total mRNA (or total RNA) is randomly degraded and subsequently sequenced. Here we refer to such datasets as mRNA-seq. To date, the majority of published Riboseq/mRNA-seq raw sequencing data have been deposited in NCBI's Sequence Read Archive (SRA) (16).
The GWIPS-viz browser (http://gwips.ucc.ie/) uses the functionality of the UCSC Genome Browser (17) to provide visualizations of Ribo-seq coupled with mRNA-seq controls so that users can freely explore pre-populated Riboseq/mRNA-seq tracks without the need to download, preprocess and align raw sequencing data to the corresponding genomes. Since its original publication (18), we have striven to expand the repertoire of Ribo-seq/mRNA-seq data hosted on GWIPS-viz. We have also incorporated additional tracks as well as improved visualizations to help users better interpret the Ribo-seq/mRNA-seq data.
In addition to individual tracks reflecting Ribo-seq data generated under different conditions for each study, we aggregate each study's data into an All track. We then aggregate the All tracks from each study into a Global Aggregate track for each genome ( Figure 1A-D). This has the effect of improving the overall Ribo-seq signal by reducing the contribution of dataset specific biases and stochastic noise due to low coverage. The increased number of datasets is expected to yield higher sensitivity. The Global Aggregate tracks are set as the default for each genome and users can turn on/off each study's data contribution to the aggregated data and then refine the visualizations by turning on/off individual tracks in each study. In addition, we provide Global Aggregate tracks through the UCSC Genome Browser for the human hg38 and hg19 assemblies.
We have also incorporated a new type of data into GWIPS-viz. Recently the Preiss group developed a technique called translation complex profile sequencing (TCPseq), where ribosome subunits are cross-linked to mRNA, allowing footprinting of both elongating and scanning ribosomes (15). We now provide this data in the Small Ribosomal Subunits (Footprints) track group for S. cerevisiae. Given that TCP-seq is a powerful tool for studying translation initiation (132), we anticipate that the experimental protocol will be adapted for other species quite soon and we will strive to incorporate new TCP-seq data into GWIPS-viz as it becomes available.
For S. cerevisiae we also generated an additional gene annotation track from transcript isoform sequencing (TIF-seq) data (133). Given that Saccharomyces Genome Database (SGD) (134) and Ensembl (135) gene annotations for S. cerevisiae do not include UTR regions, and given that Ribo-seq has shown that extensive translation occurs in 5 leader regions (1) and to a lesser extent in 3 trailer regions (97), we thought it useful to integrate the 5 leader and 3 trailer gene annotations as interpreted by (133) from their TIF-seq data for S. cerevisiae.
With respect to 5 leader region delineation, we incorporated Riken 5 cap analysis gene expression (CAGE) data (136) as a permanent track in the Annotations Tracks and External Data track group for the human hg19 assembly. At the time of the incorporation, Riken CAGE data tracks for other GWIPS-viz assemblies were not available for permanent track integration. However, as GWIPS-viz also now includes the UCSC Genome Browser's Track Hub functionality (137), we provide Riken's FANTOM5 tracks for hg38, hg19, mm9 and rn6 as public track hubs. While these tracks are hosted and managed by the Riken group on their own server, a simple connection makes it easy to explore their CAGE data in conjunction with our Ribo-seq/mRNA-seq tracks in GWIPS-viz.
Initially we did not provide UCSC Genome Browser's custom track feature (138) in GWIPS-viz. The custom track is only accessible to the user who uploads it, i.e. it is not a publicly available track. Many GWIPS-viz users, however, expressed an interest in the custom track feature as a means to explore their own Ribo-seq data in the context of published data and so we now include it. The custom track feature is also particularly useful for users of RiboGalaxy (139), a Galaxy based platform (140) that we have developed specifically for processing, mapping and analysing Ribo-seq data. Researchers can use the GWIPS-viz suite of tools in RiboGalaxy to generate Ribo-seq profiles that infer either the A-site (elongating ribosomes) or P-site (initiating ribosomes) from either the 5 end or the 3 end of Ribo-seq reads and the resulting profiles can be directly visualised as custom tracks in GWIPS-viz. The direct interface between GWIPS-viz and RiboGalaxy also allows data from GWIPS-viz to be retrieved into RiboGalaxy. We also provide a direct link to RiboGalaxy (http://ribogalaxy.ucc.ie/) from the GWIPS-viz homepage.

Improvements in data visualizations
Previously for bacterial genomes (E. coli K12, B. subtilis) our Ribo-seq and mRNA-seq profiles on GWIPS-viz did not provide strand orientation information. Our Ribo-seq density plots also used the center-weighted approach (141) to infer ribosome A-sites. Since then, several studies have shown that inferring the ribosome decoding center from the 3 ends of bacterial Ribo-seq data is more accurate (124,142,143). We decided to carry-out an overhaul of our bacterial tracks and now provide strand orientation information using the UCSC Genome Browser overlay functionality (144) in addition to A-site inference using a fixed offset from 3 footprint ends ( Figure 1A,B). We have extended these improvements to the new bacterial genomes we now host in GWIPS-viz (E. coli BW25113, C. crescentus, S. coelicolor, S. aureus NCTC 8325, S. aureus USA300 FPR3757).
Recently we integrated the multi-region exon-only view (17), which is particularly useful for displaying Ribo-seq data for higher eukaryotes where exonic regions may be interrupted by long intronic regions ( Figure 1C,D). For higher eukaryotic datasets, we also now provide characteristics obtained with RUST (145). RUST utilizes Riboseq Unit Step Transfomation to normalize ribosome profiling data. It further provides characteristics of ribosome profiling datasets among which is a metafootprint profile which shows the difference between observed (experimental) and expected (equiprobable) frequencies of specific sequences (commonly codons) in the vicinity of a ribosome footprint. The expectation is that the highest variation in codon frequencies should occur at the ribosome decoding center (Asite) ( Figure 1E). A high variation at the end of footprints would occur due to sequencing biases. Thus, metafootprint profiles can be used for assessing the level of sequencing biases in individual datasets. Clicking on each study link in the GWIPS-viz genome page will open a new page with the link to the RUST quality plots which include the RUST metafootprint profile as well as a plot showing triplet periodicity for reads of different lengths. The RUST plots also include a panel that shows the relative inferred A-site dwell times for each amino acid.

Downloading Ribo-seq and mRNA-seq alignments
Following user requests, we added the functionality to download our genomic alignments in bigWig format. While the Table Browser provides the option to download our Ribo-seq and mRNA-seq alignments in bedGraph format, many users requested our original alignment files. Hence, we built a separate Downloads page ( Figure 1F) for this purpose. For each Ribo-seq study hosted on GWIPS-viz, users can download (1) ribosome profiles of elongating ribosomes (number of footprints whose inferred A-site match a specific coordinate, (2) Ribo-seq and (3) mRNA-seq coverage plots that provide the number of reads that map to each coordinate. Where available, data enriched with footprints of initiating ribosomes, represented as coordinates of inferred Psite codons, can also be downloaded. In addition, footprints of small ribosome subunits generated by TCP-seq (15) are available for download for S. cerevisiae as coverage plots.

FUTURE PLANS
The development of an automated computational data integration pipeline has greatly helped us to keep pace with the flux of new Ribo-seq data for genomes already existing in GWIPS-viz. We do, however, still have some backlog in terms of Ribo-seq data generated for genomes that we need to manually add to GWIPS-viz. We are examining ways in how we can improve our capacity to add new genomes particularly genomes that are not hosted on the UCSC Genome Browser.
We also aim to continue to improve the visualizations of Ribo-seq data. We wish to extend the overlay functionality with strand-specific display to all the data tracks in GWIPSviz (currently provided for bacterial genomes only). Conversely, we want to provide RUST characteristics for Riboseq data generated for bacteria. This requires adapting the RUST programming code to using 3 end offsetting for Asite inference which we plan to do soon. In addition, we want to use RUST parameters and other Ribo-seq specific parameters such as the triplet periodicity, to develop a quality scoring method and provide this information on GWIPS-viz for each dataset. This will allow users to compare Ribo-seq datasets in terms of their quality across studies. It will also help us and our user-base to determine what governs Ribo-seq data quality and how it can be maintained and improved in future Ribo-seq data.
We also wish to avail of new functionality as it becomes available on the UCSC Genome Browser. For this reason, we have carried out several upgrades of the GWIPS-viz browser to keep in-line with the UCSC Genome Browser. The recent exon-only view is one such example which is particularly beneficial for eukaryotic organisms that extensively use RNA splicing, permitting genomic alignments of Riboseq data to be explored without intervening introns.