A combinatorial and integrated method to analyse RNA-seq reads

Nicolas Philippe 1 Mikael Salson 2, 3 Thérèse Commes 4 Eric Rivals 1, *
* Auteur correspondant
1 MAB - Méthodes et Algorithmes pour la Bioinformatique
LIRMM - Laboratoire d'Informatique de Robotique et de Microélectronique de Montpellier
2 BONSAI - Bioinformatics and Sequence Analysis
LIFL - Laboratoire d'Informatique Fondamentale de Lille, Inria Lille - Nord Europe
Abstract : RNA sequencing enables a complete investigation covering the full dynamic spectrum of a tran- scriptome. It thus paves the way to a better understanding of the function of gene expression in different tissues, during development or pathological states. However, the splicing process, which generates both co-linear and non co-linear RNAs, the inclusion of sequencing errors, somatic muta- tions, polymorphisms, and rearrangements make the reads differ from the reference genome in a variety of ways. This complicates the task of comparing reads with a genome. Currently, the analysis paradigm consists in: 1. mapping the reads to a reference genome contiguously allowing as many differences as one expects to be necessary to accommodate sequence errors and small polymorphisms; 2. using uniquely mapped reads to determine covered genomic regions, either for computing a local coverage to predict mutations and filter out sequence errors (cf. program ERANGE), or for delimiting expressed exons approximately (cf. program TopHat); 3. re-aligning unmapped reads, which were not mapped contiguously at step one, to reveal splicing junctions. Limitations of this approach include lack of precision, redundant computations due to multi- mapping steps, error propagation due to heuristics and the absence of back-tracking. We propose a novel, integrated approach to analyze nowadays longer reads (> 50 bp). The idea is to adopt a k-mer approach that combines the genomic positions and local coverage to perform a complex analysis of each read and detect in a single step, mutations, indels, errors, as well as both normal and chimeric splice junctions. Comparisons with other tools demonstrate the feasibility of this ap- proach, which yields both sensitive and highly specific inferences.
Liste complète des métadonnées

Littérature citée [1 références]  Voir  Masquer  Télécharger

https://hal-lirmm.ccsd.cnrs.fr/lirmm-00757979
Contributeur : Eric Rivals <>
Soumis le : mardi 27 novembre 2012 - 18:04:02
Dernière modification le : jeudi 11 janvier 2018 - 06:26:13
Document(s) archivé(s) le : samedi 17 décembre 2016 - 16:26:51

Fichier

290-2698-1-PB.pdf
Fichiers éditeurs autorisés sur une archive ouverte

Identifiants

  • HAL Id : lirmm-00757979, version 1

Citation

Nicolas Philippe, Mikael Salson, Thérèse Commes, Eric Rivals. A combinatorial and integrated method to analyse RNA-seq reads. EMBnet.journal, EMBnet, 2011, 17 (Supplement B), pp.1. 〈http://journal.embnet.org/index.php/embnetjournal/article/view/290〉. 〈lirmm-00757979〉

Partager

Métriques

Consultations de la notice

513

Téléchargements de fichiers

269