Introduction à l'analyse des SNPs

Introduction

Un SNP (Single Nucleo­tid Poly­mor­phism) est la muta­tion d'un seul nucléo­tide à une posi­tion don­née, entre indi­vi­dus d'une même popu­la­tion. Ces sub­sti­tu­tions ponc­tuelles per­mettent d'identifier des sous-popu­la­tions.

Vous avez sûre­ment étu­dié en cours l'exemple de la dré­pa­no­cy­tose, une mala­die cau­sée par la muta­tion ponc­tuelle d'un gène. L'origine de cette mala­die géné­tique peut être obser­vée en ana­ly­sant les SNPs dans la popu­la­tion humaine. Ana­ly­ser ces der­niers peut ain­si per­mettre de trou­ver des pré­dis­po­si­tions à cer­taines mala­dies, mais peut éga­le­ment être utile entre autres, en géné­tique des popu­la­tions ou en can­cé­ro­lo­gie. L'arrivée des séquen­ceurs de deuxième géné­ra­tion (les fameux NGS) a faci­li­té leur étude.

Pour la suite, nous par­ti­rons du prin­cipe que vous êtes fami­liers avec les concepts liés à l'analyse de don­nées NGS (reads, map­ping, for­mat SAM/​BAM, fastq, pro­fon­deur de séquen­çage etc.) sinon la suite risque d'être un peu ardue.

L'analyse

Pour ana­ly­ser vos SNPs, je vous pro­pose le pipe­line ci-des­sous. Ce der­nier a pour but d'identifier les SNPs pré­sents dans vos échan­tillons puis de les fil­trer. À la fin de l'article nous ver­rons com­ment les anno­ter par rap­port à une base de don­nées.

L'analyse com­mence après ali­gne­ment de vos reads contre votre génome de réfé­rence via votre map­per pré­fé­ré. À la suite de cette étape vous obte­nez clas­si­que­ment un fichier au for­mat BAM. La pre­mière étape propre à l'analyse des SNPs consiste à uti­li­ser la com­mande

mpileup
, fournie avec l'outil
samtools
. Elle s'utilise de la manière suivante :
samtools mpileup -uf ref.fasta fichier.bam > output.bcf
D'autres options peuvent bien sûr être rajoutées (cf. la doc).
Si vous possédez plusieurs fichiers BAM (correspondant à plusieurs échantillons mappés contre le même génome), vous pouvez utiliser cette commande sur tous vos fichier d'un coup, il est d'ailleurs fortement conseillé de procéder ainsi. D'une part vous n'aurez qu'un fichier à gérer par la suite et d'autre part, l'algorithme de 
mpileup
utilise les reads de tous vos échantillons pour ses calculs. La commande 
mpileup
produit un fichier au format pileup (plus trop utilisé actuellement). Plus généralement, l'option
-u
est utilisée afin de produire un fichier binaire (format BCF). Ce dernier est le format compressé de VCF : Variant Call Format (en gros BCF est à VCF ce que BAM est à SAM).
 

 
Ci-dessus un exemple de fichier VCF. Ce format se compose de 3 parties :

  Les premières lignes, commençant par "##", contiennent la description des champs. Dans l'exemple ci-dessus on peut lire "##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Read Depth">" ce qui indique que le champ DP ne doit contenir qu'un seul chiffre ("Number=1") entier ("Type=Integer") et ce chiffre correspond à la profondeur du read ("Description="Reads Depth"").
L'en-tête (header) : Le nom des colonnes, commençant par "#"
Les données elles-mêmes contenant différents champs séparés par une tabulation.

 
Plus d'infos sur ce format sur la page wikipedia d'où provient la capture d'écran. La documentation officielle du format VCF est disponible ici.
Revenons à notre pipeline, 
mpileup
ne fait qu' "entasser" les reads alignés sur votre génome (du verbe anglais to pileup). Afin de supprimer les faux positifs dus à des erreurs de séquençage (c'est-à-dire un mismatch entre un read et le génome de référence sans pour autant que cela soit un SNP), il est possible d'utiliser l'outil
bcftools
. Celui-ci vous permettra également d'obtenir un fichier au format VCF qui sera utile pour la suite.
Par exemple, pour convertir le fichier
output.bcf
  obtenu précédemment, il suffit d'écrire :
bcftools view -v output.bcf > output.vcf
L'option
-v
sert à afficher uniquement les potentiels variants.
D'une manière générale 
bcftools
permet de manipuler les fichier VCF et BCF (
bcftools
  est aux formats VCF et BCF ce que 
samtools
est aux formats SAM et BAM).
À noter : 
La commande 
view
est faite pour fonctionner avec 
samtools
≤ 0.1.19 (la version de 
samtools
disponible dans les paquets Ubuntu à l'heure actuelle). Si vous avez une version de
samtools
et 
bcftools
plus récente il faudra passer par la commande
call
(l'étape de détection des SNPs se nomme SNP calling).  Ainsi la commande sera :
bcftools call -v -m output.bcf > output.vcf
Pourquoi
-m
? Eh ben d'après le site officiel de bcftools :
 Users are now required to choose between the old samtools calling model (-c/--consensus-caller) and the new multiallelic calling model (-m/--multiallelic-caller). The multiallelic calling model is recommended for most tasks.
L'option
-m
est donc la nouvelle recommandation des auteurs. L'option
-c
est plutôt utilisée pour rechercher des sites bialléliques.
La différence entre 
mpileup
et 
bcftools
peut vous sembler floue. Pour résumer, 
mpileup
 rassemble tous les reads à la même position sur le génome et
bcftools
 détecte les SNPs de manière plus fine en supprimant les faux positifs. Dans la pratique, ces deux commandes sont utilisées conjointement via le symbole pipe (
|
) sous Unix. C'est d'ailleurs cette façon de faire qui est indiquée dans la doc de mpileup.
Ainsi si vous ne deviez retenir qu'une seule commande, retenez celle-là :
samtools mpileup -uf ref.fasta fichier.bam  | bcftools call -v -m > output.vcf
 
Filtres
Maintenant que vous avez identifié vos SNPs, vous pouvez les compter, les regarder, mais vous pouvez aussi vous demander que faire (de plus) avec tous ces SNPs ? Comment ne garder que les plus "pertinents" ? (La difficulté ici est de définir "pertinent").
Vous pouvez filtrer vos SNPs selon :

Le nombre de fois qu'ils sont lus;
La qualité de l'alignement;
Le taux d’allèles à une position donnée

 
C'est par exemple ce qui est effectué dans Seo et. al (2012)  sur des SNPs humains (cf. paragraphe "Somatic single nucleotide and short-indel discovery").
Afin de vous aider à choisir vos filtres, vous pouvez vous poser les deux questions suivantes :

Avez-vous des contrôles positifs ? Si vous êtes sûrs que vous devez retrouver tel ou tel SNPs, est-ce que vous les retrouvez après avoir appliqué vos filtres ?;
Quel est le contexte de votre étude ? Est-ce que vous vous attendez à avoir "beaucoup" ou "peu" de SNPs ?

 
Pour filtrer vos SNPs, vous pouvez par exemple utiliser
snpsift
.
Utilisation :
java -jar SnpSift.jar filter [options] 'expression‘
Par exemple, dans la publication citée ci-dessus, les auteurs préconisent de ne garder que les SNPs lus avec une profondeur  ≥ 3. Ce qui se traduit par la commande suivante :
java -jar SnpSift.jar filter "(DP >= 3)" -f Outputfile.vcf > Outputfile_filter.vcf
 
Plus d'info sur http://snpeff.sourceforge.net/SnpSift.html#filter
Sachez que certains outils, tels que GATK, offrent des filtres plus fin pour analyser les SNPs. Par contre ce dernier est plus compliqué à utiliser : l'installation n'est pas simple et avant de commencer à analyser vos SNPs, prévoyez de nombreuses heures pour lire la doc.
L'annotation
Le pipeline présenté tout au long de cet article s'arrête après l'étape de filtrage. Ceci dit, une fois vos SNPs obtenus, vous pouvez également les comparer à des bases de données afin de les annoter.
SnpSift
permet de le faire grâce à la commande
annotate
. L'idée est de regarder si vos SNPs sont présents dans la base de données choisie et si oui de leur rajouter certaines informations. Prenons par exemple la base de données dbSNP qui est une base de données de SNPs humains hébergée par le NCBI. Mettons que si l'un de vos SNPs est présent chez dbSNP alors vous voulez rajouter l'information à propos de la fréquence de l'allèle basée sur l'étude 1 000 genomes (cette information est contenue dans le champ CAF d'après le glossaire). Pour réaliser ceci vous pouvez utiliser la commande suivante :
java -jar SnpSift.jar annotate 00-All.vcf myFile.vcf -info 'CAF' > annotatedSample.vcf
Dans la commande ci-dessus
00-All.vcf
correspond à la base de données.
Conclusion
Pour résumer, une analyse de SNPs commence toujours par un mapping, suivi d'un SNP calling (avec samtools et bcftools) afin d'obtenir un fichier VCF. Ce dernier peut éventuellement être filtré dans le but d'éliminer les allèles moins fréquents par exemple. Enfin, la partie annotation des SNPs est la plus compliquée car elle est dépendante des données et plus particulièrement du génome et du contexte de l'étude.
Comme évoqué dans l'introduction, les NGS ont grandement facilité l'étude des SNPs. La probabilité de détecter un SNP rare dépend du nombre d'individus séquencés. Ainsi, grâce au projet 100,000 genomes (séquençage de 100,000 génomes humains pour étudier des maladies rares) on peut s'attendre à découvrir de nouveaux SNPs chez l'humain. Pour finir, sachez que les SNPs sont le plus souvent étudiés dans le cadre des maladies génétiques mais rien ne vous empêche de les analyser dans d'autres domaines comme la phylogénie par exemple.
 
Merci à Julien Fouret pour son aide ainsi qu'à Lelouar, Bérénice, et NiGoPol pour les commentaires et discussions lors de l’édition de cet article.
      
          
        Vous avez aimé ? Dites-le nous !      
            
      
                  
                      
                  
                      
                  
                      
                  
                      
                  
                      
              
    
    
          Soumettre        
      Moyenne : 0 / 5. Nb de votes : 0    
    
      Pas encore de vote pour cet article.