The recent development of various omics technologies in biology has increased the systems level analytical methods. It includes complete genome, transcriptome, proteome and metabolome analysis with various high-throughput technologies, including microarrays, Next-Generation Sequencing (NGS) and mass spectrometry. Multi-omics data integration provides an opportunity to understand the cellular or tissue level environment as a whole. IntOmix is an application specifically geared towards integrated network analysis of multi-omics with a focus on the identification of significantly changing metabolic network modules between two experimental conditions. The integration of the data generated by these two approaches can give simultaneous information about metabolomic changes along with the transcriptional changes.
Scope of the App
- Go from differentially expressed data to visualisations of pathway-level changes in less than 5 minutes.
- Process single or multi-omics data.
- Search for publications based on keywords added and significantly perturbed genes or metabolites by adding/editing keywords before or after optimization.
- Compare insights obtained from the data with existing literature.
- Compare results obtained from the same data at different parameters.
- Download results as publication-quality figures in various formats such as .png, .jpeg, .pdf or .xgmml.
IntOmix requires two types of .csv input files:
- Metabolite differential expression file
- Gene differential expression file
Both the input files consist of information about log2 fold change and the p-value of each gene or metabolite. In case of comparisons between multiple cohorts, you can also provide the name of the cohorts in the columns state1 and state2 as shown in Figure 2. You can also use files without the cohort information as shown in Figure 3.
Steps in data processing
- Upload differentially expressed metabolomics/transcriptomics/genomics files.
- Set parameters to perform the optimization of the network.
- Perform Pathway Enrichment analysis.
- Perform cross-cohort comparison to find common genes and metabolites across the various analysis.
- You can run the app using any single file as input if required.
- KEGG or HMDB ID for metabolites and gene id for genes is required to map these on the pathway map.
- If state 1 and state 2 columns are not present in the input files, the entire data will be treated as if from one comparison with the default cohort names, Cohort1 and Cohort2.
- In case there are blank cells for p-value or log2 fold change in the input files, the corresponding metabolites or genes are dropped before optimization and a notice is shown to that effect.
- Currently, we only support reactions from Homo sapiens, Mus musculus and Drosophila melanogaster.
- At a time only one keyword can be annotated as the primary keyword.
- Publication Search module requires at least one keyword to be annotated as the primary keyword.
- To utilize the Cross Cohort Comparison functionality, run your analysis at least twice at different parameters.
Select IntOmix from the dashboard under the Multi-omic Data Tab as shown in Figure 4. Create a New Workspace or choose an existing one from thr drop-down and provide a Name of the Session to be redirected to IntOmix's upload page.
The file upload interface allows you to upload input files to be processed in IntOmix as shown in Figure 5.
Click on Upload .csv Files to upload the differential expression files in the format mentioned above and then click on Continue to Dashboard.
This interface allows you to adjust parameters to be used to generate the optimized network.
- You can select cohorts as per your requirements.
- The organism will be detected automatically based on data input if the gene file has been uploaded.
- You can either proceed with the default parameters or specify the values. To know more about the significance of the parameters, refer to the Details of the App section.
- To isolate keywords from a paragraph, copy it on the left side of the Add Keywords section and click on Filter Keywords. The isolated keywords will be displayed below.
- To add specific keywords, write them on the right side of the Add Keywords section separated by a comma and click on Add Keywords.
- Click on any particular keyword to annotate it as the primary keyword. The primary keyword selection can be changed on the next interface as well as we can see in Figure 6.
- After specifying all the required parameters, click on Get Started to process your data.
- Processing of data requires as much time as specified to optimize the significantly perturbed networks present in the data.
After data processing, a KEGG map is generated based on the input parameters. The most significantly changing genes and metabolites between the two conditions are captured in the form of a network module. This network module is overlaid on the KEGG pathway to effectively visualize the overall topology of the variations between the two conditions. Hovering on any node or edge will display a dialog box with information about the selected gene or metabolite as shown in Figure 8. You can also zoom in, zoom out, select nodes or edges on the interactive map or download it in one of the several formats available.
You can add all metabolites and genes which do not satisfy the parameters specified for the optimized network.
Option and Information
In the Option & Information pop-up window as shown in Figure 9, you can select All Others option to add all genes and metabolites in the data to the generated map. The genes and metabolites have color ranging from red-green scale . All upregulated metabolites or genes are represented by a shade of red and downregulated metabolites or genes as a shade of green. The missing metabolites and genes are represented by shade of blue.
Analyze your data
After the KEGG network is created, the predefined set of queries can be selected from Analyze your Data option which gives you more insights on the results derived.
What are the key pathways?
This option enables you to find the key pathways that are perturbed in the data. The results of this query are based on two scoring methods; pathway enrichment and topology analysis. Pathway enrichment analysis refers to the quantitative approach based on the concentration of metabolites present. Gene set enrichment analysis is one of the famous algorithms used for enrichment analysis. In PollyTM IntOmix, R packages such as globaltest and GlobalAncova are used to analyze the pathway enrichment. Similarly, for topology analysis, betweenness centrality and degree centrality are considered as key parameters to find the most influential and important pathways in the whole subnetwork.
How do publicly available results correlate with my results?
This question enables you to select any gene or metabolite that is significantly differentially expressed in their data along with the specified primary keyword. Using both of these as the input, IntOmix will display a list of relevant publications that you can read more to gain insights about their data and generate hypotheses.
How do my other comparisons compare with my results?
Using this question you can compare the analysis of the same data using different parameters. The output of this question is displayed as a Venn diagram wherein the intersection contains the common genes and metabolites found in both the set of analyses. This helps you understand how a single parameter or a group of parameters affect their analysis.
Watch the demo video below for more information.
Details about the app
IntOmix is based on CoMBI-T i.e., Concordant Metabolomics Integration with Transcription pipeline to utilize the multiple omics techniques (metabolomics and transcriptomics) to focus on the identification of significantly altering metabolic network modules between two experimental conditions under investigation. This network includes both metabolites and enzymes. The most co-ordinately responsive sub-network is identified using the integrated network analysis and also, the novel points of metabolic rewiring are revealed.
Finding a module
After making a network, one can find a connected module that contains the most significantly changed genes and metabolites. Internally, this is done by first scoring nodes and edges based on their p-values in such a way that positive scores correspond to significant p-values and negative scores correspond to insignificant changes. Then the problem of finding a connected subgraph with maximum summary weight is solved.
The FDR values and Score for absent metabolites options control the size of the module. Increasing/decreasing the FDR for reactions (FDR for metabolites) value makes adding reactions (metabolites) to a module easier/harder.
It is recommended to start from the default values and then gradually change them depending on the results. Default FDR values are generated so that the module will be the size of ~100 reactions.
- Sergushichev, Alexey A., et al. "GAM: a web-service for integrated transcriptional and metabolic network analysis." Nucleic acids research 44.W1 (2016): W194-W200.
- Goeman, Jelle J., and Peter Bühlmann. "Analyzing gene expression data in terms of gene sets: methodological issues." Bioinformatics 23.8 (2007): 980-987.
- Hummel, Manuela, Reinhard Meister, and Ulrich Mansmann. "GlobalANCOVA: exploration and assessment of gene group effects." Bioinformatics 24.1 (2007): 78-85.
- Jha, Abhishek K., et al. "Network integration of parallel metabolic and transcriptional data reveals metabolic modules that regulate macrophage polarization." Immunity 42.3 (2015): 419-430.