Visualising ribosome profiling and using it for reading frame detection and exploration of eukaryotic translation initiation

Ribosome profiling (ribo-seq) is a recently developed technique that provides genomewide information on protein synthesis (GWIPS) in vivo. The high resolution of ribo-seq is one of the exciting properties of this technique. In Chapter 2, I present a computational method that utilises the sub-codon precision and triplet periodicity of ribosome profiling data to detect transitions in the translated reading frame. Application of this method to ribosome profiling data generated for human HeLa cells allowed us to detect several human genes where the same genomic segment is translated in more than one reading frame. Since the initial publication of the ribosome profiling technique in 2009, there has been a proliferation of studies that have used the technique to explore various questions with respect to translation. A review of the many uses and adaptations of the technique is provided in Chapter 1. Indeed, owing to the increasing popularity of the technique and the growing number of published ribosome profiling datasets, we have developed GWIPS-viz (http://gwips.ucc.ie), a ribo-seq dedicated genome browser. Details on the development of the browser and its usage are provided in Chapter 3. One of the surprising findings of ribosome profiling of initiating ribosomes carried out in 3 independent studies, was the widespread use of non-AUG codons as translation initiation start sites in mammals. Although initiation at non-AUG codons in mammals has been documented for some time, the extent of non-AUG initiation reported by these ribo-seq studies was unexpected. In Chapter 4, I present an approach for estimating the strength of initiating codons based on the leaky scanning model of translation initiation. Application of this approach to ribo-seq data illustrates that initiation at non-AUG codons is inefficient compared to initiation at AUG codons. In addition, our approach provides a probability of initiation score for each start site that allows its strength of initiation to be evaluated.

The relevance of bacterial isoprenoid biosynthetic pathways for host-microbe interactions

This thesis was undertaken to investigate the relevance of two bacterial isoprenoid biosynthetic pathways (Mevalonate (MVAL) and 2-C-methyl-D-erythritol 4-phosphate (MEP)) for host-microbe interactions. We determined a significant reduction in microbial diversity in the murine gut microbiota (by next generation sequencing) following oral administration of a common anti-cholesterol drug Rosuvastatin (RSV) that targets mammalian and bacterial HMG-CoA reductase (HMG-R) for inhibition of MVAL formation. In tandem we identified significant hepatic and intestinal off-target alterations to the murine metabolome indicating alterations in inflammation, bile acid profiles and antimicrobial peptide synthesis with implications on community structure of the gastrointestinal microbiota in statin-treated animals. However we found no effect on local Short Chain Fatty Acid biosynthesis (metabolic health marker in our model). We demonstrated direct inhibition of bacterial growth in-vitro by RSV which correlated with reductions in bacterial MVAL formation. However this was only at high doses of RSV. Our observations demonstrate a significant RSV-associated impact on the gut microbiota prompting similar human analysis. Successful deletion of another MVAL pathway enzyme (HMG-CoA synthase (mvaS)) involved in Listeria monocytogenes EGDe isoprenoid biosynthesis determined that the enzyme is non-essential for normal growth and in-vivo pathogenesis of this pathogen. We highlight potential evidence for alternative means of synthesis of the HMG-CoA substrate that could render mvaS activity redundant under our test conditions. Finally, we showed by global gene expression analysis (Massive Analysis of cDNA Ends (MACE RNA-seq) a significant role for the penultimate MEP pathway metabolite (E)-4-hydroxy-3-methyl-2-but-2-enyl pyrophosphate (HMBPP) in significant up regulation of genes of immunity and antigen presentation in THP-1 cells at nanomolar levels. We infected THP-1 cells with wild type or HMBPP under/over-producing L. monoctyogenes EGDe mutants and determined subtle effects of HMBPP upon overall host responses to Listeria infection. Overall our findings provide greater insights regarding bacterial isoprenoid biosynthetic pathways for host-microbe/microbe-host dialogue.

Comparative survey of the relative impact of mRNA features on local ribosome profiling read density

Ribosome profiling (Ribo-seq), a promising technology for exploring ribosome decoding rates, is characterized by the presence of infrequent high peaks in ribosome footprint density and by long alignment gaps. Here, to reduce the impact of data heterogeneity we introduce a simple normalization method, Ribo-seq Unit Step Transformation (RUST). RUST is robust and outperforms other normalization techniques in the presence of heterogeneous noise. We illustrate how RUST can be used for identifying mRNA sequence features that affect ribosome footprint densities globally. We show that a few parameters extracted with RUST are sufficient for predicting experimental densities with high accuracy. Importantly the application of RUST to 30 publicly available Ribo-seq data sets revealed a substantial variation in sequence determinants of ribosome footprint frequencies, questioning the reliability of Ribo-seq as an accurate representation of local ribosome densities without prior quality control. This emphasizes our incomplete understanding of how protocol parameters affect ribosome footprint densities.

Identification of RNA editing in the human exome and development of DARNED DatabaseSF

RNA editing is a biological phenomena that alters nascent RNA transcripts by insertion, deletion and/or substitution of one or a few nucleotides. It is ubiquitous in all kingdoms of life and in viruses. The predominant editing event in organisms with a developed central nervous system is Adenosine to Inosine deamination. Inosine is recognized as Guanosine by the translational machinery and reverse-transcriptase. In primates, RNA editing occurs frequently in transcripts from repetitive regions of the genome. In humans, more than 500,000 editing instances have been identified, by applying computational pipelines on available ESTs and high-throughput sequencing data, and by using chemical methods. However, the functions of only a small number of cases have been studied thoroughly. RNA editing instances have been found to have roles in peptide variants synthesis by non-synonymous codon substitutions, transcript variants by alterations in splicing sites and gene silencing by miRNAs sequence modifications. We established the Database of RNA EDiting (DARNED) to accommo-date the reference genomic coordinates of substitution editing in human, mouse and fly transcripts from published literatures, with additional information on edited genomic coordinates collected from various databases e.g. UCSC, NCBI. DARNED contains mostly Adenosine to Inosine editing and allows searches based on genomic region, gene ID, and user provided sequence. The Database is accessible at http://darned.ucc.ie RNA editing instances in coding region are likely to result in recoding in protein synthesis. This encouraged me to focus my research on the occurrences of RNA editing specific CDS and non-Alu exonic regions. By applying various filters on discrepancies between available ESTs and their corresponding reference genomic sequences, putative RNA editing candidates were identified. High-throughput sequencing was used to validate these candidates. All predicted coordinates appeared to be either SNPs or unedited.