WashU EpiGenome Browser使用教程

WashU EpiGenome Browser 是我用过最赞的浏览器,没有之一。希望大伙跟着教程好好学习下!

还有更多教程见:http://epigenomegateway.wustl.edu/support/workshop2015keystone.html

在这个浏览器里面可以任意加载已有的公共数据,选择基因,选择坐标

这个教程总共讲解了12个操作范例,如下:

This tutorial has 12 demos:

[1-3] Loading the EpiGenome Browser and public data hubs

[4-5] Exploring Metadata

[6] Genome navigation

[7] Sessions

[8-9] Track visualization & customization

[10] Adding new tracks

[11] Using EpiGenome Browser apps

[12] Loading human genetic variation tracks

Loading the EpiGenome Browser and public data hubs

  1. To access the Browser, go to epigenomegateway.wustl.edu/browser. Browse the menu to select the genome you want to use.
  • For today, click “human hg19”.
  • Then click “Public hubs” to go access publicly available data.
  1. Load the public data hub.
  • A data hub is a collection of publicly available data. There are several available, including hubs of data from the Roadmap Epigenome Project and the ENCODE Consortium. The “Reference human epigenomes from Roadmap Epigenomics Consortium” is a collection of 4 hubs.
  • Click the “Reference human epigenomes from Roadmap Epigenomics Consortium” box to access these hubs.
  • Then load the “Roadmap Data from GEO” hub by clicking the “Load” button.
  • Click the “X” at the top of any floating window or press “Esc” to close it.
  1. Orientation to the Wash U Epigenome browser view.

Exploring Metadata

  1. Use the colormap panel and information box to find metadata terms associated with tracks.
  2. The information box for each track displays metadata terms associated with that track and information about the experiment that generated that track.
  • [A] Right-click a track to access the information box to learn more about the sample. [B] See how the metadata are organized hierarchically.

[A] [B]

  1. The colormap organizes metadata terms for quick display.
  • [A] Add metadata terms to the colormap by right-clicking on one of the titles, “Assay” or “Sample” and select “+ Add metadata terms”.
  • [B] Mark the “Track type” check box to add that category to the metadata colormap.
  • [C] Drag the RefSeq genes track above the chromosome ideogram. See that it now acquires a box in the metadata heatmap. Mouse-over to see that the “Track type” is “hammock”.

[A] [B] [C]

  1. The facet table is an intuitive way to browse large amounts of data by using metadata to stratify datasets.
  • The green numbers in each cell represent the number of tracks available for that Assay/Sample. The red number is the number of tracks for that Assay/Sample currently displayed.
  • Access the facet table through the Tracks button. Then click the track number icon.

Genome navigation

  1. Genome navigation in the Wash U EpiGenome browser.
  • Use the navigation box to relocate to EDC locus. Enter the locus coordinates in the “coordinate” search box: chr1:151880250-153605000.

Sessions

  1. The sessions app allows you to save a browsing status to revisit later. Each session can have several versions.
  • This workshop has a pre-saved session associated with it that includes different versions throughout the tutorial. If you get lost during the workshop, you can catch up by loading the pre-saved session. To start again at any point, load the saved session status for the section just completed, indicated by the ✰ and green highlighted-text.
  • [A] Click the Apps button to enter the Sessions function from the apps menu.
  • [B] Enter the session ID in the “Retrieve” text box and click “Retrieve”. The session ID for this workshop is ddSQFzBWvS.

[A] [B]

  • [C] Then select the 1) EDC locus version of this session.

[C]

Track visualization & customization

  1. Track navigation & customization is processed by right-clicking a track to enter the Configure menu.
  • [A] Enter the configure menu for H3K9me3 fibroblasts by right-clicking on the H3K9me3 track.
  • [B] Click the “positive” button to change the color to blue. To exit the menu, click anywhere outside the configure menu box.

[A] [B]

  • You can change the y-axis scale, track height, and other rendering features in the Configure menu.
  1. Use the matplot feature to compare two numerical tracks on the same y-axis scale.
  • [A] Use the multiple select function to select both H3K9me3 tracks. Hold the Shift button and right-click on both H3K9me3 tracks. They will be highlighted in yellow.
  • [B] Then right-click to open the options menu and choose Apply matplot.

[A] [B]

  • Using matplot, the two tracks can share the same y-axis. This makes the data easy to compare, in this case for the same histone modification ChIP-seq in two different cell types.
  • To catch up, open the saved session and the status: 2) H3K9me3 matplot.

Adding new tracks

  1. Use the dataset search box to find data sets using keywords.
  • To test our hypothesis that the EDC genes are cell type-specifically expressed, we will add keratinocyte and fibroblast RNA-seq datasets.
  • [A] Click the Tracks button to access the dataset search box. Search for “keratinocyte AND RNA”.
  • Note: the search function is case-sensitive and plural-sensitive!
  • [B] Then choose the result “RNA-seq of Penis Foreskin Keratinocyte Primary Cells”. It will turn green. Click the green “Add 1 track” button to add the track.

[A] [B]

  • Repeat the search for “fibroblast AND RNA”. Add the “RNA-seq of Penis Foreskin Fibroblast Primary Cells” track.
  • Both tracks appear at the bottom of the browser view.
  • To catch up, open the saved session and the status: 3) RNA-seq tracks.

Using EpiGenome Browser apps

  1. Use applications found in the browser’s Apps menu to generate quantitative tests of our hypothesis.
  2. The Gene set app allows the user to submit a list of genes or genomic positions to the browser for analysis.
  • To examine our hypothesis that the SPRR genes are differentially expressed between epidermal keratinocytes and dermal fibroblasts, we will use the SPRR gene set list.
  • [A] Open the Apps menu and navigate to the Gene & region set app.

[A]

  • The “SPRR genes” gene set has been pre-loaded in the session version for this workshop.
  • [B] Click the “SPRR genes” box and the “edit” button to view the gene set configurations.

[B]

  • In the genes panel, the double green arrow button can be used to find more gene models for genes in the gene set. This workshop uses all RefSeq gene models.
  • Below the genes panel is the gene-region specifications. This tells the browser what gene-associated regions to analyze. We want to analyze the entire gene plus the 5kb upstream region.

[C] Configure the gene set region:

  • (1) Click the “change >>” button to access the default gene-regions.
  • (2) Ensure that the radio button for “5’ and 3’ flanking” is selected.
  • (3) Move the green cursor to the “5 kb” mark.
  • (4) Then select “entire gene or interval” radio button.

[C]

  • [D] To trigger the gene set view, click on the banner box for “SPRR genes”. Then click the app “gene set view”.

[D]

  • The browser is now displaying the specified gene and gene-associated regions contiguously in the browser.
  • To catch up, open the saved session and the status: 4) Gene set view.
  1. To test the hypothesis that the DNase hypersensitivity regions are proximal enhancers, we will use the Scatterplot app to plot DNase vs H3K4me1 signal in the keratinocyte sample.
  • [A] Navigate to the Scatterplot app in the Apps menu.

[A]

  • [B] In the Scatterplot app window:
  • (1) Click “Choose a gene set” button and select the “SPRR genes”.
  • (2) Click “choose track>>” to select the x-axis track.
  • (3) Select “H3K4me1 of Penis Foreskin Keratinocyte Primary Cells”.
  • (4) Click “choose track>>” to select the y-axis track.
  • (5) Select “DNase hypersensitivity of Penis Foreskin Keratinocyte Primary Cells”.
  • (6) Click “SUBMIT” to run the app.

[B]

  • Now the app window is displaying the scatterplot with specified x- and y-axes. The result shows a general positive correlation, where an increase in DNase hypersensitivity corresponds with an increase in H3K4me1 signal in keratinocytes. Configure the scatterplot using the options at the below the graph.
  1. Last use the Gene plot app to quantify the RNA signal over SPRR genes in both cell types.
  • [A] Navigate to the Apps menu and use the app search box to search for the “Gene plot” app. Select “Gene plot” to enter the gene plot app window.

[A]

  • [B] Configure the gene plot data:
  • (1) In the first box of the app window, click the “Choose a gene set” button and select “SPRR genes”.
  • In the “Data track” box, click “Select numerical track >>” (2) and choose the “RNA-seq of Penis Foreskin Keratinocyte Primary Cells” (3).

[B]

  • [C] Configure the gene plot graph:
  • (4) In the “Graph type” box, select the “gene parts” plot.
  • (5) Click the orange button to “Make gene plot”.

[C]

  • [D] The gene plot show the data for the RNA-seq track chosen over the SPRR gene bodies. Data are averaged, binned, and plotted on a normalized “metagene” along the x-axis. The y-axis is the RNA-seq values.
  • [D] To generate the same plot for fibroblast RNA-seq data, click the “Go back” button.

[D]

  • [E] Then change the numerical track to “RNA-seq of Penis Foreskin Fibroblast Primary Cells” and click the “Make gene plot” button again.
  • [F] Now the gene plot is plotting data from the fibroblast RNA-seq track over the same set of genes. We can see they have much lower expression in this cell type.

[F]

Loading human genetic variation tracks

  1. View human genetic variation tracks in the EpiGenome Browser.
  • The EpiGenome Browser allows visualization of human genetic variation. HapMap and dbSNP tracks are found under the “Variation” category of Annotation tracks.
  • [A] Click the “X” to exit the gene set view to return to the linear genome.

[A]

  • [B] Click the Tracks button and select “Annotation tracks”.
  • [C] Scroll to “Population variation,” then chose “dbSNP release 137”.

[B] [C]

  • [D] The dbSNP track is displayed below the genes track. SNPs are color-coded according to their type of mutation, e.g. deletion or insertion.

[D]

大家可以先跟着英文教程试着摸索一下,可以在我们论坛上面分享自己的经验哈。

如果搞不定,我们会抽空录制一个视频给大家,敬请期待!

原文发布于微信公众号 - 生信技能树(biotrainee)

原文发表时间:2017-03-29

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏数据库新发现

VCS学习笔记

VERITAS Cluster Server(VCS) connects, or clusters, multiple, independent systems...

512
来自专栏Android先生

快给你的app上锁吧(android图案解锁)

思路 这里又是一个九宫格布局,布局可以参考上一篇快给你的app上锁吧(android数字解锁),只不过这里的九宫格上我们画的是图片(bitmap)。onDraw...

1232
来自专栏李智的专栏

python pandas教程

1182
来自专栏FreeBuf

Office”组合”式漏洞攻击样本分析

by hcl, nine8 of code audit labs of vulnhunt.com 1 概述 网上公开一个疑似CVE-2014-1761的RTF样...

2009
来自专栏函数式编程语言及工具

Akka(9): 分布式运算:Remoting-远程构建式

   上篇我们讨论了Akka-Remoting。我们说Akka-Remoting是一种点对点的通讯方式,能使两个不同JVM上Akka-ActorSystem上的...

1969
来自专栏开发与安全

muduo网络库学习之muduo_inspect 库涉及到的类

muduo inspect 库通过HTTP方式为服务器提供监控接口, 现在只实现进程相关信息的监控,通过成员ProcessInspector 实现。 Pro...

1835
来自专栏机器学习入门

LWC 63:748. Shortest Completing Word

LWC 63:748. Shortest Completing Word 传送门:748. Shortest Completing Word Problem: ...

1919
来自专栏iOSDevLog

Xcode 9 制作 iOS 通用 Framework创建 Framework模拟器和真机通用 FrameworkBundle使用 Framework

3397
来自专栏落影的专栏

使用AudioToolbox编码AAC

前言 使用VideoToolbox硬编码H.264 使用VideoToolbox硬解码H.264 这次在编码H.264视频流的同时,录制并编码AAC音频流。...

3637
来自专栏C/C++基础

计算机基础知识

攻击过程如下: 1. 攻击者通过MIM(比如arp欺骗等)劫持server与客户端浏览器之间的http包; 2. 攻击者生成一对伪造的RSA密钥: fa...

601

扫码关注云+社区