This repo contains a series of functions and scripts to facilitate the execution of limma::roast (rotation gene set tests) as part of an MSstats workflow and generate tabular outputs and visualizations.
- Run roast against four gene set databases (GO, Reactome, KEGG and MSigDB),
- Select specific subcategories within each database (i.e.
MF
ontology for GO, Hallmark subcategory for MSigDB), - Setup a p-value cutoff for the tabular output,
- Select the size of the gene sets to be evaluated,
- Filter-out redundant (Semantically similar) GO terms or non-interesting ones (Blacklist) from Reactome or KEGG.
-
Download this repo into your computer. You should have two scripts and two folders.
-
Initialize the repo in your local computer as an RStudio project.
-
Add the required input data files into the
Data
folder:proteinGroups.txt
file from the MaxQuant output.evidence.txt
file from the MaxQuant output.annotation.tsv
file with the MSstats/MSstatsTMT experimental design (You can consult their documentation here: https://msstats.org/).
-
In RStudio, open
prepMSstats_format_script_vx.x.R
-
Modify the parameters in lines 14 (type of MSstats used) and 19 (the names of the two experimental conditions to compare as they appear in the
annotation.tsv
file. -
Click
Source
. -
The script will run MSstats/MSstatsTMT for summarization and normalization and will produce three new files:
annotation_to_ROAST.tsv
: a reduced version of your initial annotation file with only the selected conditions for comparison.ProcessedMSstatsData.rds
: an R object file containing the MSstats/MSstatsTMT processed data.ProcessedMSstatsData_to_ROAST.rds
: an R object file containing the MSstats/MSstatsTMT processed data, reduced to only contain information related to the two experimental conditions you want to compare.
-
Open the main script
enrichRoast_script_vx.x.R
. -
Modify the parameters between lines 6 and 100, by answering the questions in the comments.
-
Click
Source
. -
A new
Outputs
folder should have been generated with figures, tabular outputs and an HTML report.
- You will need a
txt
file namedinput_limma.txt
. The first column should be nameID
and each other column should correspond to different samples/MS runs. - Download this repo into your computer. You should have two scripts and two folders.
- Initialize the repo in your local computer as an RStudio project.
- Add
input_limma.txt
into theData
folder. - Open the main script
enrichRoast_script_limma_input_vx.x.R
- Modify the parameters between lines 8 and 122, by answering the questions in the comments.
- Click
Source
. - A new
Outputs
folder should have been generated with figures, tabular outputs and an HTML report.
The visualizeKEGG.R
script would help you map your enriched KEGG pathways to the quantified proteins in your data set. This script should only be executed after the ROAST enrichment analysis (after executing any of the previous enrichRoast_**
scripts).
- Open the
visualizeKEGG.R
script. - Execute the lines between 6 and 21 to check for the proper installation and loading of the required packages.
- Modify the input in lines 25 and 31 according to the questions in the comments.
- Execute the rest of the script.
- A number of
png
files with the protein to KEGG pathway mapping should have been created in your working directory.
- The
R
folder contains all the functions required to run this analysis in a modular way. The idea is to create a package from them so people in the lab could install it and implement these functions into their own scripts/workflows as they require it. Everyone is encouraged to test the functions individually too. I will soon write individual documentation for each one. - The processing and summarization with MSstats/MSstatsTMT takes time, this is why the script is generating an R object (
ProcessedMSstatsData.rds
/ProcessedMSstatsData_to_ROAST.rds
) to avoid running MSstats/MSstatsTMT every time we need to run roast. If you have already processed data through MSstats/MSstatsTMT you can save it using the functionsaveRDS()
with the nameProcessedMSstatsData_to_ROAST.rds
if it already contains only two experimental conditions or with the nameProcessedMSstatsData.rds
if you would require to use theprepMSstats_format_script_vx.x.R
for selecting the desired two conditions to compare. - If you have a
ProcessedMSstatsData_to_ROAST.rds
/ProcessedMSstatsData.rds
file, you can include it in theData
folder as an input to avoid running MSstats/MSstatsTMT every time you want to do a roast enrichment using these scripts/functions.