- SearchUI.exe is unclear why you need it, well, that is, I can’t say exactly what it does. However, we know that it is connected with the work of voice helper Cortana (the game is still that) and it also seems to be like a graphical search engine.
- Hi John, It's interesting that the permissions didn't seem to be maintained following the upgrade. However, I did a clean install then attached my machine to our corporate domain and Cortana/searchui.exe stopped working.
In this tutorial, we perform proteogenomic database searching using the Mass Spectrometry data. The inputs for performing the proteogenomic database searching are the peaklist MGF files and the FASTA database file. The FASTA database is obtained by running the first workflow “Uniprot_cRAP_SAV_indel_translatedbed.FASTA”. The second workflow focuses on performing database search of the peak list files (MGFs).
The SearchGui program is a subset of the HTMLtools Java program. For details, see the full HTMLtools reference manual ReferenceManual.html.
Agenda
In this tutorial, we will deal with:
- Match peptide sequences
- Create a SQLite database for peptide, protein and genomic annotation visualization
hands_on Hands-on: data upload and organization
- Create a new history and name it something meaningful (e.g. Proteogenomics DB search)
- Import the four MGF MS/MS files and the Trimmed_ref_5000_uniprot_cRAP.FASTA sequence file from Zenodo.
Tip: Importing via links
- Copy the link location
- Open the Galaxy Upload Manager (galaxy-upload on the top-right of the tool panel)
- Select Paste/Fetch Data
- Paste the link into the text field
- Press Start
- Close the windowBy default, Galaxy uses the URL as the name, so rename the files with a more useful name.
- Rename the datasets to something more recognizable (strip the URL prefix)
Build a Dataset list for the three MGF files, name it as “Mo_Tai_MGFs”
Tip: Creating a dataset collection
- Click on Operations on multiple datasets (check box icon) at the top of the history panel
- Check all the datasets in your history you would like to include
- Click For all selected.. and choose Build dataset list
- Enter a name for your collection
- Click Create List to build your collection
- Click on the checkmark icon at the top of your history again
The search database labelled Uniprot_cRAP_SAV_indel_translatedbed.FASTA
is the input database thatwill be used to match MS/MS to peptide sequences via a sequence database search.
For this, the sequence database-searching program called SearchGUI will be used.The generated dataset collection of the three MGF files in the history is used as the MS/MS input. We will walk through a number of these settings in order to utilize SearchGUI on these example MGF files.
SearchGUI
hands_on Hands-on: SearchGUI
- SearchGUItool with the following parameters:
param-file“Protein Database”:
Uniprot_cRAP_SAV_indel_translatedbed.FASTA
(Or however you named theFASTA
file)comment Comment:
The “Uniprot_cRAP_SAV_indel_translatedbed” FASTA database is obtained when you run the first proteogenomics workflow. Please make sure to run the 1st workflow.
param-files“Input Peak lists (mgf)”:
MGF files
dataset collection.tip Tip: Select dataset collections as input
Click the Dataset collection icon on the left of the input field:
Select the appropriate dataset collection from the list
- Section Search Engine Options:
param-check“DB-Search Engines”:
X!Tandem
comment Comment:
The section Search Engine Options contains a selection of sequence database searching programs that are available in SearchGUI. Any combination of these programs can be used for generating PSMs from MS/MS data. For the purpose of this tutorial, X!Tandem we will be used.
- Section Protein Digestion Options:
- param-select“Enzyme”:
Trypsin
- param-text“Maximum Missed Cleavages”:
2
- param-select“Enzyme”:
- Section Precursor options:
- param-select“Precursor Ion Tolerance Units”:
Parts per million (ppm)
- param-text“Precursor Ion Tolerance”:
10
- param-text“Fragment Tolerance (Daltons)”:
0.05
(this is high resolution MS/MS data) - param-text“Minimum charge”:
2
- param-text“Maximum charge”:
6
- param-select“Forward Ion”:
b
- param-select“Reverse Ion”:
y
- param-text“Minimum Precursor Isotope” :
0
- param-text“Maximum Precursor Isotope” :
1
- param-select“Precursor Ion Tolerance Units”:
- Section Protein Modification Options:
- param-select“Fixed Modifications”:
Carbamidomethylation of C, ITRAQ-4Plex of K, ITRAQ-4Plex of Ntermini
param-select“Variable modifications”:
Oxidation of M, ITRAQ-4Plex of Y
tip Tip: Search for options
For selection lists, typing the first few letters in the window will filter theavailable options.
- param-select“Fixed Modifications”:
- Section Advanced Options:
- param-select“X!Tandem Options”:
Advanced
- param-check“X!Tandem: Quick Acetyl”:
No
- param-check“X!Tandem: Quick Pyrolidone”:
No
- param-check“X!Tandem: Protein stP Bias”:
No
- param-text“X!Tandem: Maximum Valid Expectation Value”:
100
- param-check“X!Tandem: Quick Acetyl”:
- param-select“X!Tandem Options”:
Once the database search is completed, the SearchGUI tool will output a file (called aSearchGUI archive file) that will serve as an input for the next section, PeptideShaker.Rename the output as “Compressed SearchGUI results”
comment Comment:
Note that sequence databases used for proteogenomics are usually much larger thanthe excerpt used in this tutorial. When using large databases, the peptide identificationstep can take much more time for computation. In proteogenomics, choosing the optimaldatabase is a crucial step of your workflow.
PeptideShaker
PeptideShaker is a post-processing software tool thatprocesses data from the SearchGUI software tool. It serves to organize the Peptide-SpectralMatches (PSMs) generated from SearchGUI processing and is contained in the SearchGUI archive.It provides an assessment of confidence of the data, inferring proteins identified from thematched peptide sequences and generates outputs that can be visualized by users to interpretresults. PeptideShaker has been wrapped in Galaxy to work in combination with SearchGUIoutputs.
comment Comment: File Formats
There are a number of choices for different data files that can be generated using PeptideShaker. A compressed file can be made containing all information needed to view them results in the standalone PeptideShaker viewer. A mzidentML
file can be generated that contains all peptide sequence matching information and can be utilized by compatible downstream software. Other outputs are focused on the inferred proteins identified from the PSMs, as well as phosphorylation reports, relevant if a phosphoproteomics experiment has been undertaken.
hands_on Hands-on: PeptideShaker
- PeptideShakertool with the following parameters:
- param-file“Compressed SearchGUI results”: The SearchGUI archive file
- param-select“Specify Advanced PeptideShaker Processing Options”:
Default Processing Options
- param-select“Specify Advanced Filtering Options”:
Default Filtering Options
- param-check“Include the protein sequences in mzIdentML”:
No
- param-check“Output options”: Select the
PSM Report
(Peptide-Spectral Match) and theCertificate of Analysis
comment Comment: Certificate of Analysis
The “Certificate of Analysis” provides details on all the parametersused by both SearchGUI and PeptideShaker in the analysis. This can be downloaded from theGalaxy instance to your local computer in a text file if desired.
- Inspect galaxy-eye the resulting files
A number of new items will appear in your History, each corresponding to the outputs selected in the PeptideShaker parameters. The Peptide Shaker’s PSM report is used as an input for the BlastP analysis. Before performing BlastP analysis. The Query Tabular tool and few test manipulation tools are used to remove spectra that belongs to the reference proteins. The output tabular file “Peptides_for_Blast-P_analysis” will contain only those spectra that did not belong to any known proteins.
The mzidentml output from the Peptide shaker is converted into an sqlite database file by using the mz to sqlite tool. This sqlite output is used to open the Multi-omics visualization platform, wherein you can view the spectra of the peptides using Lorikeet parameters. To open the MVP viewer, click on the “Visualize in MVP Application” icon ( this will pop-open the interactive multi-omics viewer in a new window/tab)
hands_on Hands-on: mz to sqlite
This tool extracts mzidentml and its associated proteomics datasets into a sqlite db
- mz to sqlitetool with the following parameters:
- param-file“Proteomics identification files”:
PeptideShaker_mzidentml
- param-file“Proteomics Spectrum files”:
Mo_Tai_MGFs
- param-file“Proteomics Search Database Fasta”:
Uniprot_cRAP_SAV_indel_translatedbed.FASTA
- param-file“Proteomics identification files”:
The next step is to remove known peptides from the list of PSM’s that we acquired from the Peptide Shaker results. For that we need to perform Query tabular to extract list of known peptides from the UniProt and cRAP database.
Query Tabular
hands_on Hands-on: Remove Reference proteins
- Query Tabulartool with the following parameters:
- param-repeatInsert Database Table (b):
PSM report
- Section Filter Dataset Input:
- param-repeatInsert Filter Tabular Input Lines
- “Filter by”:
skip leading lines
- “Skip lines”:
1
- “Filter by”:
- param-repeatInsert Filter Tabular Input Lines
- Section Table Options:
- “Specify Name for Table”:
psms
- “Use first line as column names” :
No
- “Specify Column Names (comma-separated list)”:
id,Proteins,Sequence
- “Only load the columns you have named into database”:
Yes
- param-repeatTable Index
- “Table Index”:
No
- “Index on Columns”:
id
- “Table Index”:
- “Specify Name for Table”:
- Section Filter Dataset Input:
param-repeatInsert Database Table (c):
PSM report
- Section Filter Dataset Input
- param-repeatInsert Filter Tabular Input Lines
- “Filter by”:
skip leading lines
- “Skip lines”:
1
- “Filter by”:
- param-repeatInsert Filter Tabular Input Lines
- “Filter by”:
select columns
- “Enter column numbers to keep”:
1,2
- “Filter by”:
- param-repeatInsert Filter Tabular Input Lines
- “Filter by”:
normalize list columns,replicate rows for each item in the list
- “Enter column numbers to normalize”:
2
- “List item delimiter in column”:
,
- “Filter by”:
- param-repeatInsert Filter Tabular Input Lines
- Section Table Options:
- “Specify Name for Table”:
prots
- “Use first line as column names” :
No
- “Specify Column Names (comma-separated list)”:
id,prot
- “Only load the columns you have named into database”:
Yes
- param-repeatInsert Table Index:
- “Table Index”:
No
- “Index on Columns”:
prot,id
- “Table Index”:
- “Specify Name for Table”:
- Section Filter Dataset Input
- param-repeatInsert Database Table (a):
Reference_Protein_Acessions
- Section Table Options:
- “Tabular Dataset for Table”:
Uniprot
- “Use first line as column names” :
No
- “Specify Column Names (comma-separated list)”:
prot
- param-repeatInsert Table Index:
- “Table Index”:
No
- “Index on Columns”:
prot
- “Table Index”:
comment Comment:
By default, table columns will be named: c1,c2,c3,…,cn (column names for a table must be unique).You can override the default names by entering a comma separated list of names, e.g.
,name1,name2
would rename the second and fifth columns.Check your input file to find the settings which best fits your needs. - “Tabular Dataset for Table”:
“Save the sqlite database in your history”:
No
comment Querying an SQLite Database
Query Tabular can also use an existing SQLite database.Activating
Save the sqlite database in your history
will store the generated database in the history, allowing to reuse it directly.- “SQL Query to generate tabular output”:
- “include query result column headers”:
Yes
- Section Table Options:
- Click Execute and inspect the query results file after it turned green.
The output from this step is that the resultant peptides would be those which doesn’t belong in the Uniprot or cRAP database. The query tabular tool is used again to create a tabular output containing peptides ready for Blast P analysis.
hands_on Hands-on: Query Tabular
Query Tabulartool: with the following parameters:
- Section Filter Dataset Input
- param-repeat“Insert Filter Tabular Input Lines”
- “Filter by”:
skip leading lines
- “Skip lines”:
1
- “Filter by”:
- param-repeat“Insert Filter Tabular Input Lines”
- Section Table Options:
- “Specify Name for Table”:
psm
- “Use first line as column names” :
No
- “Specify Column Names (comma-separated list)”:
id,Proteins,Sequence
- “Only load the columns you have named into database”:
Yes
- “Specify Name for Table”:
- “SQL Query to generate tabular output”:
“include query result column headers”:
Yes
- Click Execute and inspect the query results file after it turned green.
BlastP search is carried out with the PSM report (output from PeptideShaker). Before, BlastP analysis the “Peptides_for_Blast-P_analysis” is first converted from Tabular format to FASTA file format which can be easily read by the BlastP algorithm. This is done with the help of “Tabular to FASTA” conversion tool.The short BlastP uses parameters for short peptide sequences (8-30 aas). Please use the rerun option to look at the parameters used.
hands_on Hands-on: Tabular to FASTA (version 1.1.1)
- Tabular-to-FASTAtool: with the following parameters:
- “Title column”:
1
- “Sequence Column”:
2
- “Title column”:
The output FASTA file is going to be subjected to BLAST-P analysis.
comment Comment: Tool Versions
The tools are subjected to changes while being upgraded.Thus, running the workflow provided with the tutorial, the user might need to make sure they are using the latest version including updating the parameters.
This completes the walkthrough of the proteogenomics database search workflow. This tutorial is a guide to perform database searching with mass spectronetry files and have peptides ready for Blast-P analysis, you can perform follow up analysis using the next GTN “Proteogenomics Novel Peptide Analysis”. Researchers can use this workflow with their data also, please note that the tool parameters, reference genomes and the workflow will be needed to be modified accordingly.
This workflow was developed by the Galaxy-P team at the University of Minnesota. For more information about Galaxy-P or our ongoing work, please visit us at galaxyp.org
keypoints Key points
With SearchGUI and PeptideShaker you can gain access to multiple search engines
Further information, including links to documentation and original publications, regarding the tools, analysis techniques and the interpretation of results described in this tutorial can be found here.
Did you use this material as an instructor? Feel free to give us feedback on how it went.
- Subina Mehta, Timothy J. Griffin, Pratik Jagtap, Ray Sajulga, James Johnson, Praveen Kumar, 2021 Proteogenomics 2: Database Search (Galaxy Training Materials). /training-material/topics/proteomics/tutorials/proteogenomics-dbsearch/tutorial.html Online; accessed TODAY
- Batut et al., 2018 Community-Driven Data Analysis Training for Biology Cell Systems 10.1016/j.cels.2018.05.012
details BibTeX
After 3 years of work, SearchGUI (SG),a tool that performs protein identification using various search engines, and PeptideShaker (PS) for protein identification (which uses SearchGUI results) have been deeply updated and released in new major versions: 4.0.7 and 2.0.5, respectively.Furthermore, two new utility tools have been added to make the overall workflow more flexible:
- Identification parameters for creating parameter files reusable by both SG and PS.
- FastaCLI, which appends decoy sequences to a FASTA file in a format readable by SG and PS.
The number of changes and improvements made is large, a few highlights:
SearchGUI
- Support for spectrum files in mzML format.
- Search engines updated to the most recent versions. Options updated accordingly.
- FASTA file is now out of the search parameters.
- Many more supported search engines parameters.
- More detailed release notes.
Searchguide.level3.com Remove
PeptideShaker
Searchguide Level 3 Mac
- Support for spectrum files in mzML format.
- Internal refactoring, new internal database and more efficient data management.
- Reduced file sizes and memory load.
- More detailed release notes.
Searchguide
Thanks, Carlos Horro Marcos for your awesome work on the Galaxy integration.