SNP calling from RNA-seq data without a reference genome: Identification, quantification, differential analysis and impact on the protein sequence
Carregando...
Arquivos
Data
2016-11-02
Autores
Lopez-Maestre, Hélène
Brinza, Lilia
Marchet, Camille
Kielbassa, Janice
Bastien, Sylvère
Boutigny, Mathilde
Monnin, David
Filali, Adil El
Carareto, Claudia Marcia [UNESP]
Vieira, Cristina
Título da Revista
ISSN da Revista
Título de Volume
Editor
Resumo
SNPs (Single Nucleotide Polymorphisms) are genetic markers whose precise identification is a prerequisite for association studies. Methods to identify them are currently well developed for model species, but rely on the availability of a (good) reference genome, and therefore cannot be applied to non-model species. They are also mostly tailored for whole genome (re-)sequencing experiments, whereas in many cases, transcriptome sequencing can be used as a cheaper alternative which already enables to identify SNPs located in transcribed regions. In this paper, we propose a method that identifies, quantifies and annotates SNPs without any reference genome, using RNA-seq data only. Individuals can be pooled prior to sequencing, if not enough material is available from one individual. Using pooled human RNA-seq data, we clarify the precision and recall of our method and discuss them with respect to other methods which use a reference genome or an assembled transcriptome. We then validate experimentally the predictions of our method using RNA-seq data from two non-model species. The method can be used for any species to annotate SNPs and predict their impact on the protein sequence. We further enable to test for the association of the identified SNPs with a phenotype of interest.
Descrição
Palavras-chave
Como citar
Nucleic Acids Research, v. 44, n. 19, 2016.