Is it possible to do linkage analysis using whole exome data only (i.e. without using micro array data)? I would like to use whole exome sequence from small pedigrees (i.e children and parents) and be able to perform linkage analysis. I am mainly interested in defining haploblocks, homozygosity mapping, linkage regions detection etc. Would PLINK be an option for this ?
We're also thinking of using exome data to map loci - the idea seems reasonable:
- sequence exomes
- map both existing and new SNPs and corresponding genotypes across all samples
- construct custom genetic map file of markers from (2)
- format genotype data as a matrix similar to that obtained from microarrays (e.g PED file)
- run standard linkage / homozygosity mapping tool (Allegro, GeneHunter, etc)
PLINK does find ROHs (regions of homozygosity) but doesn't seem to do standard linkage. This paper describes how PLINK does "population-based linkage" based on analysis of shared IBS segments.
I have found this article and website interesting... http://bioinf.wehi.edu.au/software/linkdatagen/#mps http://genomebiology.com/2011/12/9/R85/
Will try it out and compare with our snpArray linkage in 4 patients (2 families)