Skip to content

Latest commit

 

History

History
138 lines (105 loc) · 6.49 KB

README.md

File metadata and controls

138 lines (105 loc) · 6.49 KB
Raw

VCF2Dis

VCF2Dis: A new simple and efficient software to calculate p-distance matrix and construct phylogeny based Variant Call Format input

1) Install and Parameter

    1. Install

The new version will be updated and maintained in hewm2008/VCF2Dis, please click below Link to download the latest version

hewm2008/VCF2Dis

Download


Just sh make.sh to compile. The executable VCF2Dis can be found in the folder of bin/VCF2Dis
For Linux /Unix and macOS

        tar -zxvf  VCF2DisXXX.tar.gz            # if Link do not work ,Try re-install [zlib]library
        cd VCF2DisXXX;                          # [zlib] and copy them to the library Dir
        sh make.sh;                             # VCF2Dis-xx/src/include/zlib
        ./bin/VCF2Dis

Note: If fail to link,try to re-install the libraries zlib
Note:: R with ape, dplyr and ggtree are recommended

    1. Main parameter description:
	Usage: VCF2Dis -InPut  <in.vcf>  -OutPut  <p_dis.mat>

		-InPut     <str>     Input one or muti GATK VCF genotype File
		-OutPut    <str>     OutPut Sample p-Distance matrix

		-InList    <str>     Input GATK muti-chr VCF Path List
		-SubPop    <str>     SubGroup SampleList of VCF File [ALLsample]
		-Rand      <float>   Probability (0-1] for each site to join Calculation [1]

		-help                Show more help [hewm2008 v1.53s]


For more details, please use -help and see the example

		-InFormat      <str>   Input File is [VCF/FA/PHY] Format,defaut: [VCF]
		-InSampleGroup <str>   InFile of sample Group info,format(sample groupA)
		-TreeMethod    <int>   Construct Tree Method,1:NJ-tree 2:UPGMA-tree [1]
		-KeepMF                Keep the Middle File diff & Use matrix

2) Example


Three examples were provided in the directory of example/Example*

1) an Example of nj-tree with no boostrap

    1. To Create the p_distance matrix and construct nj-tree newick tree
# 1.1) To new all the sample p_distance matrix and newick tree based VCF, run VCF2Dis directly
      ./bin/VCF2Dis	-InPut	in.vcf.gz	-OutPut p_dis.mat
      #  ./bin/VCF2Dis     -InPut  in.fa.gz	-OutPut p_dis.mat -InFormat  FA

# 2.2) To new sub group sample p_distance matrix and and newick tree ; put their sample name into File sample.list
      ./bin/VCF2Dis	-InPut	chr1.vcf.gz chr2.vcf.gz	-OutPut p_dis.mat  -SubPop  sample.list
    1. Simple tree visualization (for advanced tree display and annotation please refer to iTOL, Evolview, MEGA)
      you will obtain the p_dis.nwk tree file and neighbor-joining tree in pdf format p_dis.pdf after VCF2Dis.


Note::if you can't get the p_dis.nwk tree file but had the p_dis.mat, here are the 3 methods to get the tree file.

2) an Example of nj-tree with boostrap

    1. Running multiple times by using a method of sampling with replacement. Users can randomly select a part of the sites [-Rand] and construct a new nj-tree as above, and Repeat NN times [recommand NN=100]. X=(1,2....NN);
#!/bin/bash
NN=100
if [ "$#"  -eq  1 ]; then
	NN=$1
fi
for X in $(seq 1 $NN)
do
	./bin/VCF2Dis -InPut in.vcf.gz -OutPut p_dis_${X}.mat -Rand 0.25
    # PHYLIPNEW-3.69.650/bin/fneighbor -datafile p_dis_${X}.mat -outfile tree.out1_${X}.txt -matrixtype s -treetype n -outtreefile tree.out2_${X}.tre
done
    1. Merge all the nj-tree and construct and display a boostrap nj-tree. (For advanced display tree and annotation please refer to iTOL, Evolview and MEGA)
#!/bin/bash
NN=100
if [ "$#"  -eq  1 ]; then
  NN=$1
fi

cat  p_*.nwk  >    alltree_merge.tre   #  cat  tree*.tre  > alltree_merge.tre
PHYLIPNEW-3.69.650/bin/fconsense   -intreefile   alltree_merge.tre  -outfile out  -treeprint Y
perl  ./bin/percentageboostrapTree.pl    alltree_merge.treefile    $NN    Final_boostrap.tre  # NN is the input number

How to Install PHYLIPNEW please Click on here or Click on here(Chinese)

4) Introduction

The formula for calculating p-distance between indivisuals from VCF SNP datasets was listed below:

            D_ij=(1/L) * [(sum(d(l)_ij))]


Where L is the length of regions where SNPs can be identified, and given the alleles at position l are A/C between sample i and sample j:

            d(l)_ij=0.0     if the genotypes of the two individuals were AA and AA;
            d(l)_ij=0.5     if the genotypes of the two individuals were AA and AC;
            d(l)_ij=0.0     if the genotypes of the two individuals were AC and AC;
            d(l)_ij=1.0     if the genotypes of the two individuals were AA and CC;
            d(l)_ij=0.0     if the genotypes of the two individuals were CC and CC;

To further know about the p_distance matrix based the VCF file, please refer to this website.

5) Results

VCF2Dis have been cited in more than 150 times by searching against google scholar.
Below were some NJ-tree images that I draw in the paper before.

  • 50 Rices NBT
  • 31 soybeans NG
    Display tree by MAGA after test Data VCF2Dis -i ALL.chr*.genotypes.vcf.gz -SubPop subsample203.list -InSampleGroup pop.info

example1.png

6) Discussing

######################swimming in the sky and flying in the sea ########################### ##