User:Daly

Revision as of 15:07, 16 February 2006 by Daly (talk | contribs) (Normalize)

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.
    • 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: decribe this aspect
  • Insert developer benefit ( plugin)


Slide1.gif

Tutorial

Welcome to the geWorkbench tutorial.

In the following tutorials, you will learn how to use geWorkbench. While reading the tutorials, when you see numbered steps like the following, it's an instruction for you to follow.

  1. Click on Submit.
  2. Close the window.


Before You Begin

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

Getting Started

  • Starting the application
  • GUI elements
    • Panels
    • Navigation

Loading Data

  • Data formats

Working with Marker and Phenotype Panels

In the this tutorial, you will:

  • Create active phenotype panels
  • Classify the panels you created
  • Visualize the data

Before you can continue, geworkbench should be running. For help with this, please refer to the Getting Started (insertlink) section.The data we will use in the this tutorial is from the Load Data 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 on 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 the 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
  • Modify 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.



Microarray View: Used to inspect each separate microarray using the scroll bar. image]])
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. image]])
Color Mosaic: Heat maps for microarray expression data, organized by phenotypic or gene groupings. image]])
Expression Profiles: This is a line graph of genes expression profiles across several arrays/ hybridizations.. image]])
Scatter Plot A pairwise (array vs. array and marker vs. marker) comparison and plotting of expression values.

.||image]])

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.


Insert image

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.


(insert image)

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 filter methods:

  • Missing Values
  • Deviation
  • Expression Threshold
  • Affy Detection
  • 2 Channel Threshold


[[[[{|style="border: 2px solid darkGray" !|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.]]) |- |-


|- |} Perform the following steps to filter out data called absent:

  1. In the Filtering Panel, select the 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 on Filter. Markers with more than 0 missing values are removed.You’ll notice the yellow values are gone


   (insert image)'

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.]])

Apply Quantile Normalizer


1 In the Normalization Panel, select the Quantile Normalizer from the selection list.


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.

   (insert image)'

Clustering Gene Expression Data

  • Hierarchical Clustering
  • Self Organizing Map (SOM)

Differential Expression

T Test

Differential expression using t-test.

In this test, the data are divided into two groups, case and control. The t-test asks, for each marker, whether there is a significant difference between the two groups.

  1. Create a new project ( DESCRIBE)
  2. Right click on the project, select Open File.( DESCRIBE)
  3. Activate the arrays to be included in the t-test, highlight the two phenotypes and right click, then select activate. Only data in active panels will be used in differential expression.
  4. Mark the cardio phenotype as case. By default, panels are marked as control.

Perform the following steps to mark case:

  1. Right-click on Cardio phenotype.
  2. Select Classification.
  3. Select Case radio button. In the Phenotype panel, items marked Case are shown with a red thumbtack icon.
  4. From the Analysis Panel, select T-Test Analysis.
  5. Populate the parameters values.
  • 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.
  1. Click on Analyze.

T-Test Results

Switch to the Gene Panel.  The markers which met the significance test are included in a new panel called “Significant Genes”. 
    • Note that this new panel of genes is active.

The values of the T-Test can be seen in the Color Mosaic panel. Perform the following steps to modify the default view to view the results:

Uncheck All Markers so as to include only markers in the active panels. Uncheck All Arrays (though here it will have no effect as all arrays have been activated in the Phenotype Panel).

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.

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

Enrichment Analysis

  • Go Term
    • Go Miner

Sequence Analysis

  • Sequence Retrieval
  • Sequence Homology Analysis
    • Blast
    • Other

Pattern Discovery

  • Position Histogram

Promoter Analysis