Differences

This shows you the differences between two versions of the page.

--- workvafin [2019/11/18 12:23] – krian
+++ workvafin [2021/01/19 17:37] (current) – [Results and interpretation] jdopazo
@@ Line 1: / Line 1: @@
 ====== Worked example Variant Interpreter ======
-===== Example 1: Fanconi Anemia study =====
+===== Diabetes study =====
 **1-** Log in into HiPathia. For further information on this step visit [[logging_in|Logging in]].
-**2-** We will work with expression data from custom tissue related to Fanconi Anemia (bone marrow derived cells). You can download the raw expression matrix used for this example from GEO repository, accession number [[https://www.ncbi.nlm.nih.gov/pubmed/19850743|GSE16334]]. You can read the original study here: Vanderwerf et al. TLR8-dependent TNF-(alpha) overexpression in Fanconi anemia group C cells. Blood. 2009 Dec 17;114(26):5290-8).
+**2-** We will work with expression data from custom tissue related to Diabetes. In this case we will work with a gene expression dataset of human islets from 63 donors (54 non-diabetic and 9 type 2 diabetic). You can download the raw expression matrix used for this example from GEO repository, accession number [[https://www.ncbi.nlm.nih.gov/pubmed/22768844|GSE38642]] You can read the original study here: https://pubmed.ncbi.nlm.nih.gov/22768844/
-Hipathia only takes as input normalized data, so you will have to correctly normalize the raw expression matrix to be uploaded.
-  * Custom tissue expression matrix after normalization: {{:onlyCexpreset_hinorm.txt|FA expression matrix}}
-**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]].
+For the purpose of this example we will only use the data of the 54 non diabetic donors. The Variant interpreter tool allows to explore the functional consequences of the variation in the expression of a gene or a combination of genes. The aim of this example is to investigate the functions and biological processes punitively/in silico affected by Loss-of-function (LoF) variants in the RAP1GAP gene. To do so, we will simulate the knock out of the gene RAP1GAP (reducing its expression value to 0), this gene is known to play a key role in inflammatory processes which are highly related to the development of type 2 diabetes mellitus.
-{{ :variantexample1.png |}}
-**4-** We will evaluate the impact of disrupting variants in gene //BRCA1// (Biallelic BRCA1 mutations are regarded either embryonically lethal or to cause Fanconi anemia (FA), a genomic instability syndrome characterized by bone marrow failure, developmental abnormalities, and cancer predisposition [[https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6732317/ | cita ]] ). We need to upload the list of genes in txt format the same way we did to upload expression data. Both gene names and entrez ids are accepted. Select Gene list as data type.
-**5-** Press the //Variant Interpreter// button.
-{{ :variantinterbtn.png |}}
-**6-** In the //Input data panel//. Press the //File browser// and select the desired file for gene list and custom tissue expression file (OPTIONAL).
-{{ :uploadgenelist.png |}}
-**7-** (OPTIONAL) Select one or more in the GTEx tissues list in order to compare between them.
-**8-** (OPTIONAL) Select a subset of signaling pathways to evaluate, all pathways are selected by default.
-{{ :pathwayslist.png |}}
-**9-** In the //Study information// panel, press the //File browser// button and select the desired output folder. In this case we will use //FA_KO_analysis//. Give a name to the study, for example "FA KO analysis".
-{{ :studyinfovafin.png |}}
-**10-** Press the //Run analysis// button. A study will be created and listed in the studies panel. You can access this panel by clicking on the //My studies// button.
-===== Example 2: Diabetes study =====
-**1-** Log in into HiPathia. For further information on this step visit [[logging_in|Logging in]].
-**2-** We will work with expression data from custom tissue related to Diabetes. You can download the raw expression matrix used for this example from GEO repository, accession number [[https://www.ncbi.nlm.nih.gov/pubmed/22768844|GSE38642]] You can read the original study here: https://www.ncbi.nlm.nih.gov/pubmed/22768844
 Hipathia only takes as input **normalized data**, so you will have to correctly normalize the raw expression matrix to be uploaded.
-  * Custom tissue expression matrix after normalization: {{:onlyCexpresetDiabetes.txt|Diabetes expression matrix}}
+  * Custom tissue expression matrix after normalization: [[http://hipathia.babelomics.org/data/exp_data_onlyControls.tsv|exp_data_onlyControls.tsv]]
 **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]].
-**4-** We will evaluate the impact of disrupting variants in each gene **independently** of these gene lists. Each gene list correspond to all the genes that compose a single circuit from related and non related pathways to the disease, circuits  are named after their effector node:
+**4-** We will evaluate the impact of disrupting variants of RAP1GAP. In this example the gene list will contain just the RAP1GAP gene
-  * {{ :listagenesKO_diabetesRap1.txt | Gene list for example 1: Rap1 signaling pathway: PRKCI PARD6A PARD3}}
+  * {{ :RAP1GAP.txt | Gene list for example 1: Rap1 signaling pathway: RAP1GAP}}
-  * {{ :listagenesKO_diabetesChemokine.txt | Gene list for example 2: Chemokynes signaling pathway: PARD3-PRKCZ-TIAM1}}
-  * {{ :listagenesKO_diabetesTNF.txt | Gene list for example 3: TNF signaling pathway:CREB3}}
+Then we will need to upload it the same way we did to upload expression data. Both gene names and entrez ids are accepted. Select Gene list as data type.
-After choosing one or more from the studied example and prepare a list of gene(s) in txt format as described [[gene_list_file_format|here]]. We need to upload it the same way we did to upload expression data. Both gene names and entrez ids are accepted. Select Gene list as data type.
 **5-** Press the //Variant Interpreter// button.
@@ Line 61: / Line 31: @@
 **9-** Press the //Run analysis// button. A study will be created and listed in the studies panel. You can access this panel by clicking on the //My studies// button.
+===== Results and interpretation =====
+When the analysis is finished at Hipathia web, the report/results will be available in “My studies” section once the study is done.\\  You can download any table or image showed in the report page by clicking on the name right before it.\\ You can also download the circuit activity values and function matrices by clicking on //Circuit values//, //GO terms values// and //Uniprot keywords values// respectively.
+As we can see in the figure below, the results after the simulating the LoF variant of gene RAP1GAP over the 54 non diabetic samples only show deregulation on the pathways where the gene RAP1GAP can be found. Since RAP1GAP is an inhibitor of the Rap1 pathway, mostly all the functions associated to the circuits that compose the pathway are UP-regulated. The functional consequences of the KO of RAP1GAP vary from UP-regulation of Inflammatory processes, Cell adhesion and Stress response among others. The results shown seem to indicate that RAP1GAP plays an important role in the control of these processes highly related to diabetes. Therefore, LoF mutations in those genes inactivating RAP1 protein would most likely result in an over-activation of inflammatory response and thus, an increased risk of developing immune-mediated diabetes ([[https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6908734/|Peña-Chilet et al., 2019]]).
+{{ ::rap1_signaling_pathway_hsa04015_1_.svg.png?nolink |}}
+Moreover, LoF variants vary considerably in their effects on human phenotype. They can represent severely deleterious disease-causing mutations, however, they will also include mildly deleterious variants with small effects on fitness, and neutral variants disrupting the function of non-essential genes. The Variant interpreter option provide a user-friendly tool for //in-silico// simulation of the consequences of LoF variants on transcriptomic data. The potential of the tool relies on its capability to foreground biological processes potentially affected by LoF variants, enabling to explore the functional map displayed by a single, or a set of, LoF mutations present in the condition studied. Subsequently, LoF variants affecting highly relevant functions could be interpret as a first line of evidence supporting the pathogenicity of that variant.