Ns126:Calendar/NOTES/2016-1-5

From ZhangLabWiki
Revision as of 21:47, 19 January 2016 by >Shicheng (→‎Evaluation of the data quality)
Jump to navigation Jump to search

Collection RRBS Dataset From Richard Myers, HAIB (Encode Project)

Data Description

  • Normal tissues from RRBS data were downloaded from the ENCODE USCS genome browser [1]. Raw sequence data from 16 tissues were collected (Supplementary Table). Sequences were mapped to the bisulfite-converted human genome (NCBI hg19) using Bismark [2], after trimming the first three nucleotides, using the following parameters: --bowtie2 --phred64-quals --fastq -L 20 -N 1.
  • 1.RRBS and ENCODE, http://hgdownload.cse.ucsc.edu/goldenPath/hg19/encodeDCC/wgEncodeHaibMethylRrbs/.
  • 2.Krueger, F. and S.R. Andrews, Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics, 2011. 27(11): p. 1571-2.
  • Main important tissues were selected and merge to prediction model
ENCFF000LWX RRBS uterus adult female 44 year 36bp Illumina Genome Analyzer IIx
ENCFF000LXB RRBS uterus adult female 44 year 36bp Illumina Genome Analyzer IIx
ENCFF000LVO RRBS lung adult female 83 year 36bp Illumina Genome Analyzer IIx
ENCFF000LVR RRBS lung adult female 83 year 36bp Illumina Genome Analyzer IIx
ENCFF000LVA RRBS kidney adult female 83 year 36bp Illumina Genome Analyzer IIx
ENCFF000LVB RRBS kidney adult female 83 year 36bp Illumina Genome Analyzer IIx
ENCFF000LUT RRBS breast adult female 21 year 36bp Illumina Genome Analyzer IIx
ENCFF000LUV RRBS breast adult female 21 year 36bp Illumina Genome Analyzer IIx
ENCFF000LUQ RRBS brain adult male 66 year 36bp Illumina Genome Analyzer IIx
ENCFF000LUU RRBS brain adult male 66 year 36bp Illumina Genome Analyzer IIx
ENCFF000LWS RRBS stomach adult female 83 year 36bp Illumina Genome Analyzer IIx
ENCFF000LWW RRBS stomach adult female 83 year 36bp Illumina Genome Analyzer IIx
ENCFF000LVU RRBS pancreas adult male 71 year 36bp Illumina Genome Analyzer IIx
ENCFF000LVW RRBS pancreas adult male 71 year 36bp Illumina Genome Analyzer IIx
ENCFF000LWO RRBS zone of skin adult female 83 year 36bp Illumina Genome Analyzer IIx
ENCFF000LWP RRBS zone of skin adult female 83 year 36bp Illumina Genome Analyzer IIx
ENCFF000LVI RRBS mononuclear cell adult unknown 54 year 36bp Illumina Genome Analyzer IIx
ENCFF000LVK RRBS mononuclear cell adult unknown 54 year 36bp Illumina Genome Analyzer IIx
ENCFF000LVJ RRBS liver adult female 83 year 36bp Illumina Genome Analyzer IIx
ENCFF000LVN RRBS liver adult female 83 year 36bp Illumina Genome Analyzer IIx

Fastq Download

xargs -n 1 curl -O -L < haib.download.files.txt

Check Phred Score

perl ~/bin/checkphred.pl *fastq

Fastq Quality Control

  • trim_galore can used gz compressed fastq file as the input (so as the bismark)
# trim_galore --phred64 --fastqc --illumina --non_directional --rrbs *.fastq
trim_galore --phred64 --fastqc --illumina --rrbs *.fastq

Fastq Alignment

  • bismark can used gz compressed fastq file as the input (so as the trim_galore)
bismark --bowtie2 --phred64-quals --fastq -L 30 -N 1 /home/shg047/db/aligndb/hg19/bismark -1 ENCFF000MLM_trimmed.fq -o ../bam

Methylation haploinfo

bismark_methylation_extractor \
--single-end \
--bedGraph \
--buffer_size 2G \
--remove_spaces \
--zero_based \
--merge_non_CpG \
--comprehensive \
--output ../methyfreq \
ENCFF000MLM_trimmed.fq_bismark_bt2.bam
bismark_methylation_extractor --single-end --bedGraph --buffer_size 2G --remove_spaces --zero_based --merge_non_CpG --comprehensive --output ../methyfreq  ENCFF000MLM_trimmed.fq_bismark_bt2.bam

Methylfreq Data Output

  • Bismark Run Perl Beginning Script
#/usr/bin/perl
use strict;
use Cwd;
my $dir=getcwd;
chdir $dir;
my @file=glob("*.fastq");
foreach my $file(@file){
open OUT,">$file.bismark.sh";
print OUT "cd $dir\n";
chomp(my $phredcheck=`perl /home/shg047/bin/checkphred.pl $file`);
my ($phred)=split /\s+/,$phredcheck;
my $phred="--phred$phred";
print OUT "trim_galore $phred --fastqc --illumina --non_directional --rrbs $file\n";
print OUT "bismark --bowtie2 --phred64-quals --fastq -L 30 -N 1 /home/shg047/db/aligndb/hg19/bismark -1 $file\_qual_trimmed.fastq -o ../bam\n";
print OUT "bismark_methylation_extractor --single-end --bedGraph --buffer_size 2G --remove_spaces --zero_based --merge_non_CpG --comprehensive --  output ../me
}

Methylation Haplotype Load Matrix

  • methylation haplotype load for Haib data were save in TSCC
 /home/shg047/oasis/Haib/haploinfo/methHapLoad.matrix.txt

Evaluation of the data quality

  • cluster analysis

Merge with our RRBS data