Noi/NOTES/2012-4-26

From ZhangLabWiki
Revision as of 21:23, 26 April 2012 by >Noi
Jump to navigation Jump to search

mQTL analysis on UCLA-SZ data set (no CpG-SNP)

  • Details will be filled later
  • Filter out SNP from methylMatrix (Note: the sample IDs were modified to be compatible with tfam file and plink)
./remove_snps.pl ./UCLA-4batches_methyl_min76_minSTD0.1-edit /t4T/GenomeDB/HsGenome/snp134_snv.txt > UCLA-4batches_methyl-noSNP_min76_minSTD0.1

61664 UCLA-4batches_methyl-noSNP_min76_minSTD0.1

  • Split the methylMatrix in to 10,000 CpG sites each
split -10000 UCLA-4batches_methyl-noSNP_min76_minSTD0.1 UCLA_noSNP_
10000 UCLA_noSNP_aa
10001 UCLA_noSNP_ab
10001 UCLA_noSNP_ac
10001 UCLA_noSNP_ad
10001 UCLA_noSNP_ae
10001 UCLA_noSNP_af
1666 UCLA_noSNP_ag
  • Run mQTL (no CpG-SNP)
nohup ../mQTL_2MB_plink_allP.pl ../UCLA.hg19SNP ../UCLA_noSNP_aa > 2MB_noNSP_aa
nohup ../mQTL_2MB_plink_allP.pl ../UCLA.hg19SNP ../UCLA_noSNP_ab > 2MB_noNSP_ab 
nohup ../mQTL_2MB_plink_allP.pl ../UCLA.hg19SNP ../UCLA_noSNP_ac > 2MB_noNSP_ac
nohup ../mQTL_2MB_plink_allP.pl ../UCLA.hg19SNP ../UCLA_noSNP_ad > 2MB_noNSP_ad
nohup ../mQTL_2MB_plink_allP.pl ../UCLA.hg19SNP ../UCLA_noSNP_ae > 2MB_noNSP_ae
nohup ../mQTL_2MB_plink_allP.pl ../UCLA.hg19SNP ../UCLA_noSNP_af > 2MB_noNSP_af
nohup ../mQTL_2MB_plink_allP.pl ../UCLA.hg19SNP ../UCLA_noSNP_ag > 2MB_noNSP_ag
  • Concatenate the out put files to the same file
cat mQTL_noSNP_a*/2MB_noNSP_a* > 2MB_mQTL-noSNP_rawoutput
awk '{if ($1 ~ /chr/) print $0}' 2MB_mQTL-noSNP_rawoutput > 2MB_mQTL-noSNP_output
  • Number of test
433,947 2MB_mQTL-noSNP_output 
  • Sort out adj. p-val and count the number of significant associated CpG-SNP roughly at different p-value
sort column7 (adj p_val)
sort -g -k7 2MB_mQTL-noSNP_output | awk '{if ($7<=0.05) print $0}' | wc -l --> 6,418
sort -g -k7 2MB_mQTL-noSNP_output | awk '{if ($7<=0.01) print $0}' | wc -l --> 2,292
sort -g -k7 2MB_mQTL-noSNP_output | awk '{if ($7<=0.001) print $0}' | wc -l --> 1,026
sort -g -k7 2MB_mQTL-noSNP_output | awk '{if ($7<=0.0001) print $0}' | wc -l --> 698
sort -g -k7 2MB_mQTL-noSNP_output | awk '{if ($7<=0.00001) print $0}' | wc -l --> 526
  • Print out the site at min p-value 0.05 before multiple test correction
sort -g -k7 2MB_mQTL-noSNP_output | awk '{if ($7<=0.05) print $0}' > sorted_2MB_mQTL-noSNP_minPVal0.05.txt
6418 sorted_2MB_mQTL-noSNP_minPVal0.05.txt





  • Note: P-value of 10%FDR is p <= 1.636E−4
Associated CpG-SNP distance Number of assiciated CpG-SNP %
10%FDR_10kb 378 50.94%
10%FDR_10-20kb 27 3.64%
10%FDR_20-30kb 21 2.83%
10%FDR_30-40kb 36 4.85%
10%FDR_40-50kb 2 0.27%
10%FDR_50kb-2Mb 278 37.47%
Total in 2Mb 742


Associated CpG-SNP distance Number of assiciated CpG-SNP %
minpval0.05_10kb 862 13.43%
minpval0.05_10-20kb 124 1.93%
minpval0.05_20-30kb 133 2.07%
minpval0.05_30-40kb 108 1.68%
minpval0.05_40-50kb 75 1.17%
minpval0.05_50kb-2Mb 5,116 79.71%
Total in 2Mb 6,418