User:Daly

Revision as of 10:37, 21 February 2006 by Daly (talk | contribs)

Overview

geWorkbench is an open-source bioinformatics platform that offers a comprehensive and extendible collection of tools for the management, analysis, visualization and annotation of biomedical data. This tool is aimed at providing researchers a centralized repository for the data analysis and visualiztion of gene expression data, gene sequences and protein sequences.

What can you do with geWorkbench?

  • Use one program that intergrates with multple existing bioinformatics modules for analysis and visualization.
  • Access remote servers and clusters for the performance of computationally intensive calculations - quicker!
  • Streamline data analysis with:
    • flexible import options that support merging files from various sources (insert possible datasource ie RMA express).
    • support for a variety of genomic data including microarrays, sequences, pathways, networks, alignments and phenotypes.
    • access analyses with biological annotations from the National Cancer Institute.
  • Community: Insert developer benefit ( plugin)


Slide1.gif

Tutorial

Welcome!

Welcome to the geWorkbench tutorial. In the following tutorials you will learn how to use geWorkbench. This an exmaple based,illustrated guide for both beginners and experienced users.


Using the Tutorial

These tutorials assume that you have successfully completed the installation instructions. If you haven't installed all the program, please go the Download and follow the installation instructions.

The tutorials are organized into a number of separate topics. Sample data sets that are used throughout the tutorials. Each tutorial is self contained and does not depend on any other tutorial for sample data.


(we should have sample data labeled Tutorial #)


  Within the tutorials there are two basis types of text:   
  Text that looks like this explains topics.
  1. Text in numbered steps which are instructions for you to follow.

==Getting Started== MV

  • Starting the application
  • GUI elements
    • Panels
    • Navigation

== Loading Data== MV

  • Data formats

Working with Marker and Phenotype Panels

In the this tutorial, you will:

  • Create active phenotype panels
  • Classify the panels you created


Before you can continue, geworkbench should be running. For help with this, please refer to the Getting Started (insertlink) section. The file cardiomyopathy.exp used in the this tutorial is from the Load Data (insertlink)tutorial.

Creating Panels

When working with microarrays, geWorkBench uses the term marker to refer to a gene probe (in other cases, it can be individual items from other data sets, such as sequences). phenotype refers to any user-defined grouping of microarrays. These microarrays will often share some common property that in most cases is phenotypic, although this is not a requirement. For example, one such “phenotype” might be a single experiment on a tumor tissue sample, with a second “phenotype” defined as a collection of experiments performed on normal tissue samples.


Assign panel

1. In the Phenotype panel, select and the following arrays in the contain samples from the congestive cardiomyopathy disease state. JB-ccmp_0120.txt, JB-ccmp_0218.txt, JB-ccmp0718.txt, JB-ccmp0811.txt, JB-ccmp1003.txt, JB-ccmp_1109.tx

2. Right-click, select Add to Panel.

3. Enter "Cardio" in the input box and click OK.


T PanelLabelCardio.png

4. Next, similarly label the follow arrays as "Normal" ( repeat steps 2 & 3 ). JB-n_0106.txt , JB-n_0821.txt, JB-n_0915.txt, JB-n_1303.txt


5. Select the checked boxes next to the panel name to indicate that these groups of arrays are "Active". Various analysis and visualization components can be set to only use/display activated arrays or markers.


T PhenotypesPriorToCase.png


Classify panel

For statistical tests such as the t-test, Case and Control groups can be specified.

  1. Left-click on the thumb-tack icon in front of the phenotype name.
  2. Selectcare to specifythe disease arrays as the "Case". The remaining "Normal" arrays are by default labeled control.

T PhenotypeSettingCase.png


A red thumbtack indicates the arrays have been specified as "Case".

T PhenotypeCaseSet.png

Visualize Gene Expression

In the this tutorial, you will:

  • Get acquiainted with the various visualization tools and view a dataset in geWorkbench
  • Modify the visualization preference settings


Before you can continue, geworkbench should be running. For help with this, please refer to the Getting Started (insertlink)section. The file cardiomyopathy.exp is used in the this tutorial is from the Load Data (insertlink)tutorial.


Visualization Tools

Visualization tools provide a view of the chip(s) under investigation and can be used for ascertaining the quality of the data. Active gene and phenotype panels (insertlink) restricts the number of markers/arrays displayed.The images created can be saved and exported. A detailed description on how to manipulate each visualization is described in online help.



Microarray View: Used to inspect each separate microarray using the scroll bar. Ema.gif
Tabular Microarray Panel: Presents the numerical values of the expression measurements in a table format. One row is created per individual marker/probe and one column per microarray. Etab.gif
Color Mosaic: Heat maps for microarray expression data, organized by phenotypic or gene groupings. ECm.gif
Expression Profiles: This is a line graph of genes expression profiles across several arrays/ hybridizations.. Eep.gif
Scatter Plot: A pairwise (array vs. array and marker vs. marker) comparison and plotting of expression values. Esp.gif



View a dataset

  1. Select a cardiomyopathy.exp in the Project Panel.
  2. Select the Microarray Panel visualization component in the View Area at the top-right section of the interface.
  3. Deaselect the All Markers checkbox to display the entire dataset.


Note: '''All Arrays''' and '''All Markers''' checkboxes determine which data points are included in the display. If neither is checked, then the entire data set is shown. The All Arrays control is useful when working with data sets comprising multiple arrays. In this case, only those arrays that are included in a currently activated phenotype panel will be displayed.

Allm.gif

Preferences

The Preferences selection in the Tools menu allows users to specify how certain aspects of the system will behave. Once your preferences are set, these preferences are persistent between application sessions and are applied at once.

Modifying settings

1. From the main menu, click on Tools>Preferences.

2. In the Preferences pop-up window, you can define settings for:

  • Text Editor: The editor selected will be used to open and inspect data sets loaded in a project. Notepad is the defeault setting.
  • Visualization: The color scheme selected be applied to color mosaic images.
    • Absolute: (default) Let M = max{|min|, |max|} over all expression measurements, across all arrays. If expression value x > 0, assign it the red spectrum x / M * 256. If expression value x is negative, assign it to the green spectrum -x / M * 256.
    • Relative: This is similar to the setting for Absolute, but each marker is mean-variance normalized first.
  • Genepix Value Computation: You can specify how compute the value displayed for Genepix array. The default setting is Option (Mean F635 - Mean B635) / (Mean F532 - Mean B532).


Select Relative for the visualization preference.


3. Click on OK.


Preferences.gif

Filter and Normalize Data

In the this tutorial, you will:

  • Get acquiainted with the various filters and normalizers
  • Apply a filter and normalizer


Before you can continue, geworkbench should be running. For help with installation, please refer to the Getting Started (INSERT LINK)section. The file cardiomyopathy.exp is used in the this tutorial is from the Load Data (insertlink)tutorial.



Filter

Filtering can be used to screen out missing data points, remove low quality data or reduce the size of the dataset by removing less interesting data.Available geWorkbench filters are as follows:


Filter Description
Affy Detection Call Applicable to Affymetrix data only. Sets all measurements whose detection status is any user-defined combination of A, P or M as missing.
Missing values Discards all markers that have “missing” measurements in at least n microarrays, where n is defined by the user. Another filter must first be applied however, in order to generate the missing values upon which this filter can operate.
Deviation Sets all markers whose measurements deviate below a given value across all microarrays as missing.
Expression Threshold Sets all markers whose measurements are inside (or outside) a user-defined range as missing.
2 Channel Applicable to 2-channel arrays (Genepix) data only. Defines applicable ranges for each channel, and sets all values for which either channel intensity is inside (or outside) the defined range as missing.


Filterpanel.gif


Perform the following steps to filter out data called absent in a Affymetrix file:

  1. In the Filtering Panel, select Affy Detection Call Filter.
  2. Select ‘A’ (Absent) checkbox and Filter.Values that were removed (marked as missing) are highlighted in yellow.
  3. In the Filtering Panel, select Missing Values Filter.
  4. Choose the maximum number of arrays that can have missing values before marker is removed – default is 0.
  5. Click Filter. Markers with more than 0 missing values are removed. You’ll notice the yellow values are gone


Affy Detection Call Filter Missing Values Filter
Filtered.gif Mvfilter.gif

Normalize

Normalization can be used to decrease the effects of systematic differences across a set of experiments. In geWorkbench, normalization results in replacing values with new values. Available geWorkbench normalization methods are as follows:


Normalizer Description
Missing value calculation Replaces every missing value with either the mean value of that marker across all microarrays or with the mean measurement of all markers in the microarray where the missing value is observed
Log2 Transformation Applies a log2 transformation to all measurements in a microarray
Threshold Normalizer All data points whose value is less than (or greater than) a user-specified minimum (maximum) value are raised (reduced) to that minimum (maximum) value
Marker-based Centering Subtracts the mean (median) measurement of a marker profile from every measurement in the profile
Array-based centering Subtracts the mean (median) measurement of a microarray from every measurement in that microarray
Mean-variance normalizer For every marker profile, the mean measurement of the entire profile is subtracted from each measurement in the profile and the resulting value is divided by the standard deviation

Normalpanel.gif

Apply Quantile Normalizer


1. In the Normalization Panel, select Quantile Normalizer.


2. Leave the default averaging method of Mean Profile Marker to indicate handling of missing values..


3. Click Normalize. The View Area is updated to reflect normalization (after the screen has been refreshed). Note: The first value in the second row was update from 41,394.6 to 55,779.26.

PRENORMALIZATION NORMALIZED
Prenormalizer ed.gif Postnormalizera.gif.

== Clustering Gene Expression Data== ken

  • Hierarchical Clustering
  • Self Organizing Map (SOM)

Differential Expression

In the this tutorial, you will:

  • Get acquiainted with the T-Test and Multi T-Test
  • Apply a Ttest and Multi T-Test


Before you can continue, geworkbench should be running. For help with installation, please refer to the Getting Started (INSERT LINK)section. The file cardiomyopathy.exp is used in the this tutorial is from the Load Data (insertlink)tutorial



T Test

In the T Test analysis, the data is divided into two groups, case and control. The t-test asks, for each marker, whether there is a significant difference between the two groups. To perform this analysis, you must classify the panels, set the analysis paramters and view the results in the visulization components. ( do I need to have why


Classify the Panels

1. Mark the Cardio phenotype a 'Case'. PBy default, panels are marked as control.

  • Right-click on Cardio phenotype.
  • Select Classification.
  • Select Case radio button. In the Phenotype panel, items marked Case are shown with a red thumbtack icon.

2. Activate the arrays Normal and Cardio by selecting the checkboxes next to the panel name.


T PhenotypeSettingCase.png

Set Analysis Parameters

  1. From the Analysis Panel, select T-Test Analysis.
  2. Populate the parameters values and click on Analyze.
  • In Alpha-corrections tab: Just Alpha.
  • In P-Value Parameters tab: p-values based on t-distribution. Note that the default alpha (critical p-value) is set to 0.01.
  • In Degree of Freedom tab: Welch approximation - unequal group variances.

Ttest.gif


T-Test Results

In the Gene Panel, the markers which met the significance test are included in a new active panel called “Significant Genes”. The values of the T-Test can be seen in the Color Mosaic panel and the Volcano Plot. (DESCRIBE BOTH)


Perform the following steps to to view the results: (TEST THIS)

  1. Uncheck All Markers so as to include only markers in the active panels.
  2. Uncheck All Arrays (though here it will have no effect as all arrays have been activated in the Phenotype Panel).
  3. Click on All to display the markers in Significant Genes Panel ( the active Gene Panel from the t-test). The p-values from the t-test are also shown.


VOLCANO PLOT COLOR MOSAIC
Vplot.png Ed cm.png
INSERT DESCRIPTION The gene height and width values can be altered to make the display more readible.



  • Multi Test
    • Volcano Plot
    • Color Mosaic

Regulatory Network

  • Reverse Engineering
  • Cytoscape

Integrated Annotation Information

In the this tutorial, you will:

  • rrrrr
  • rrrrr

Before you can continue, geworkbench should be running. For help with this, please refer to the Getting Started (insertlink) section. The file (INSERT FILE NAME) used in the this tutorial.


Perform the following steps to view annotations of a dataset(Not recommended for large sets of genes):


1. In the Marker Annotations panel in the View pane, click Use Panels.

A table of gene names and links to pathway diagrams is returned if available from the NCI CGAP database.


2. Clicking on gene names returns annotations from the NCI CGAP database in your default web browser window.


3. Clicking on the pathway diagram name will display a pathway diagram in the caBIO Pathways panel in the View pane. View pathway display in caBIO pathways viewer - pathway genes symbols can in turn be clicked to display their gene information in the default web browser window.

Enrichment Analysis

In the this tutorial, you will:

  • Go Term


Before you can continue, geworkbench should be running. For help with this, please refer to the Getting Started (insertlink) section. The file (INSERT FILE NAME) used in the this tutorial.

== Sequence Analysis== ken

  • Sequence Retrieval
  • Sequence Homology Analysis
    • Blast
    • Other

==Pattern Discovery== ken

  • Position Histogram

Promoter Analysis