Team III Gene Prediction Group

Introduction

Gene prediction is the process of identifying the regions of genomic DNA that encode genes which primarily include protein-coding and non-coding genes. Gene prediction is an important process that aids in the identification of fundamental and essential elements of the genome. With the overall goal to investigate a foodborne outbreak caused by a prokaryotic organism, our team developed a pipeline for the prediction of coding and non-coding genes in prokaryotes. Our final objective was to carry out a thorough and exhaustive prediction of all coding and non-coding genes of the 50 assembled genomes provided by the Genome Assembly team. In prokaryotic genomes, DNA sequences that encode proteins are transcribed into mRNA, and then RNA is usually translated directly into proteins without significant modification. They have a higher gene density in comparison to eukaryotes.

There are two gene prediction methods - Homology methods and ab-initio methods. Homology-based gene prediction methods rely on extrinsic evidence and makes predictions by comparing with sequences of previously known genes. Homology-based approaches are helpful as they are used to validate ab-initio methods but are limited by existing knowledge, sometimes computationally expensive, and require an extensive and vast database as all genes are not expressed at the same time. Ab-initio gene prediction methods, on the other hand, rely on intrinsic evidence of the genome. These methods detect promoter sequences, start and stop codons and GC content to predict ORFs. The disadvantages of ab-initio methods are the high false-positive rate, no experimental verification and not robust as homology-based tools to name a few. The main algorithms used in various ab-initio tools are Hidden Markov Modelling (HMM), Interpolated Markov Modelling (IMM) and Dynamic Programming.

Non-Coding homology - Tools

ARAGORN

Aragorn is a computer program identifies tRNA and tmRNA genes. The program employs heuristic algorithms to predict tRNA secondary structure, based on homology with recognized tRNA consensus sequences and ability to form a base-paired cloverleaf. tmRNA genes are identified using a modified version of the BRUCE program.

Infernal

Infernal is for searching DNA sequence databases for RNA structure and sequence similarities. It is an implementation of a special case of profile stochastic context-free grammars called covariance models (CMs). A CM is like a sequence profile, but it scores a combination of sequence consensus and RNA secondary structure consensus, so in many cases, it is more capable of identifying RNA homologs that conserve their secondary structure more than their primary sequence.

Non-Coding Ab initio - Tools

RNAmmer

RNAmmer predicts ribosomal RNA genes in full genome sequences by utilizing two levels of Hidden Markov Models: An initial spotter model searches both strands. The spotter model is constructed from highly conserved loci within a structural alignment of known rRNA sequences.

Barrnap

Barrnap predicts the location of ribosomal RNA genes in genomes by using HMMER 3.1 for HMM searching in RNA:DNA style. It supports bacteria (5S,23S,16S), archaea (5S,5.8S,23S,16S), metazoan mitochondria (12S,16S) and eukaryotes (5S,5.8S,28S,18S).

Benchmarking of non-coding + Results inferred

rRNA Comparison

tRNA Comparison

Results: Non-Coding

Aragorn & RNAmmer

$ aragorn -l -t -gc1 -w input.fasta -fo –o output.fasta
$ aragorn -l -m -gc1 -w input.fasta -fo –o output.fasta

• Average of tRNA: 40.9
• Average of tmRNA: 1
• Average of rRNA: 2.2

Infernal

$ cmscan --cut_ga --rfam --nohmmonly --tblout $output/$(basename $filename .fasta).tblout --fmt 2 --clanin Rfam.clanin Rfam.cm $filename > $output/$(basename $filename . fasta).cmscan

• Average of tRNA: 50.5
• Average of tmRNA: 1
• Average of rRNA: 3.34