#acl All:read
{{attachment:Software/EnrichmentMap/UserManual/enrichmentmap_logo3.png|Enrichment Map Logo|align="right"}}<<BR>> 
= Enrichment Map GREAT Tutorial =
<<TableOfContents(3)>>

== Outline ==

This quick tutorial will guide you through the generation of an Enrichment Map for an analysis performed using [[http://great.stanford.edu/great/public/html/| Genomic Region Enrichment Annotation Tool (GREAT)]],  

To run this tutorial: 
    * You need to have Cytoscape installed : minimally 3.1.0 must be installed but preferable to have the latest version of Cytoscape  (e.g. version 3.1.1)
    * Install the Enrichment Map App from the Cytoscape App manager in cytoscape 3. 
    * You need to download the test data: [[attachment:GREATTutorial.zip]]

Description of the tutorial files contained in the GREATTutorial folder:
    * TestRegions_ForGREAT.bed : Example GREAT binary genomic region input file.
    * GreatExprotAll.tsv : Example of download GREAT output file.
    * 20140919-public-2.0.2-3vD5MB-hg19-all-gene.txt : Example downloaded GREAT gene to region association file.


== Instructions ==
=== Step 1: Generate GREAT output files ===
  1. GO to GREAT website - [[http://great.stanford.edu/great/public/html/]]
  1. Select ''Species Assembly'' associated with your data.  For this tutorial select ''Human: GRCh37''
  1. In ''Test regions'' click on ''Choose File''
  1. Navigate to files provided and select ''TestRegions_forGreat.bed''
  1. Click on ''Submit''
  1. Once the results page has loaded download all the results - in the ''Global controls'' heading click on the down arrow next to ''Global Export''
  1.  Select ''All data as tsv'' - greatExportAll.tsv will automatically be downloaded to your default Downloads directory. --> This is the file you can use in Enrichment Map (Dataset 1 or 2:Enrichment Results)

{Optional} - Download the Gene-to-region used by GREAT and modify it to be used in EM as an expression file.
  1. In the ''Global controls'' heading click on the down arrow next to ''Global Export''
  1. Select ''view all region-gene associations''
  1. Next to ''Gene->genomic region association table'' [The table on the right hand side of the page] click on ''Download table as text''.
  1. File will automatically downloaded into your default Downloads directory (file name is similar to DATE-public-2.0.2-3vD5MB-hg19-all-gene.txt where DATE is the date of download.  Name will also change depending on the version of GREAT and genome selected).


=== Stpe 2: Generate Enrichment Map with GREAT Output ===
{{attachment:Great_inputpanel.png|Screenshot Great Input Panel|align="right"}}
   1. Open Cytoscape
   1. Click on Apps / Enrichment Maps / Load Enrichment Results
   1. Make sure the Analysis Type is set to DAVID/BiNGO/GREAT
   1. Please select the following files by clicking on the respective (...) button and selecting the file in the Dialog:
      * '''NO GMT file is required for GREAT Analysis'''
      * Dataset 1 / Expression: `GeneToRegion.txt` (OPTIONAL)
      * Dataset 1 / Enrichments: `GreatExportAll.tsv`     
   1. Tune Parameters
      * P-value cut-off `0.001`
      * Q-value cut-off `0.05`
      * Check Overlap Coefficient
         * Overlap coefficient cut-off `0.6`
   1. Build Enrichment Map
   1. Once the network starts to build a dialog will pop up asking you how you would like to filter the GREAT results.  There are four options:
      1. Hypergeometric test p-values and FDR.
      1. Binomial test p-values and FDR
      1. Both
      1. Either
   1. Select ''Both''
 
=== Step 3: Examining Results ===
{{attachment:GREAT_EM.png|GREAT EM Result}}<<BR>>
'''Legend''':
  1. Node (inner circle) size corresponds to the number of genes in dataset 1 within the geneset
  1. Node border (outer circle) size corresponds to the number of genes in dataset 2 within the geneset
  1. Colour of the node (inner circle) and border(outer circle) corresponds to the significance of the geneset for dataset 1 and dataset 2, respectively.
  1. Edge size corresponds to the number of genes that overlap between the two connected genesets.  Green edges correspond to both datasets when it is the only colour edge.  When there are two different edge colours, green corresponds to dataset 1 and blue corresponds to dataset 2.
    * '''NOTE''': if you are using two enrichment sets you will see two different colours of edges in the enrichment map.  When the set of genes in the two datasets are different (for example, when you are comparing two different species or when you are comparing results from two different platforms) the overlaps are computed for each dataset separately as there is a different set of genes that the enrichments were calculated on.  In this case, since the enrichments were reduced to only a subset of most differentially expressed at each time point the set of genes the enrichments are calculated on are different and overlap are calculated for each set separately.