Convert SNPs in VCF format to PHYLIP, NEXUS, binary NEXUS, or FASTA alignments for phylogenetic analysis
MIN_SAMPLES_LOCUS：在每一个record中，基因型不为./.的样本数目，例如1/1、./.、./.、1/1、1/1五个样本中有2个为./.，这样的话MIN_SAMPLES_LOCUS就是3小于minimum for phylogenetics的4，因此就被忽略了- 该脚本具体的做法：如果基因型为
1/1，那么就将该样本对应的fasta中写入alt对应的AGCT；如果是0/1(ref/alt=A/G)，则改位点可能是A、G，按照IUPAC的对应关系将会写入R；多个样本多行按照这个规则写入。- 程序使用的IUPAC的规则见下
- 原始的IUPAC规则可以参考这个网站
# Dictionary of IUPAC ambiguities for nucleotides
# '*' means deletion for GATK (and other software?)
# Deletions are ignored when making the consensus
# Dictionary to translate IUPAC ambiguities, lowercase letters are used when "*" or "N" were present for a position,
# however, software like Genious for example are case insensitive and will imply ignore capitalization
amb = {"*":"-", "A":"A", "C":"C", "G":"G", "N":"N", "T":"T",
"*A":"a", "*C":"c", "*G":"g", "*N":"n", "*T":"t",
"AC":"M", "AG":"R", "AN":"a", "AT":"W", "CG":"S",
"CN":"c", "CT":"Y", "GN":"g", "GT":"K", "NT":"t",
"*AC":"m", "*AG":"r", "*AN":"a", "*AT":"w", "*CG":"s",
"*CN":"c", "*CT":"y", "*GN":"g", "*GT":"k", "*NT":"t",
"ACG":"V", "ACN":"m", "ACT":"H", "AGN":"r", "AGT":"D",
"ANT":"w", "CGN":"s", "CGT":"B", "CNT":"y", "GNT":"k",
"*ACG":"v", "*ACN":"m", "*ACT":"h", "*AGN":"r", "*AGT":"d",
"*ANT":"w", "*CGN":"s", "*CGT":"b", "*CNT":"y", "*GNT":"k",
"ACGN":"v", "ACGT":"N", "ACNT":"h", "AGNT":"d", "CGNT":"b",
"*ACGN":"v", "*ACGT":"N", "*ACNT":"h", "*AGNT":"d", "*CGNT":"b",
"*ACGNT":"N"}This script takes as input a VCF file and will use the SNP genotypes to create a matrix for phylogenetic analysis in the PHYLIP (relaxed version), FASTA, NEXUS, or binary NEXUS formats. For heterozygous SNPs the consensus is made and the IUPAC nucleotide ambiguity codes are written to the final matrix(ces), any ploidy level is allowed and automatically detected. The code is optimized for large VCF matrices (hundreds of samples and millions of genotypes), for example, in our tests it processed a 20GB VCF (~3 million SNPs x 650 individuals) in ~27 minutes. The initial version of the script just produced a PHYLIP matrix but now we have added other popular formats, including the binary NEXUS file to run SNPs analysis with the SNAPP plugin in BEAST (only for diploid samples).
Additionally, you can choose a minimum number of samples per SNP to control the final amount of missing data. Since phylogenetic software usually root the trees at the first sequence in the alignment (e.g. RAxML, IQTREE, and MrBayes), the script also allows you to specify an OUTGROUP sequence that will be written in the first place in the alignment.
Compressed VCF files can be directly analyzed but the extension must be .vcf.gz.
The script has been tested with VCF files produced by pyrad v.3.0.66, ipyrad v.0.7.x, Stacks v.1.47, dDocent, GATK, and freebayes.
Please don't hesitate to open an Issue if you find any problem or suggestions for a new feature.
Just type python vcf2phylip.py -h to show the help of the program:
usage: vcf2phylip.py [-h] -i FILENAME [-m MIN_SAMPLES_LOCUS] [-o OUTGROUP]
[-p] [-f] [-n] [-b]
The script converts a collection of SNPs in VCF format into a PHYLIP, FASTA,
NEXUS, or binary NEXUS file for phylogenetic analysis. The code is optimized
to process VCF files with sizes >1GB. For small VCF files the algorithm slows
down as the number of taxa increases (but is still fast).
Any ploidy is allowed, but binary NEXUS is produced only for diploid VCFs.
optional arguments:
-h, --help show this help message and exit
-i FILENAME, --input FILENAME
Name of the input VCF file, can be gzipped
-m MIN_SAMPLES_LOCUS, --min-samples-locus MIN_SAMPLES_LOCUS
Minimum of samples required to be present at a locus,
default=4 since is the minimum for phylogenetics.
-o OUTGROUP, --outgroup OUTGROUP
Name of the outgroup in the matrix. Sequence will be
written as first taxon in the alignment.
-p, --phylip-disable A PHYLIP matrix is written by default unless you
enable this flag
-f, --fasta Write a FASTA matrix, disabled by default
-n, --nexus Write a NEXUS matrix, disabled by default
-b, --nexus-binary Write a binary NEXUS matrix for analysis of biallelic
SNPs in SNAPP, disabled by default
In the following examples you can omit python if you change the permissions of vcf2phylip.py to executable.
Example 1: Use default parameters to create a PHYLIP matrix with a minimum of 4 samples per SNP:
python vcf2phylip.py --input myfile.vcf
# Which i equivalent to:
python vcf2phylip.py -i myfile.vcf
# This command will create a PHYLIP called myfile_min4.phyExample 2: Create a PHYLIP and a FASTA matrix using a minimum of 60 samples per SNP:
python vcf2phylip.py --input myfile.vcf --fasta --min-samples-locus 60
# Which is equivalent to:
python vcf2phylip.py -i myfile.vcf -f -m 60
# This command will create a PHYLIP called myfile_min60.phy and a FASTA called myfile_min60.fastaExample 3: Create all output formats, and select "sample1" as outgroup:
python vcf2phylip.py --input myfile.vcf --outgroup sample1 --fasta --nexus --nexus-binary
# Which is equivalent to:
python vcf2phylip.py -i myfile.vcf -o sample1 -f -n -b
# This command will create a PHYLIP called myfile_min4.phy, a FASTA called myfile_min4.fasta, a NEXUS called myfile_min4.nexus, and a binary NEXUS called myfile_min4.bin.nexusExample 4: If, for example, you wish to disable the creation of the PHYLIP matrix and only create a NEXUS matrix:
python vcf2phylip.py --input myfile.vcf --phylip-disable --nexus
# Which is equivalent to:
python vcf2phylip.py -i myfile.vcf -p -n
# This command will create only a NEXUS matrix called myfile_min4.nexus- Code: Edgardo M. Ortiz
- Data and testing: Juan D. Palacio-Mejía
Ortiz, E.M. 2019. vcf2phylip v2.0: convert a VCF matrix into several matrix formats for phylogenetic analysis. DOI:10.5281/zenodo.2540861