[Pub] Takashima-2014 Track Settings
 
Resetting Transcription Factor Control Circuitry Towards Ground State Pluripotency In Human

Maximum display mode:       Reset to defaults   
Select dataType (Help):
partially methylated domain       hypomethylated regions       methylation level ▾       coverage ▾       allele specific methylated regions      
Select subtracks by celltype:
Celltype
Human H9 primed R2
Human H9 primed R3
Human H9 primed R1
Human H9 R3
Human H9 R2
Human H9 R1
Human H9 primed
Human H9
List subtracks: only selected/visible    all    ()
  Celltype↓1 dataType↓2   Track Name↓3  
hide
 Human H9 primed R2  hypomethylated regions  Human_H9-primed_R2_HMR   Data format 
hide
 Configure
 Human H9 primed R2  methylation level  Human_H9-primed_R2_Meth   Data format 
hide
 Human H9 primed R3  hypomethylated regions  Human_H9-primed_R3_HMR   Data format 
hide
 Configure
 Human H9 primed R3  methylation level  Human_H9-primed_R3_Meth   Data format 
hide
 Human H9 primed R1  hypomethylated regions  Human_H9-primed_R1_HMR   Data format 
hide
 Configure
 Human H9 primed R1  methylation level  Human_H9-primed_R1_Meth   Data format 
hide
 Configure
 Human H9 R3  methylation level  Human_H9_R3_Meth   Data format 
hide
 Configure
 Human H9 R2  methylation level  Human_H9_R2_Meth   Data format 
hide
 Configure
 Human H9 R1  methylation level  Human_H9_R1_Meth   Data format 
    
Assembly: Human Feb. 2009 (GRCh37/hg19)

Takashima_Human_2014
We are planning to introduce the new version of methylone track hubs sometime between February 7 and February 14 2024. The following assemblies will be updated: mm39, gorGor6, canFam6, GCF_000001735.3, rn7, panTro6, hg38.

Description

Sample BS rate* Methylation Coverage %CpGs #HMR #AMR #PMD
Human_H9-primed_R2 Primed H9 ESC Rep 2 0.991 0.852 3.409 0.881 33794 235 3658 non-WGBS; LowCov; Download
Human_H9-primed_R3 Primed H9 ESC Rep 3 0.991 0.849 3.343 0.879 33478 192 3659 non-WGBS; LowCov; Download
Human_H9-primed_R1 Primed H9 ESC Rep 1 0.990 0.849 3.567 0.889 33831 240 3732 non-WGBS; LowCov; Download
Human_H9_R3 Naive H9 ESC Rep 3 0.995 0.309 2.553 0.848 0 79 4397 non-WGBS; LowCov; Download
Human_H9_R2 Naive H9 ESC Rep 2 0.995 0.307 2.916 0.871 0 96 4722 non-WGBS; LowCov; Download
Human_H9_R1 Naive H9 ESC Rep 1 0.996 0.301 2.853 0.864 0 75 4287 non-WGBS; LowCov; Download
Human_H9-primed Primed H9 ESC 0.991 0.850 10.319 0.949 41713 0 0 non-WGBS; Download
Human_H9 Naive H9 ESC 0.995 0.305 8.322 0.948 55264 0 0 non-WGBS; Download

* see Methods section for how the bisulfite conversion rate is calculated
Sample flag:
non-WGBS:  sample is not generated with whole-genome bisulfite sequencing (WGBS);
LowCov:  sample has low mean coverage (<6.0)

Terms of use: If you use this resource, please cite us! The Smith Lab at USC has developed and is owner of all analyses and associated browser tracks from the MethBase database (e.g. tracks displayed in the "DNA Methylation" trackhub on the UCSC Genome Browser). Any derivative work or use of the MethBase resource that appears in published literature must cite the most recent publication associated with Methbase (see "References" below). Users who wish to copy the contents of MethBase in bulk into a publicly available resource must additionally have explicit permission from the Smith Lab to do so. We hope the MethBase resource can help you!

Display Conventions and Configuration

The various types of tracks associated with a methylome follow the display conventions below. Green intervals represent partially methylated region; Blue intervals represent hypo-methylated regions; Yellow bars represent methylation levels; Black bars represent depth of coverage; Purple intervals represent allele-specific methylated regions; Purple bars represent allele-specific methylation score; and red intervals represent hyper-methylated regions.

Methods

All analysis was done using a bisulfite sequnecing data analysis pipeline MethPipe developed in the Smith lab at USC.

Mapping bisulfite treated reads: Bisulfite treated reads are mapped to the genomes with the rmapbs program (rmapbs-pe for pair-end reads), one of the wildcard based mappers. Input reads are filtered by their quality, and adapter sequences in the 3' end of reads are trimmed. Uniquely mapped reads with mismatches below given threshold are kept. For pair-end reads, if the two mates overlap, the overlapping part of the mate with lower quality is clipped. After mapping, we use the program duplicate-remover to randomly select one from multiple reads mapped exactly to the same location.

Estimating methylation levels: After reads are mapped and filtered, the methcounts program is used to obtain read coverage and estimate methylation levels at individual cytosine sites. We count the number of methylated reads (containing C's) and the number of unmethylated reads (containing T's) at each cytosine site. The methylation level of that cytosine is estimated with the ratio of methylated to total reads covering that cytosine. For cytosines within the symmetric CpG sequence context, reads from the both strands are used to give a single estimate.

Estimating bisulfite conversion rate: Bisulfite conversion rate is estimated with the bsrate program by computing the fraction of successfully converted reads (read out as Ts) among all reads mapped to presumably unmethylated cytosine sites, for example, spike-in lambda DNA, chroloplast DNA or non-CpG cytosines in mammalian genomes.

Identifying hypo-methylated regions: In most mammalian cells, the majority of the genome has high methylation, and regions of low methylation are typically more interesting. These are called hypo-methylated regions (HMR). To identify the HMRs, we use the hmr which implements a hidden Markov model (HMM) approach taking into account both coverage and methylation level information.

Identifying hyper-methylated regions: Hyper-methylated regions (HyperMR) are of interest in plant methylomes, invertebrate methylomes and other methylomes showing "mosaic methylation" pattern. We identify HyperMRs with the hmr_plant program for those samples showing "mosaic methylation" pattern.

Identifying partially methylated domains: Partially methylated domains are large genomic regions showing partial methylation observed in immortalized cell lines and cancerous cells. The pmd program is used to identify PMDs.

Identifying allele-specific methylated regions: Allele-Specific methylated regions refers to regions where the parental allele is differentially methylated compared to the maternal allele. The program allelicmeth is used to allele-specific methylation score can be computed for each CpG site by testing the linkage between methylation status of adjacent reads, and the program amrfinder is used to identify regions with allele-specific methylation.

For more detailed description of the methods of each step, please refer to the reference by Song et al. For instructions on how to use MethPipe, you may obtain the MethPipe Manual.

Credits

The raw data were produced by Takashima Y et al. The data analysis were performed by members of the Smith lab.

Contact: Benjamin Decato and Liz Ji

Terms of Use

If you use this resource, please cite us! The Smith Lab at USC has developed and is owner of all analyses and associated browser tracks from the MethBase database (e.g. tracks displayed in the "DNA Methylation" trackhub on the UCSC Genome Browser). Any derivative work or use of the MethBase resource that appears in published literature must cite the most recent publication associated with Methbase (see "References" below). Users who wish to copy the contents of MethBase in bulk into a publicly available resource must additionally have explicit permission from the Smith Lab to do so. We hope the MethBase resource can help you!

References

MethPipe and MethBase

Song Q, Decato B, Hong E, Zhou M, Fang F, Qu J, Garvin T, Kessler M, Zhou J, Smith AD (2013) A reference methylome database and analysis pipeline to facilitate integrative and comparative epigenomics. PLOS ONE 8(12): e81148

Data sources

Takashima Y, Guo G, Loos R, Nichols J, Ficz G, Krueger F, Oxley D, Santos F, Clarke J, Mansfield W, et al Resetting transcription factor control circuitry toward ground-state pluripotency in human. Cell. 2014 158(6):1254-69