Chernobyl Thyroid Cancer - Methylation

I. Description

This workflow was used for the detection of differentially methylated cpg islands in the Chernobyl Thyroid Cancer Study.

Major steps in the workflow are:

1. Reading in normalized beta value data and sample table and preparing related data for the following analysis
2. Detecting differentially methylated cpg islands using the tool minfi
3. Detecting differentially methylated cpg islands using the tool ChAMP
4. Reporting the overlap of differntially methylated cpg islands detected by the two tools

II. Dependencies

• Snakemake
• R:
  • rmakrdown
  • data.table
  • minfi
  • ChAMP
  • ggfortify
  • ggplot2
  • pheatmap
  • plotly

III. Input

• config.yaml: select the type of phenotype data
• data/normBeta.txt: normalized beta value data stored in the directory: data/
• data/pheno.csv: sample table stored in the direcotry: data/
• data/MethylationEPIC_v-1-0_B2_anno.csv: probe annotation file stored in the directory: data/

IV. Output

• Results from minfi in directory: minfi/
• Results from ChAMP in directory: champ/
• Final results in directory: overlap/

V. Working directory structure

.
├── champ                                               
│   ├── champ_data.Rdata
│   ├── champ.html
│   ├── phenotype_champ_dmp_q0.01.csv
│   ├── phenotype_champ_dmp_q0.05.csv
│   ├── phenotype_champ_GSEA_path_sig_dmp_q0.01.csv
│   ├── phenotype_champ_GSEA_path_sig_dmr_q0.05.csv
│   └── phenotype_champ_sig_dmr_q0.05.csv
├── config.yaml 
├── data 
│   ├── data_read_in.html
│   ├── inital_data.Rdata
│   ├── MethylationEPIC_v-1-0_B2_anno.csv 
│   ├── normBeta.txt 
│   └── pheno.csv 
├── log 
│   └── log files
├── minfi 
│   ├── minfi_cat.html
│   ├── minfi_data.Rdata
│   └── phenotype_minfi_dmp_q0.01_annot.csv
├── overlap 
│   ├── check_overlap.html
│   ├── overlap_phenotype_dmp_q0.01.csv
│   ├── overlap_phenotype_dmp_q0.01fc2.csv
│   ├── overlap_phenotype_sig_dmp_q0.01fc2_beta_data.csv
│   └── overlap.Rdata
├── R_code 
│   ├── champ.Rmd
│   ├── check_overlap.Rmd
│   ├── data_read_in.Rmd
│   ├── minfi_cat.Rmd
│   └── minfi_con.Rmd
└── Snakefile