worked_example_differential_example
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
worked_example_differential_example [2019/12/20 18:22] – krian | worked_example_differential_example [2021/01/28 21:23] (current) – krian | ||
---|---|---|---|
Line 1: | Line 1: | ||
====== Worked example Differential signaling ====== | ====== Worked example Differential signaling ====== | ||
- | {{ : | ||
- | **1-** Log in into HiPathia. For further information on this step visit [[logging_in|Logging in]]. | ||
- | **2-** We will work with a Breast Cancer dataset from the repository The Cancer Genome Atlas. | + | ===== The oestrogen receptor (ER+) positive group Vs negative group (ER-) ===== |
+ | ==== Differential signaling inputs ==== | ||
+ | **Note:** You can also launch | ||
- | | + | **1-** Log into HiPathia. For further information on this step visit [[logging_in|logging in]]. |
- | | + | |
- | **3-** Upload the data to HiPathia in the data panel by clicking on //My data//. For further information on this step visit [[upload_your_data|Upload your data]]. | + | **2-** Collection of data. We will work with a Breast Cancer dataset from the repository The Cancer Genome Atlas (TCGA) [[https://portal.gdc.cancer.gov/projects/TCGA-BRCA |
+ | More information on the proposed dataset is available here: | ||
+ | * https:// | ||
+ | * https:// | ||
- | **4-** Press the //Differential signaling// button. | + | Before use in HiPathia, |
- | {{ :hipathia_bar_diffsig_5.png? | + | We have randomly selected 178 samples of breast cancer tumors from the dataset, annotated as ER+ and ER-(the molecular annotations come from [[https://www.ncbi.nlm.nih.gov/ |
- | **5-** In the //Input data panel//. Press the //File browser// of the // | + | The expression |
- | {{ :hipathia_work1.png? | + | * Expression matrix: [[http:// |
+ | * Experimental design: | ||
- | **6-** In the //Design | + | **3-** Upload |
- | {{ : | + | **4-** Click // |
- | **7-** Select Human as species (by default). | + | {{ : |
- | **8-** In the //Function level analysis//, check both checkboxes, | + | **5-** In the //Input data// panel, click the //File browser// button of the section |
- | **9-** In the // | + | {{ : |
- | **10-** In the //Study information// | + | As mentioned before, the number of samples is 356 (178 ER+ and 178 ER-), and this information appears after selecting |
- | {{ :hipathia_work3_5.png?600 |}} | + | {{ :hipathia_work1-2_6.png?600 |}} |
- | **11-** | + | **6-** In the //Design data// panel, select //Two group comparison// |
+ | |||
+ | {{ : | ||
+ | |||
+ | Automatically, | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | **7-** Select //Human (Homo sapiens)// as species (default). | ||
+ | {{ :: | ||
+ | |||
+ | **8-** In //Effector annotation source//, check both the //Gene ontology// and //Uniprot keywords// boxes. | ||
+ | |||
+ | **9-** In the // | ||
+ | |||
+ | **10-** In the //Study information// | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Give a name to the study, for example, " | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | **11-** | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | A study will be created and listed in the studies panel. You can access this panel by clicking the //My studies// button. | ||
+ | {{ : | ||
+ | |||
+ | This a video tutorial for differential signaling study: {{: | ||
+ | |||
+ | ==== Differential signaling results and interpretations ==== | ||
+ | Once the analysis is finished at HiPathia web, the report/ | ||
+ | |||
+ | For further information about the differential signaling report please visit [[differential_signaling# | ||
+ | /* | ||
+ | Furthermore, | ||
+ | {{ :: | ||
+ | The downloaded RAR file has to be extracted then you can open the //" | ||
+ | |||
+ | The results are divided in different panels: | ||
+ | === Study Information === | ||
+ | There is more information about the study available here: | ||
+ | {{ :: | ||
+ | * **Name**: the study name. | ||
+ | * **Description**: | ||
+ | * **Tool**: the name of the used tool (in this case, is Hipathia). | ||
+ | * **Date**: the launching date of the study (MM/ | ||
+ | === Input parameters === | ||
+ | Here you can visualize the parameters with which the current study was launched. | ||
+ | {{ :: | ||
+ | * **Expression file**: The name of the expression file that has been used in the current study. | ||
+ | * **Design file**: The name of the design file that has been used in the current study. | ||
+ | * **Comparison**: | ||
+ | * **Paired analysis**: Have the input data been paired? **No** or **Yes**. | ||
+ | * **Species**: | ||
+ | === Pathways === | ||
+ | {{ : | ||
+ | A pathway viewer has been developed to facilitate the visualization of the comparison results. It shows a graphical representation of the activation levels of the " | ||
+ | {{ :: | ||
+ | The information summarized in the pathway viewer includes: | ||
+ | |||
+ | Node shape: | ||
+ | * **Ellipses** correspond to genes. | ||
+ | * **Circles** represent metabolites. | ||
+ | * **Rectangles** are cellular functions. | ||
+ | |||
+ | Node color: | ||
+ | * **Blue** for down-regulated genes, **red** for up-regulated genes and **grey** for non-significantly differentially expressed genes. | ||
+ | * The intensity of the color depends on the significance of the differential expression. | ||
+ | |||
+ | Edge shape: Edges represent either activations or inhibitions. | ||
+ | * **Activations** are represented by classic arrows. | ||
+ | * **Inhibitions** are represented by T arrows. | ||
+ | |||
+ | Edge color: | ||
+ | * Circuits which are **up-regulated** in the ER- condition have their arrows colored in **red**. | ||
+ | * Circuits which are **down-regulated** in the ER- condition (condition 1 ) have their arrows colored in **blue**. | ||
+ | * Circuits which are **not significant** for the test have their arrows colored in **grey**. | ||
+ | * **When an edge belongs to different circuits which should be in different colors, different arrows are depicted, each one for each necessary color.** | ||
+ | The pathway viewer has four principal parts in order to facilitate the visualization of the results: | ||
+ | {{ : | ||
+ | In the upper-right part of the visualization tool, all selected pathways from the differential signaling form are shown along with one or two arrows. These arrows indicate whether in one of the " | ||
+ | ---- | ||
+ | {{ : | ||
+ | In the lower-right part of the tool will appear all the circuits in which a pathway (previously selected on the upper part -1-) can be decomposed.\\ Once any of the circuit is clicked the nodes and interactions (edges) that form part of this circuit are highlighted in the pathway viewer (), One example might be the red-highlighted circuit on the figure below. | ||
+ | {{ :: | ||
+ | ---- | ||
+ | In the visualization part there are two types of objects, the nodes and the edges. | ||
+ | {{ : | ||
+ | * The nodes represent the different proteins or metabolites that are responsible on passing the signal from the previous node to the next one. As commented before, the nodes can be plain or complexes. In the former it is only necessary the presence of one protein (although that same function can be done by several proteins), while in the complex ones it is necessary the presence of two or more proteins to pass the signal through.\\ The color on the nodes are the result of a differential expression analysis on the gene/s involved on performing that node's function. | ||
+ | * The edges represent how the interactions between the different nodes are.\\ If the edge is an arrow then the previous node will be activating the next one, while if it ends with a vertical bar is the former node will inhibit the functionality of the following node. This interactions may be depicted in red or blue depending on the circuit they form part of (whether they are up-regulated or down-regulated). It may occur, that two or even three colors for a same edge are shown, but this only happens when a circuit is not yet selected on the lower-right part of the tool.\\ Once an certain circuit is selected all its edges will be colored in the same color depending on the result of the differential signaling activation analysis.\\ | ||
+ | ---- | ||
+ | The top part contains the title of selected pathway,by clicking this button {{:: | ||
+ | {{ : | ||
+ | You can also find other buttons like: | ||
+ | * {{:: | ||
+ | {{ :: | ||
+ | * {{:: | ||
+ | * {{:: | ||
+ | * {{:: | ||
+ | * {{:: | ||
+ | === Circuit values === | ||
+ | Here you can find additional tables and plots: | ||
+ | * The first section below contains a link to download the computed circuit activity values. This matrix file indicates for each effector circuit the level of activation calculated using the HiPathia method for each sample.{{ :: | ||
+ | * The Heatmap plot represented on the hipathia results page is a heatmap of the activation values from the most differentiated effector circuits (rows) between groups along with a clustering of the samples (columns). This plot allows to observe if its possible to differentiate the groups that are compared according to its effector circuit activation values.\\ The colors depicted here indicates the level of activation for the different circuit in each sample, being the bluish ones those with lower activation levels and the reddish the highest ones (e.g. if a given cell of the heatmap is blue would indicate that this particular effector circuit is poorly activated in that determined sample) {{ :: | ||
+ | * The next plot on the report corresponds to a Principal Component Analysis plot. This figure is useful to determine if the activation levels of the pathways calculated with Hipathia are able to differentiate between the two groups that are being compared. {{ :: | ||
+ | * In the section below there is a table (that can be downloaded) with the results of the differential activation analysis. This table indicates for each effector circuit whether or not there is a different level of activation depending on the group the samples belong. {{ :: | ||
+ | * **circuit/ | ||
+ | * **UP/ | ||
+ | * **statistic**: | ||
+ | * **p.value**: | ||
+ | * **FDRp.value**: | ||
+ | * **Fold Change**: is a measure describing how much the ** circuit activity** changes between conditions.[[https:// | ||
+ | * **logFC**: The logarithm to base 2 of calculated Fold change. | ||
+ | === Function based analysis === | ||
+ | |||
+ | For each of the functional analysis selected to be performed, the same descriptive statistics as in the Path values case are shown: | ||
+ | |||
+ | * **GO terms/ | ||
+ | * **Heatmap**: | ||
+ | * **PCA**: A Principal Components Analysis is performed in order to see if the matrix contains information enough to separate between the two groups. | ||
+ | * **GO terms/ | ||
+ | * **path/ | ||
+ | * **UP/ | ||
+ | * **statistic**: | ||
+ | * **p.value**: | ||
+ | * **FDRp.value**: | ||
+ | * **Fold Change**: is a measure describing how much the ** circuit activity** changes between conditions.[[https:// | ||
+ | * **logFC**: The logarithm to base 2 of calculated Fold change. | ||
+ | */ | ||
+ | |||
+ | /* | ||
+ | ==== Results interpretations ==== | ||
+ | */ | ||
+ | |||
+ | **Estrogen receptor (ER) pathway (hsa04915): | ||
+ | |||
+ | If we look at the estrogen receptor pathway, we can see it is mainly downregulated in ER-negative (ER-) tumors. Considering that ER- tumors are characterized by lacking estrogen receptors, these results are consistent with prior knowledge that the estrogen pathway is downregulated in tumors that are not sensitive to estrogen [1]. | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Despite gene expression being upregulated in ER- for CREB3, HiPathia shows its activity to be decreased according to the signaling cascade that activates it. The role of CREB3 in cancer progression and metastasis is controversial, | ||
+ | |||
+ | |||
+ | **Focal adhesion pathway (hsa04510): | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | |||
+ | In the focal adhesion pathway, the activity of proteins involved in actin polymerization, | ||
+ | |||
+ | |||
+ | ==== References ==== | ||
+ | |||
+ | [1] Frasor J, Danes JM, Komm B, Chang KC, Lyttle CR, Katzenellenbogen BS. Profiling of estrogen up- and down-regulated gene expression in human breast cancer cells: insights into gene networks and pathways underlying estrogenic control of proliferation and cell phenotype. Endocrinology. 2003 Oct; | ||
+ | |||
+ | [2] Howley, B. V., Link, L. A., Grelet, S., El-Sabban, M., & Howe, P. H. (2018). A CREB3-regulated ER-Golgi trafficking signature promotes metastatic progression in breast cancer. Oncogene, 37(10), 1308–1325. DOI: 10.1038/ | ||
+ | |||
+ | [3] Kim, HC., Choi, KC., Choi, HK. et al. HDAC3 selectively represses CREB3-mediated transcription and migration of metastatic breast cancer cells. Cell. Mol. Life Sci. 67, 3499–3510 (2010). DOI: 10.1007/ | ||
+ | |||
+ | [4] Padilla-Rodriguez, | ||
+ |
worked_example_differential_example.1576866172.txt.gz · Last modified: 2019/12/20 18:22 by krian