Improving the Caenorhabditis elegans Genome Annotation Using Machine Learning

For modern biology, precise genome annotations are of prime importance, as they allow the accurate definition of genic regions. We employ state-of-the-art machine learning methods to assay and improve the accuracy of the genome annotation of the nematode Caenorhabditis elegans. The proposed machine learning system is trained to recognize exons and introns on the unspliced mRNA, utilizing recent advances in support vector machines and label sequence learning. In 87% (coding and untranslated regions) and 95% (coding regions only) of all genes tested in several out-of-sample evaluations, our method correctly identified all exons and introns. Notably, only 37% and 50%, respectively, of the presently unconfirmed genes in the C. elegans genome annotation agree with our predictions, thus we hypothesize that a sizable fraction of those genes are not correctly annotated. A retrospective evaluation of the Wormbase WS120 annotation [1] of C. elegans reveals that splice form predictions on unconfirmed genes in WS120 are inaccurate in about 18% of the considered cases, while our predictions deviate from the truth only in 10%–13%. We experimentally analyzed 20 controversial genes on which our system and the annotation disagree, confirming the superiority of our predictions. While our method correctly predicted 75% of those cases, the standard annotation was never completely correct. The accuracy of our system is further corroborated by a comparison with two other recently proposed systems that can be used for splice form prediction: SNAP and ExonHunter. We conclude that the genome annotation of C. elegans and other organisms can be greatly enhanced using modern machine learning technology

Automated eukaryotic gene structure annotation using EVidenceModeler and the Program to Assemble Spliced Alignments

EVidenceModeler (EVM) is an automated annotation tool that predicts protein-coding regions, alternatively spliced transcripts and untranslated regions of eukaryotic genes

A first glimpse at the transcriptome of Physarum polycephalum

This is an Open Access article distributed under the terms of the Creative Commons Attribution Licens

A Machine Learning Model for Discovery of Protein Isoforms as Biomarkers

Prostate cancer is the most common cancer in men. One in eight Canadian men will be diagnosed with prostate cancer in their lifetime. The accurate detection of the disease’s subtypes is critical for providing adequate therapy; hence, it is critical for increasing both survival rates and quality of life. Next generation sequencing can be beneficial when studying cancer. This technology generates a large amount of data that can be used to extract information about biomarkers. This thesis proposes a model that discovers protein isoforms for different stages of prostate cancer progression. A tool has been developed that utilizes RNA-Seq data to infer open reading frames (ORFs) corresponding to transcripts. These ORFs are used as features for classificatio. A quantification measurement, Adaptive Fragments Per Kilobase of transcript per Million mapped reads (AFPKM), is proposed to compute the expression level for ORFs. The new measurement considers the actual length of the ORF and the length of the transcript. Using these ORFs and the new expression measure, several classifiers were built using different machine learning techniques. That enabled the identification of some protein isoforms related to prostate cancer progression. The biomarkers have had a great impact on the discrimination of prostate cancer stages and are worth further investigation

Deep RNA sequencing analysis of readthrough gene fusions in human prostate adenocarcinoma and reference samples

<p>Abstract</p> <p>Background</p> <p>Readthrough fusions across adjacent genes in the genome, or transcription-induced chimeras (TICs), have been estimated using expressed sequence tag (EST) libraries to involve 4-6% of all genes. Deep transcriptional sequencing (RNA-Seq) now makes it possible to study the occurrence and expression levels of TICs in individual samples across the genome.</p> <p>Methods</p> <p>We performed single-end RNA-Seq on three human prostate adenocarcinoma samples and their corresponding normal tissues, as well as brain and universal reference samples. We developed two bioinformatics methods to specifically identify TIC events: a targeted alignment method using artificial exon-exon junctions within 200,000 bp from adjacent genes, and genomic alignment allowing splicing within individual reads. We performed further experimental verification and characterization of selected TIC and fusion events using quantitative RT-PCR and comparative genomic hybridization microarrays.</p> <p>Results</p> <p>Targeted alignment against artificial exon-exon junctions yielded 339 distinct TIC events, including 32 gene pairs with multiple isoforms. The false discovery rate was estimated to be 1.5%. Spliced alignment to the genome was less sensitive, finding only 18% of those found by targeted alignment in 33-nt reads and 59% of those in 50-nt reads. However, spliced alignment revealed 30 cases of TICs with intervening exons, in addition to distant inversions, scrambled genes, and translocations. Our findings increase the catalog of observed TIC gene pairs by 66%.</p> <p>We verified 6 of 6 predicted TICs in all prostate samples, and 2 of 5 predicted novel distant gene fusions, both private events among 54 prostate tumor samples tested. Expression of TICs correlates with that of the upstream gene, which can explain the prostate-specific pattern of some TIC events and the restriction of the <it>SLC45A3-ELK4 </it>e4-e2 TIC to <it>ERG</it>-negative prostate samples, as confirmed in 20 matched prostate tumor and normal samples and 9 lung cancer cell lines.</p> <p>Conclusions</p> <p>Deep transcriptional sequencing and analysis with targeted and spliced alignment methods can effectively identify TIC events across the genome in individual tissues. Prostate and reference samples exhibit a wide range of TIC events, involving more genes than estimated previously using ESTs. Tissue specificity of TIC events is correlated with expression patterns of the upstream gene. Some TIC events, such as <it>MSMB-NCOA4</it>, may play functional roles in cancer.</p

Transcript Diversity In The Protozoan Parasite Toxoplasma Gondii

Technological advances have made possible to sequence RNA transcripts at unprecedented depth, enabling deep profiling of abundance and diversity under a variety of conditions. Such information permits refinement of draft genome annotation originally generated in the absence of transcript coverage data, and provides new insights into organismal biology and regulatory mechanisms. This dissertation provides an extensive analysis of mRNA-seq data from the obligate intracellular protozoan parasite Toxoplasma gondii, a ubiquitous pathogen of humans and other vertebrates. We produced and sequenced 24 strand-specific RNA libraries from several parasite strains and developmental stages, and examined these in con�junction with 45 additional mRNA-seq libraries produced by other groups. The current reference genome annotation for T. gondii, generated using de novo methods informed by cDNA sequencing prior to mRNA-seq, identifies ~8300 protein-coding genes, fragmented by ~40K introns. Untranslated regions are incompletely defined, few alternatively-spliced transcripts are described, and non-coding transcripts remain largely unexplored. mRNA-seq datasets presented in this dissertation define a total of 2.7M introns, most observed at vanishingly low abundance. Using current annotation to define parameters minimizing false discovery yields ~60K likely splice junctions. Comparing the frequency of intron-spanning reads to the abundance of transcripts to which introns belong provides a reliable metric for estimating intron excision, readily distinguishing introns that are (i) universally used, (ii) alternatively-spliced, or (iii) likely insignificant. Genome-wide analysis suggests ~3000 annotated introns that should be deleted from the reference genome, ~1400 to be added as alternative isoforms, ~3100 as additions to existing annotation (often within UTRs) and ~3400 associated with novel transcripts. Transcriptomic expression is consis�tent with biological and phenotypic variation across the complex parasite life cycle, including undescribed differences in gene expression during intracellular tachyzoite replication. Strong circumstantial evidence also suggests that lncRNAs may play an important role in regulating stage-specific expression during sexual differentia�tion and sporogony. These results provide the basis for revising the reference T. gondii genome annotation available at ToxoDB.org and GenBank. Strategies developed in this dissertation also provide the basis for defining annotation criteria for other species, including related parasites responsible for malaria and conceivably other eukaryotes as well

Spliced alignment and its application in Arabidopsis thaliana

This thesis describes the development and biological applications of GeneSeqer, which is a homology-based gene prediction program by means of spliced alignment. Additionally, a program named MyGV was written in JAVA as a browser to visualize the output of GeneSeqer. In order to test and demonstrate the performance, GeneSeqer was utilized to map 176,915 Arabidopsis EST sequences on the whole genome of Arabidopsis thaliana, which consists of five chromosomes, with about 117 million base pairs in total. All results were parsed and imported into a MySQL database. Information that was inferred from the Arabidopsis spliced alignment results may serve as valuable resource for a number of projects of special scientific interest, such as alternative splicing, non-canonical splice sites, mini-exons, etc. We also built AtGDB (Arabidopsis thaliana Genome DataBase, http://www.plantgdb.org/AtGDB/) to interactively browse EST spliced alignments and GenBank annotations for the Arabidopsis genome. Moreover, as one application of the Arabidopsis EST mapping data, U12-type introns were identified from the transcript-confirmed introns in the Arabidopsis genome, and the characteristics of these minor class introns were further explored

Rapidly evolving protointrons in Saccharomyces genomes revealed by a hungry spliceosome.

Introns are a prevalent feature of eukaryotic genomes, yet their origins and contributions to genome function and evolution remain mysterious. In budding yeast, repression of the highly transcribed intron-containing ribosomal protein genes (RPGs) globally increases splicing of non-RPG transcripts through reduced competition for the spliceosome. We show that under these "hungry spliceosome" conditions, splicing occurs at more than 150 previously unannotated locations we call protointrons that do not overlap known introns. Protointrons use a less constrained set of splice sites and branchpoints than standard introns, including in one case AT-AC in place of GT-AG. Protointrons are not conserved in all closely related species, suggesting that most are not under positive selection and are fated to disappear. Some are found in non-coding RNAs (e. g. CUTs and SUTs), where they may contribute to the creation of new genes. Others are found across boundaries between noncoding and coding sequences, or within coding sequences, where they offer pathways to the creation of new protein variants, or new regulatory controls for existing genes. We define protointrons as (1) nonconserved intron-like sequences that are (2) infrequently spliced, and importantly (3) are not currently understood to contribute to gene expression or regulation in the way that standard introns function. A very few protointrons in S. cerevisiae challenge this classification by their increased splicing frequency and potential function, consistent with the proposed evolutionary process of "intronization", whereby new standard introns are created. This snapshot of intron evolution highlights the important role of the spliceosome in the expansion of transcribed genomic sequence space, providing a pathway for the rare events that may lead to the birth of new eukaryotic genes and the refinement of existing gene function