Hi, I have got exome seq data from Ion Torrent platform. I am looking to analyze these data for GWAS studies. Please suggest pipelines and workflow.
•
written
by
Kumar
•
written
by
Kumar
Hey Manoj,
You could produce VCFs (or BCFs) and then import these to PLINK, where you could then perform standard GWAS analyses. There are 2 related answers, here: Convert vcf phased data to plink
If you have Ion Torrent data, I see no reason why you could not use GATK ( a related question, with no answer: NGS preprocessing pipleine for ion torrent data )
Kevin
Hey Manoj, well, it really depends on what you want to do. Usually, people are studying some disease and they want to find inherited genetic variants that are statistically significantly associated with the disease. A typical GWAS is actually just a basic chi-squared test on a 2 x 2 contingency table, comprising the 'mutant' allele and reference allele tallies in the cases and controls. I go through the basic calculations for this here: A: SNP dataset and Z Score (technically, you don't need PLINK).
If you are studying a condition like arthritis, heart disease, or something along these lines, it's imperative that you have a solid study design and that your study is Powered (statistically) so that you can make inferences to the wider population.
Hi, my apologies, you have DNA-seq data and you want to identify germline variants in this, correct?:
Okay, then, yes, you can follow the GATK
Neither. Use bwa mem if your reads are >=70bp in length. If shorter, use bowtie2.
Yes, the difference relates to the data-type: RNA- versus DNA-seq
Hi, I used PLINK for the vcfs, it generates .bed,.bim, .fam with the following command. Please let me know if I am doing correct.
./plink --vcf <filename> --keep-allele-order --allow-extra-chr 0 --make-bed --out <filename>
Next I used king with following command:
./king -b <filename.bed> --ibs <out-filename>
However, it is showing following ERROR: FATAL ERROR -
Too many first alleles as the major allele (~35.5%). Please use plink1.9 --make-bed to regenerate the genotype data again.
Furthermore, I tried following command:
./plink --bfile <filenames> --recode --tab --out <output filename>
Hi, I have 90 VCF files. I am looking to use fastSTRUCTURE, please give some idea how to prepare input file for fastSTRUCTURE. I am new in this analysis. Initially, I got some idea that these VCF files need to be merged in one VCF file using BCFtools and then PLINK should be used to generate .bed file using the merged file. Please let me know if it is correct way.
Log in to answer this question.
