Note this software is not being maintained any longer.

IDPASE

Software is to be run on a linux system. To run IDPASE

1. Install julia
2. From julia's command line run

Pkg.clone("git://github.com/bdeonovic/IDPASE.jl.git")

3. Install Bedtools and make sure it is in your PATH.
4. Prepare necessary input files

• VCF file
• GPD file in the Extended format
• PSL alignment file of your hybrid-Seq data
• FASTQ file of your hybrid-Seq data

5. Prepare Gene level data

mkdir temp/; mkdir gene_files; mkdir isoform_files; mkdir gene_out; mkdir isoform_out;
julia ~/.julia/v0.5/IDPASE/src/prep_runs.jl -a ~/.julia/v0.5/IDPASE/test_data/SGS.psl ~/.julia/v0.5/IDPASE/test_data/TGS.psl -g ~/.julia/v0.5/IDPASE/test_data/TDRKH.gpd -v ~/.julia/v0.5/IDPASE/test_data/sim.vcf -q ~/.julia/v0.5/IDPASE/test_data/SGS.fq ~/.julia/v0.5/IDPASE/test_data/TGS.fq -d temp/ -c 1 -f 1 1 -o gene_files/ -p sim

where flag -a is space separated list of PSL files, -g is GPD file, -v is VCF file, -q is space separated list of FASTQ files, -d is a directory for intermediate output, -c is a space separated list of chromosomes of interest, -f is a space separated list of FASTQ formats corresponding to the FASTQ files of -q, where 1 indicates PHRED+33, and 2 indicates PHRED+64, and -o indicates output directory and output prefix (so in example /out/ is directory and output files will be prefixed by sim). Make sure the directories specified exist such as ~/temp and ~/gene_files. 6. Get commands to run each gene individually

julia  ~/.julia/v0.5/IDPASE/scripts/phase_by_loci.jl -a gene_files/ -o gene_out/ -n SGS TGS -m 1 0 0 1 1 1 -d ~/.julia/v0.5/IDPASE/scripts/ -p sim > to_run_curr.sh

where -a is the -o flag from command in step 5, -o is an output directory, -n are unique names corresponding to PSL files, -m is a vector indicating which combinations of the seq data to use with IDPASE. In the above example 3 runs of IDPASE will be issued where 1 0 indicates SGS only, 0 1 indicates TGS only, and 1 1 indicates hybrid-Seq. p is the prefix specified in step 5. The output is a list of commands to run for each gene individually. The flag -d is the directory where the IDPASE scripts are stored. 7. Run each gene.

bash to_run_curr.sh

8. Concatenate all gene level results

 find gene_out/ -name "REAL*" | xargs cat > gene_out/gene_results.txt

The output columns are Number of input files, comma seperated list of read counts from each input file, total reads used in analysis, number of SNV sites, gene name, isoform names, true haplotype (if VCF phased), MCMC estimated haplotype, most often occuring haplotype in MCMC iterations, mean rho conditional on MCMC haplotype, posterior mode of rho, posterior median of rho, poseterior mean of rho, Lower 95% credible interval, Upper 95% credible interval, Lower 95% HPD interval, Upper 95% HPD interval, effective sample size for rho, run time, Gelman-Diag diagnostic, proportion rho < 0.5, proportion rho > 0.5, number of iterations, number of burnin, SNV locations, data type name, down-sample proportion, comma separated list of reads used from each input file. 9. Prepare isoform level data

julia ~/.julia/v0.5/IDPASE/src/prep_runs.jl -a ~/.julia/v0.5/IDPASE/test_data/SGS.psl ~/.julia/v0.5/IDPASE/test_data/TGS.psl -g ~/.julia/v0.5/IDPASE/test_data/TDRKH.gpd -v ~/.julia/v0.5/IDPASE/test_data/sim.vcf -q ~/.julia/v0.5/IDPASE/test_data/SGS.fq ~/.julia/v0.5/IDPASE/test_data/TGS.fq -d temp/ -c 1 -f 1 1 -o isoform_files/ -p sim -l 100 -i -s -e -r gene_out/gene_results.txt

where -l is read length for short reads, -s to skip file pre-processing (if using same GPD/VCF files are gene level), -e to use estimated haplotypes from gene leve, otherwise will use information from VCF, asumming it is phased, -r is the gene level results file. 10. Get commands to run each isoform individually

julia ~/.julia/v0.5/IDPASE/scripts/phase_isoforms_by_loci.jl -i isoform_files/ -o isoform_out/ -b ~/.julia/v0.5/IDPASE/scripts -a -p sim > to_run_isofs.sh

11. Run each isoform

bash to_run_isofs.sh  

12. concatenate all isoform level results

find isoform_out/ -name "EXTRA*" | xargs cat > isoform_out/isoform_results.txt

where the columns are: Gene name, gene level rho, comma separated isoform names, comma separated tau estimates, comma separated isoform abundances, which isoform each abundance correspond to, total reads, isoform lengths