Noi/NOTES/2012-4-26: Difference between revisions

From ZhangLabWiki
Jump to navigation Jump to search
>Noi
No edit summary
>Noi
No edit summary
Line 4: Line 4:
== mQTL analysis on UCLA-SZ data set (no CpG-SNP) ==
== mQTL analysis on UCLA-SZ data set (no CpG-SNP) ==
* Details will be filled later
* 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
* Note: P-value of 10%FDR is '''p <= 1.636E−4'''
* Note: P-value of 10%FDR is '''p <= 1.636E−4'''
{| {{table}}
{| {{table}}

Revision as of 21:15, 26 April 2012

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 





  • 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