This workflow, developed in R markdown files as well as and Jupyter notebook files, includes differential gene expression analysis, statistical analysis of metabolomics data, as well as pathway enrichment analysis for both transcriptomics and metabolomics data followed by integration of this data through network analysis to identify disease-related processes and visualization of multi-omics data. A publicly available (https://ibdmdb.org/) gut-transcriptomic and stool-metabolome dataset of the gut microbial ecosystem in inflammatory bowel diseases was used to test the proposed workflow.

Each analysis section is given with corresponding script and the topic.

To perform the workflow:


Transcriptomics Analysis

TutorialTopic
1-Data preprocessing Preprocess transcriptomics data to be ready for analysis
2-Differential gene expression analysisPerforming differential gene expression analysis
3-Pathway AnalysisPerform pathway enrichment analysis to differential expressed genes
4-Heatmap creationHeatmap visualization for enriched pathways
5-Overlapped genes extractionExtracting overlapped differential expressed genes between Chron`s disease and ulcerative colitis
6-Network analysisCreating PPI-pathway network, applying MCL (Markov Clusetring) algorithm to the created networks

Metabolomics Analysis

TutorialTopic
7-Data preprocessing Preprocess metabolomics data to be ready for analysis
8-Significantly changed metabolites analysisPerforming statistical analysis
9-Metabolite pathway AnalysisPerforming pathway enrichment analysis to metabolite data

Multi-Omics visualization

TutorialTopic
10-Identifier mapping Identifier mapping to be ready for visualization
11-VisualizationEnriched pathway visualization