专栏首页生信技能树WashU EpiGenome Browser使用教程

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),作者:生信菜鸟团

原文出处及转载信息见文内详细说明,如有侵权,请联系 yunjia_community@tencent.com 删除。

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

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

我来说两句

0 条评论
登录 后参与评论

相关文章

  • Agilent芯片表达矩阵处理(学徒作业)

    Agilent的芯片同样也是扫描得到图片,然后图像处理(主要是Agilent Feature Extraction (AFE) 软件)得到信号值,但是值得注意的...

    生信技能树
  • Variant 分析阶段小结1-基础碎碎念

    所谓遗传变异是生物体内遗传物质发生变化而造成的可以遗传给后代的变异,这些变异导致了生物在不同水品上体现出遗传的多样性。生物信息学中各种基因组研究的基础就是遗传变...

    生信技能树
  • MATH值代表的肿瘤异质性在乳腺癌与生存关系不显著

    但是今天要分享的文章Breast Cancer Research and Treatment February 2017 , 题目是:Clinical and ...

    生信技能树
  • POJ-1953 World Cup Noise(线性动规)

    World Cup Noise Time Limit: 1000MS Memory Limit: 30000K Total Submissio...

    ShenduCC
  • CodeForces 665B Shopping

    B. Shopping time limit per test 1 second memory limit per test 256 megabyt...

    ShenduCC
  • SAP OData性能分析工具

    As mentioned by title, this blog does not introduce the OData trace functionalit...

    Jerry Wang
  • CC1310空中升级笔记03 CC1350_BLE_OAD_Project0

    如果打不开这个TI的链接,我上传了整个网页,http://download.csdn.net/detail/iotisan/9750301。

    twowinter
  • 什么是SAP物料主数据里的Batch

    Materials are produced and theoretically have the same properties. Nevertheless ...

    Jerry Wang
  • 如何启用SAP Business by design里的Correction Invoice功能

    Subject: [Tip] How to enable the function Correction Invoice for customer invoic...

    Jerry Wang
  • 循环热管的动态状态空间建模与基于模型的控制设计(CS SY)

    对于航空航天、汽车或服务器系统中电子元件的热控制,散热器通常远离热源。因此,热传导系统是有效冷却电子元件所必需的。循环热管(LHPs)就是这样的传热系统,它利用...

    用户6853689

扫码关注云+社区

领取腾讯云代金券