We thank Vivian Delcourt for his help, discussions and advice on this work. X.R. is a member of the Fonds de Recherche du Qu&#233;bec Sant&#233; (FRQS)-supported Centre de Recherche du Centre Hospitalier Universitaire de Sherbrooke. This research was supported by a Canada Research Chair in Functional Proteomics and Discovery of Novel Proteins to X.R. and CIHR grant MOP-137056. We thank the team at Calcul Qu&#233;bec and Compute Canada for their support with the use of the supercomputer mp2 from Universit&#233; de Sherbrooke. Operation of the mp2 supercomputer is funded by the Canada Foundation of Innovation (CFI), le minist&#232;re de l&#39;&#201;conomie, de la science et de l&#39;innovation du Qu&#233;bec (MESI) and les Fonds de Recherche du Qu&#233;bec - Nature et technologies (FRQ-NT). The Galaxy server that was used for some proteomics calculations is in part funded by Collaborative Research Centre 992 Medical Epigenetics (DFG grant SFB 992/1 2012) and German Federal Ministry of Education and Research (BMBF grants 031 A538A/A538C RBC, 031L0101B/031L0101C de.NBI-epi, 031L0106 de.STAIR (de.NBI)).

1. OpenProt database download
NOTE: Custom databases based on RNA-seq data for example can also be obtained and the procedure is detailed in the second section of this protocol. If a custom database is needed, please skip to the next section.
<ol>
	<li>Go to the OpenProt website:&#160;www.openprot.org and open the Downloads page using the link from the top page menu.</li>
	<li>Click on the species of interest based on the analyzed experimental data.</li>
	<li>Click on the protein type desired. 
	NOTE: OpenProt offers three classifications: RefProt, Isoforms and AltProt. As shown in Figure 1, this parameter will vary based on the research objective.
	<ol>
		<li>Click on RefProt alone to generate files containing only known proteins.</li>
		<li>Click on AltProt and Isoforms to generate files containing only novel proteins - either novel isoforms of known proteins (Isoforms) or coded by an alternative ORF (AltProts). Please note that OpenProt enforces a minimum ORF length of 30 codons15.</li>
		<li>Click on AltProts, Isoforms and RefProts to generate files containing all protein types present in the OpenProt database - known and novel proteins.</li>
	</ol>
	</li>
	<li>If available, click on the annotation from which protein sequences are drawn. 
	NOTE: OpenProt offers a more exhaustive proteomic landscape by combining multiple annotations. Transcriptome annotations have a minimal overlap; thus, the selected annotation can substantially affect the visualized proteomic profile15,23.</li>
	<li>Click on the level of supporting evidence necessary for protein consideration. As shown in Figure 1, this parameter will vary based on the research objective.
	<ol>
		<li>Click on minimum of two unique peptides detected to generate files containing only the most confident proteins. 
		NOTE: A criterion of two unique peptides is currently considered a gold standard in proteomics for protein expression. If the experimental aim is to detect known and well-supported proteins, the use of this parameter is recommended.</li>
		<li>Click on minimum of one unique peptides detected to generate files containing proteins that have already been seen at least once among the mass spectrometry experiments re-analyzed by OpenProt. 
		NOTE: This allows for consideration of the shorter length of AltProts and the probability that some of them may contain only one unique tryptic peptide8,11.</li>
		<li>Click on all predicted to generate files containing all of OpenProt predictions. 
		NOTE: This setting is recommended only if the experimental aim is to discover novel proteins (Figure 1). The subsequent substantial increase in the search space calls for an adapted analysis pipeline as discussed below7,15.</li>
	</ol>
	</li>
	<li>Click on the desired file format to download. For proteomic analyses, choose the Fasta (protein) file. The readme file contains all necessary information on the file format.</li>
</ol>
2. Custom OpenProt database download
NOTE: This section details how to obtain a custom database. If no custom database is needed, skip to the next section.
<ol>
	<li>Go to the OpenProt website (www.openprot.org) and open the Search page using the link from the top page menu.</li>
	<li>Click on the species of interest based on the experimental data analyzed.</li>
	<li>Enter a list of genes or transcripts of interest.
	<ol>
		<li>When using a list of genes, enter it in the Gene query box.</li>
		<li>When using a list of transcripts, enter it in the Transcript query box.</li>
	</ol>
	</li>
	<li>Tick any box that applies to the desired database.
	<ol>
		<li>Do not click on any box to obtain a table containing all types of protein supported by OpenProt: RefProt, Isoforms and AltProts.</li>
		<li>Click on Show only proteins with experimental evidence to obtain a table containing all types of proteins (RefProts, Isoforms and AltProts) that have been detected at least once by MS and/or for which translation evidence has been collected from ribosome profiling data.</li>
		<li>Similarly, click on Show only proteins detected by MS or on Show only proteins detected by ribosome profiling to obtain a table containing all types of proteins that have been detected at least once by MS or by ribosome profiling respectively.</li>
		<li>Click on Show only AltProts or on Show only isoforms to obtain a table containing only AltProts or only Isoforms respectively.</li>
		<li>Click on both Show only AltProts and Show only Isoforms to obtain a table containing both types of proteins. 
		NOTE: All combinations of filters are possible.</li>
	</ol>
	</li>
	<li>Once all desired parameters are set, click on Search. The table output will appear below the search query fields.</li>
	<li>Click on the Download Fasta button at the right top corner of the output table. This will generate a Fasta file containing all proteins resulting from the queried list of genes or transcripts.</li>
	<li>Please note that for computational reasons, OpenProt holds a maximum of 2,000 elements to be queried (genes or transcripts) at a time. In the event of a list above that limit, several fasta can be generated and then concatenated (as detailed below); or simply download the whole OpenProt database and filter the obtained file as desired.
	<ol>
		<li>Bin the whole list of genes or transcripts into sub-lists of 2,000 entries or less. For each sub-list, download a Fasta file as described above (step 3.3 to 3.6).</li>
		<li>Log in to the European Galaxy instance (or any other instance where proteomics tools are available),&#160;https://usegalaxy.eu/.</li>
		<li>Create a new history and import all of the downloaded OpenProt databases (one per sub-list of genes or transcripts) by clicking on the upload logo at the left top of the screen.</li>
		<li>Use the Fasta Merge Files and Filter Unique Sequences tool developed by the GalaxyP developers (https://github.com/galaxyproteomics/). Select the Merge all Fasta option and input all of the imported OpenProt databases. 
		NOTE: Each tool can be searched by using the query box on the left side of the screen</li>
		<li>Select the accession only option to assess sequence unicity and copy the OpenProt identifier parse rule (&#62;(.*)\|), then click on Execute.</li>
		<li>Note that all files have been concatenated into a unique Fasta file with no redundancy that now appears in the history panel on the right side of the screen. This constitutes the working database.</li>
	</ol>
	</li>
</ol>
3. Database handling
NOTE: From now on, the Galaxy platform will be used, but the same principles can be applied to other proteomic software.
<ol>
	<li>Log in to the European Galaxy instance (or any other instance where proteomics tools are available),&#160;https://usegalaxy.eu/.</li>
	<li>Create a new history and import the downloaded OpenProt database by clicking on the upload logo at the left top of the screen.</li>
	<li>Go to the workflow page and import the Database Handling workflow (Supplementary Material S1) by clicking on the upload logo at the left top of the middle panel.</li>
	<li>Click on Run the workflow and select the imported OpenProt database as input. 
	NOTE: This workflow will append the CRAPome repository to the OpenProt fasta and generate decoy sequences (reverse sequences)24. If a shuffle decoy list is desired, it can be done by changing this parameter on the DecoyDatabase tool.</li>
	<li>Rename the obtained Fasta file to something meaningful. The database is ready to be used for proteomics analyses.</li>
</ol>
4. Mass spectrometry file preparation
NOTE: Most of the proteomics tools available on Galaxy instances use the mzML format, and peptide search engines prefer data in centroid mode.
<ol>
	<li>Open the freely available MSConvert tool from the ProteoWizard suite and upload the data file to be analyzed25.</li>
	<li>Choose the directory for the output and the desired file format to mzML.</li>
	<li>Set a peak picking filter using the wavelet based algorithm (CWT) on MS1 and MS2 levels, and start the conversion26.</li>
</ol>
5. Peptide and protein identification/quantification
NOTE: This part of the pipeline uses tools from the OpenMS suite, a versatile and easy-to-use framework18.
<ol>
	<li>Log in to the European Galaxy instance (or any other instance where proteomics tools are available),&#160;https://usegalaxy.eu/.</li>
	<li>Create a new history and transfer the previously created database (step 3.5) to this new history with a drag-and-drop.</li>
	<li>Import the transformed mzML data file (step 4.3) by clicking on the Upload logo at the left top of the screen.</li>
	<li>Go to the workflow page and import the desired workflow by clicking on the upload logo at the left top of the middle panel. 
	NOTE: MS experiments are differently designed based on the desired final output. Workflows are provided here for two frequent designs: protein identification and protein quantification based on stable isotope labeling (SIL). However, the Galaxy instance contains many other tools that will support other types of proteomic analyses27,28.
	<ol>
		<li>For a protein identification design, import the workflow provided in Supplementary Material S2.&#160;When using this workflow, please do not use the zlip compression when converting your files (step 4.2)</li>
		<li>For a protein quantification based on stable isotope labeling design, import the workflow provided in Supplementary Material S3.</li>
	</ol>
	</li>
	<li>Select run the workflow and review the different parameters.
	<ol>
		<li>Select the imported mzML data file as input, and the previously created database (step 3.5) as the database Fasta file.</li>
		<li>Since the workflow uses the X!Tandem search engine, import the X!Tandem default configuration file (provided in Supplementary Material S4)29 by clicking on the upload logo at the left top of the screen.</li>
		<li>The workflow uses multiple search engines (MS-GF+ and X!Tandem). Append other search engines or choose a single one simply by adding or removing the tools from the workflow30,31. 
		NOTE: Using multiple search engines is recommended as it increases sensibility and sensitivity of the analysis32.</li>
		<li>In order to account for the substantial increase in size when using the whole OpenProt database, use a stringent FDR15. By default, the provided workflow is set for a 0.001% FDR, adequate for the use of the whole OpenProt database. For other databases, this can be edited to any desired value. 
		NOTE: Be sure to adapt the parameters of the different tools depending on the mass spectrometer used and the experimental protocol (precursor ion and fragment error, fixed and variable modifications, used enzyme, etc.).</li>
	</ol>
	</li>
	<li>Optionally, download output for each step of the workflow for storage or quality control analysis by clicking on the chosen step from the history panel, then clicking on the Save logo that will appear underneath.</li>
</ol>
6. Quality control
NOTE: Because MS-based proteomics is the result of a complex process where each step needs to be optimized to produce reproducible results, quality control is a necessary procedure in the workflow33.
<ol>
	<li>Several metrics are common benchmark of performance, such as the number of peptide-spectrum matches (PSM), the number of identified peptides and proteins. Run the File Info tool on the IDFilter output (indicated in green in Figure 2) to provide such metrics.</li>
	<li>Although not applicable to every identification, especially with large datasets, reports of novel proteins should always be carefully evaluated. Inspection of the protein score, the sequence coverage, and the spectra supporting the finding is of vital importance. Use the TOPPview tool from the OpenMS framework to do this; it is freely available and well documented18,34,35.</li>
</ol>
7. OpenProt database mining
NOTE: Once a confident identification of a novel protein predicted by OpenProt (accession numbers starting with IP_ for AltProts and II_ for novel Isoforms) has been made, more biological information can be gathered from the OpenProt website15.
<ol>
	<li>Go to the OpenProt website:&#160;www.openprot.org and open the Search page using the link at the top page menu.&#160;</li>
	<li>Click on the species of interest (same as the one in which the protein was identified) and enter the protein accession number in the Protein query box.</li>
	<li>Click on search and a table containing basic information on the queried protein will appear. The table features: the protein length (in amino acid), its molecular weight (kDa) and isoelectric point, supporting experimental evidence by MS or ribosome profiling (Translation Evidence, TE), and functional predictions such as predicted domains and protein orthology (across the 10 species supported by OpenProt, v1.3). The table also contains information about the related gene and transcript and the localization of the protein within the transcript.</li>
	<li>Click on the Details link to gather further information. The newly opened page contains a genome browser which is centered on the queried protein, and information such as the genomic and transcriptomic coordinates and the presence of a Kozak or high-efficiency translation initiation site (TIS) motif36,37.</li>
	<li>Click on the Protein or DNA links from the info tab, to obtain protein or DNA sequences respectively.</li>
	<li>Browse detailed information about MS evidence, ribosome profiling detection, conservation and identified protein domains by clicking onto the top tabs15.</li>
</ol>

<ol>
	<li>Kim, M. S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+draft+map+of+the+human+proteome.">A draft map of the human proteome.</a> Nature. 509 (7502), 575-581 (2014).</li><li>Wilhelm, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Mass-spectrometry-based+draft+of+the+human+proteome.">Mass-spectrometry-based draft of the human proteome.</a> Nature. 509 (7502), 582-587 (2014).</li><li>Hein, M. Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+human+interactome+in+three+quantitative+dimensions+organized+by+stoichiometries+and+abundances.">A human interactome in three quantitative dimensions organized by stoichiometries and abundances.</a> Cell. 163 (3), 712-723 (2015).</li><li>Huttlin, E. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+BioPlex+Network:+A+Systematic+Exploration+of+the+Human+Interactome.">The BioPlex Network: A Systematic Exploration of the Human Interactome.</a> Cell. 162 (2), 425-440 (2015).</li><li>Huttlin, E. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Architecture+of+the+human+interactome+defines+protein+communities+and+disease+networks.">Architecture of the human interactome defines protein communities and disease networks.</a> Nature. 545 (7655), 505-509 (2017).</li><li>Kumar, D., Yadav, A. K., Dash, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Choosing+an+Optimal+Database+for+Protein+Identification+from+Tandem+Mass+Spectrometry+Data.">Choosing an Optimal Database for Protein Identification from Tandem Mass Spectrometry Data.</a> Proteome Bioinformatics. , 17-29 (2017).</li><li>Jeong, K., Kim, S., Bandeira, N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=False+discovery+rates+in+spectral+identification.">False discovery rates in spectral identification.</a> BMC Bioinformatics. 13 (Suppl 16), (2012).</li><li>Brunet, M. A., Levesque, S. A., Hunting, D. J., Cohen, A. A., Roucou, X. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Recognition+of+the+polycistronic+nature+of+human+genes+is+critical+to+understanding+the+genotype-phenotype+relationship.">Recognition of the polycistronic nature of human genes is critical to understanding the genotype-phenotype relationship.</a> Genome Research. , (2018).</li><li>Brent, M. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Genome+annotation+past,+present,+and+future:+how+to+define+an+ORF+at+each+locus.">Genome annotation past, present, and future: how to define an ORF at each locus.</a> Genome Research. 15 (12), 1777-1786 (2005).</li><li>Harrow, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=GENCODE:+The+reference+human+genome+annotation+for+The+ENCODE+Project.">GENCODE: The reference human genome annotation for The ENCODE Project.</a> Genome Research. 22 (9), 1760-1774 (2012).</li><li>Samandi, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Deep+transcriptome+annotation+enables+the+discovery+and+functional+characterization+of+cryptic+small+proteins.">Deep transcriptome annotation enables the discovery and functional characterization of cryptic small proteins.</a> eLife. 6, e27860 (2017).</li><li>Saghatelian, A., Couso, J. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Discovery+and+characterization+of+smORF-encoded+bioactive+polypeptides.">Discovery and characterization of smORF-encoded bioactive polypeptides.</a> Nature Chemical Biology. 11 (12), 909-916 (2015).</li><li>Delcourt, V., Staskevicius, A., Salzet, M., Fournier, I., Roucou, X. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Small+Proteins+Encoded+by+Unannotated+ORFs+are+Rising+Stars+of+the+Proteome,+Confirming+Shortcomings+in+Genome+Annotations+and+Current+Vision+of+an+mRNA.">Small Proteins Encoded by Unannotated ORFs are Rising Stars of the Proteome, Confirming Shortcomings in Genome Annotations and Current Vision of an mRNA.</a> Proteomics. , (2017).</li><li>Plaza, S., Menschaert, G., Payre, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=In+Search+of+Lost+Small+Peptides.">In Search of Lost Small Peptides.</a> Annual Review of Cell and Developmental Biology. 33 (1), (2017).</li><li>Brunet, M. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=OpenProt:+a+more+comprehensive+guide+to+explore+eukaryotic+coding+potential+and+proteomes.">OpenProt: a more comprehensive guide to explore eukaryotic coding potential and proteomes.</a> Nucleic Acids Research. , (2018).</li><li>Afgan, E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+Galaxy+platform+for+accessible,+reproducible+and+collaborative+biomedical+analyses:+2016+update.">The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2016 update.</a> Nucleic Acids Research. 44 (W1), W3-W10 (2016).</li><li>Afgan, E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+Galaxy+platform+for+accessible,+reproducible+and+collaborative+biomedical+analyses:+2018+update.">The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update.</a> Nucleic Acids Research. 46, W537-W544 (2018).</li><li>Sturm, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=OpenMS+&#8211;+An+open-source+software+framework+for+mass+spectrometry.">OpenMS &#8211; An open-source software framework for mass spectrometry.</a> BMC Bioinformatics. 9 (1), 163 (2008).</li><li>Carithers, L. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+Novel+Approach+to+High-Quality+Postmortem+Tissue+Procurement:+The+GTEx+Project.">A Novel Approach to High-Quality Postmortem Tissue Procurement: The GTEx Project.</a> Biopreservation and Biobanking. 13 (5), 311-319 (2015).</li><li>Aebersold, R., Mann, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Mass+spectrometry-based+proteomics.">Mass spectrometry-based proteomics.</a> Nature. 422 (6928), 6928 (2003).</li><li>Domon, B., Aebersold, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Mass+Spectrometry+and+Protein+Analysis.">Mass Spectrometry and Protein Analysis.</a> Science. 312 (5771), 212-217 (2006).</li><li>Hu, J., Coombes, K. R., Morris, J. S., Baggerly, K. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+importance+of+experimental+design+in+proteomic+mass+spectrometry+experiments:+Some+cautionary+tales.">The importance of experimental design in proteomic mass spectrometry experiments: Some cautionary tales.</a> Briefings in Functional Genomics. 3 (4), 322-331 (2005).</li><li>Wu, P. Y., Phan, J. H., Wang, M. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Assessing+the+impact+of+human+genome+annotation+choice+on+RNA-seq+expression+estimates.">Assessing the impact of human genome annotation choice on RNA-seq expression estimates.</a> BMC Bioinformatics. 14 (11), S8 (2013).</li><li>Mellacheruvu, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+CRAPome:+a+contaminant+repository+for+affinity+purification-mass+spectrometry+data.">The CRAPome: a contaminant repository for affinity purification-mass spectrometry data.</a> Nature Methods. 10 (8), 730-736 (2013).</li><li>Adusumilli, R., Mallick, P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Data+Conversion+with+ProteoWizard+msConvert.">Data Conversion with ProteoWizard msConvert.</a> Proteomics: Methods and Protocols. , 339-368 (2017).</li><li>French, W. R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Wavelet-Based+Peak+Detection+and+a+New+Charge+Inference+Procedure+for+MS/MS+Implemented+in+ProteoWizard&#8217;s+msConvert.">Wavelet-Based Peak Detection and a New Charge Inference Procedure for MS/MS Implemented in ProteoWizard&#8217;s msConvert.</a> Journal of Proteome Research. 14 (2), 1299-1307 (2015).</li><li>Kuenzi, B. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=APOSTL:+An+Interactive+Galaxy+Pipeline+for+Reproducible+Analysis+of+Affinity+Proteomics+Data.">APOSTL: An Interactive Galaxy Pipeline for Reproducible Analysis of Affinity Proteomics Data.</a> Journal of Proteome Research. 15 (12), 4747-4754 (2016).</li><li>Hoekman, B., Breitling, R., Suits, F., Bischoff, R., Horvatovich, P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=msCompare:+a+framework+for+quantitative+analysis+of+label-free+LC-MS+data+for+comparative+candidate+biomarker+studies.">msCompare: a framework for quantitative analysis of label-free LC-MS data for comparative candidate biomarker studies.</a> Molecular &amp; Cellular Proteomics: MCP. 11 (6), (2012).</li><li>Bjornson, R. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=X!!Tandem,+an+improved+method+for+running+X!tandem+in+parallel+on+collections+of+commodity+computers.">X!!Tandem, an improved method for running X!tandem in parallel on collections of commodity computers.</a> Journal of Proteome Research. 7 (1), 293-299 (2008).</li><li>Kim, S., Pevzner, P. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=MS-GF++makes+progress+towards+a+universal+database+search+tool+for+proteomics.">MS-GF+ makes progress towards a universal database search tool for proteomics.</a> Nature Communications. 5, 5277 (2014).</li><li>Vaudel, M., Barsnes, H., Berven, F. S., Sickmann, A., Martens, L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=SearchGUI:+An+open-source+graphical+user+interface+for+simultaneous+OMSSA+and+X!Tandem+searches.">SearchGUI: An open-source graphical user interface for simultaneous OMSSA and X!Tandem searches.</a> Proteomics. 11 (5), 996-999 (2011).</li><li>Shteynberg, D., Nesvizhskii, A. I., Moritz, R. L., Deutsch, E. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Combining+results+of+multiple+search+engines+in+proteomics.">Combining results of multiple search engines in proteomics.</a> Molecular &amp; Cellular Proteomics: MCP. 12 (9), 2383-2393 (2013).</li><li>Bittremieux, W., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Quality+control+in+mass+spectrometry-based+proteomics.">Quality control in mass spectrometry-based proteomics.</a> Mass Spectrometry Reviews. 37 (5), 697-711 (2018).</li><li>Bertsch, A., Gr&#246;pl, C., Reinert, K., Kohlbacher, O. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=OpenMS+and+TOPP:+Open+Source+Software+for+LC-MS+Data+Analysis.">OpenMS and TOPP: Open Source Software for LC-MS Data Analysis.</a> Data Mining in Proteomics: From Standards to Applications. , 353-367 (2011).</li><li>Pfeuffer, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=OpenMS+&#8211;+A+platform+for+reproducible+analysis+of+mass+spectrometry+data.">OpenMS &#8211; A platform for reproducible analysis of mass spectrometry data.</a> Journal of Biotechnology. 261, 142-148 (2017).</li><li>Kozak, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Pushing+the+limits+of+the+scanning+mechanism+for+initiation+of+translation.">Pushing the limits of the scanning mechanism for initiation of translation.</a> Gene. 299 (1-2), 1-34 (2002).</li><li>Noderer, W. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Quantitative+analysis+of+mammalian+translation+initiation+sites+by+FACS-seq.">Quantitative analysis of mammalian translation initiation sites by FACS-seq.</a> Molecular Systems Biology. 10, 748 (2014).</li><li>Eyckerman, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Intelligent+Mixing+of+Proteomes+for+Elimination+of+False+Positives+in+Affinity+Purification-Mass+Spectrometry.">Intelligent Mixing of Proteomes for Elimination of False Positives in Affinity Purification-Mass Spectrometry.</a> Journal of Proteome Research. 15 (10), 3929-3937 (2016).</li><li>Vizca&#237;no, J. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=2016+update+of+the+PRIDE+database+and+its+related+tools.">2016 update of the PRIDE database and its related tools.</a> Nucleic Acids Research. 44 (D1), D447-D456 (2016).</li><li>Bateman, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=UniProt:+the+universal+protein+knowledgebase.">UniProt: the universal protein knowledgebase.</a> Nucleic Acids Research. 45 (D1), D158-D169 (2017).</li><li>The Gene Ontology Consortium Expansion of the Gene Ontology knowledgebase and resources. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Expansion+of+the+Gene+Ontology+knowledgebase+and+resources.">Expansion of the Gene Ontology knowledgebase and resources.</a> Nucleic Acids Research. 45 (D1), D331-D338 (2017).</li><li>O&#8217;Leary, N. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Reference+sequence+(RefSeq)+database+at+NCBI:+current+status,+taxonomic+expansion,+and+functional+annotation.">Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation.</a> Nucleic Acids Research. 44, D733-D745 (2016).</li><li>Zerbino, D. R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Ensembl+2018.">Ensembl 2018.</a> Nucleic Acids Research. 46 (D1), D754-D761 (2018).</li><li>Andreev, D. E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Translation+of+5&#8217;+leaders+is+pervasive+in+genes+resistant+to+eIF2+repression.">Translation of 5&#8217; leaders is pervasive in genes resistant to eIF2 repression.</a> eLife. 4, e03971 (2015).</li><li>Jackson, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+translation+of+non-canonical+open+reading+frames+controls+mucosal+immunity.">The translation of non-canonical open reading frames controls mucosal immunity.</a> Nature. 564, 434-438 (2018).</li></ol>

The authors declare no conflict of interests.

When analyzing data from mass spectrometers, the quality of protein identification partly relies on the accuracy of the used database6,20. Current approaches traditionally use UniProtKB databases, yet these support the genome annotation model of a single ORF per transcript and a minimum length of 100 codons (with the exception of previously demonstrated examples)40. Multiple studies relate the shortcomings of such databases with the discovery of functional ORFs from allegedly non-coding regions8,11,12,13. Now, OpenProt allows for more exhaustive protein identification as it draws protein sequences from multiple transcriptome annotations. OpenProt retrieves NCBI RefSeq (GRCh38.p7) and Ensembl (GRCh38.83) transcriptomes and UniProtKB annotations (UniProtKB-SwissProt, 2017-09-27)40,42,43. As current annotations present little overlap, OpenProt thus displays a more exhaustive view of the potential proteomic landscape than when limited to one annotation15.
Furthermore, as OpenProt enforces a polycistronic model, it allows for multiple protein annotations per transcript. For statistical and computational reasons, OpenProt still holds a minimum length threshold of 30 codons15. Yet, it predicts thousands of novel protein sequences, thereby widening the scope of possibilities for protein identification. With this approach, OpenProt supports proteomic discoveries in a more systematic manner.
The quality of protein identification can also be affected by the parameters that are used. MS-based proteomics analyses typically hold a 1% protein FDR. However, the whole OpenProt database contains about 6 times more entries (Figure 1). To account for this substantial increase in the search space, we recommend using a more stringent FDR of 0.001%. This parameter was optimized using benchmark studies and manual evaluation of randomly selected spectra15. False positive are still a possibility, though, and we encourage thorough inspection and validation of supporting evidence for a novel protein. A recommended standard could be the identification of a protein from two different MS runs, as background data and false positives vary between datasets15.
The pipeline provided here and used for the case study can be modified as pleased to fit the experimental design and parameters. We would recommend using multiple search engines as it increases sensibility and sensitivity of peptide identification32. Furthermore, we encourage using the database corresponding best to the experimental aim (Figure 1). As using the whole OpenProt database comes with a stringent FDR, true identifications may be lost. Thus, the whole database should be intended for discovery of novel proteins, whilst classical proteomics profiling should be using the smaller OpenProt databases (such as OpenProt_2pep used in the case study above).
OpenProt currently predicts sequences starting with an ATG codon, whereas several studies highlighted translation initiation at other codons44,45. When a novel protein is identified by one or several unique peptides, it is possible the true initiation codon is not the presumed ATG. Users can look for translation evidence on the OpenProt website. Currently, OpenProt only reports translation events if they concern the entire predicted protein sequence (100% overlap)15. Thus, absence of translation evidence would not mean the protein is not translated, but that the start codon may not be the alleged ATG.
Despite its current limitations, OpenProt offers a more exhaustive view of eukaryotic genomes&#39; coding potential. OpenProt databases foster proteomic discoveries and the understanding of proteomic functions and interactions. Future developments of the OpenProt database will include annotation of other species, translation evidence from non-ATG start codon and development of a pipeline to include novel proteins in whole genome and exome sequencing studies.

Over the past decades, mass spectrometry (MS-)based proteomics has become the golden technique to decipher proteomes of eukaryotic cells1,2,3,4,5. This method relies on current genome annotations to generate a reference protein sequence database that outlines the scope of possibilities6,7,8. However, genome annotations hold arbitrary criteria for ORF annotation, such as a minimum length of 100 codons and a single ORF per transcript9,10. An increasing number of studies challenge the current annotation model and report discoveries of unannotated functional ORFs in eukaryotic genomes8,11,12,13,14. These novel proteins are found encoded in allegedly non-coding RNAs, in the 5&#39; or 3&#39; untranslated regions (UTR) of mRNAs, or overlapping the canonical coding sequence (cCDS) in an alternative frame. Although most of these discoveries have been serendipitous, they demonstrate the caveats of current genome annotations and the polycistronic nature of eukaryotic genes8.
Here, we highlight the use of OpenProt databases for MS-based proteomics. OpenProt is the first database to hold a polycistronic annotation model for eukaryotic transcriptomes. It is freely available at www.openprot.org15. A proportion of these predicted ORFs would be random and non-functional, which is why OpenProt cumulates experimental and functional evidence to increase confidence. Experimental evidence include protein expression (by MS) and translation evidence (by ribosome profiling)15. Functional evidence include protein orthology (with an In-Paranoid like approach) and functional domain prediction15.
OpenProt offers the possibility to download several databases, from containing only well-supported proteins to custom-made databases. Here, we will present a pipeline for the use of OpenProt databases and will offer insights into which database to choose considering the experimental aim. The proteomics analysis pipeline presented here is supported by the Galaxy framework as it is open-access and easy-to-use, but the databases can work with any workflow16,17,18. We will also present how to use the OpenProt website for gathering further information on novel proteins detected by MS. Using OpenProt databases will provide a more exhaustive view of the proteomic landscape and will foster proteomics and biomarkers discoveries in a more systematic way than current methods.
This protocol highlights the use of OpenProt databases15 when interrogating MS datasets; it will not review the design of the experiment itself, which has been thoroughly reviewed elsewhere20,21,22. In an effort to remain fully open-source, the protocol is freely available (Supplementary Material S1-S4). For easier reading, all terms used in OpenProt and hereby throughout this protocol are defined in Table 1.

<table border="0" cellpadding="0" cellspacing="0" style="width:686px;" width="685">
	<tbody>
		<tr height="21">
			<td height="21" style="height:21px;width:223px;">OpenProt website</td>
			<td style="width:126px;">open source</td>
			<td style="width:166px;">n/a</td>
			<td style="width:172px;">www.openprot.org</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:223px;">Galaxy Server</td>
			<td style="width:126px;">open source</td>
			<td style="width:166px;">n/a</td>
			<td>https://usegalaxy.eu/</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;">TOPPview software</td>
			<td style="width:126px;">open source</td>
			<td style="width:166px;">n/a</td>
			<td>www.openms.de</td>
		</tr>
	</tbody>
</table>

mass spectrometry-based proteomics analyses using the openprot database to unveil novel proteins translated from non-canonical open reading frames

Genome annotation is central to today's proteomic research as it draws the outlines of the proteomic landscape. Traditional models of open reading frame (ORF) annotation impose two arbitrary criteria: a minimum length of 100 codons and a single ORF per transcript. However, a growing number of studies report expression of proteins from allegedly non-coding regions, challenging the accuracy of current genome annotations. These novel proteins were found encoded either within non-coding RNAs, 5' or 3' untranslated regions (UTRs) of mRNAs, or overlapping a known coding sequence (CDS) in an alternative ORF. OpenProt is the first database that enforces a polycistronic model for eukaryotic genomes, allowing annotation of multiple ORFs per transcript. OpenProt is freely accessible and offers custom downloads of protein sequences across 10 species. Using OpenProt database for proteomic experiments enables novel proteins discovery and highlights the polycistronic nature of eukaryotic genes. The size of OpenProt database (all predicted proteins) is substantial and need be taken in account for the analysis. However, with appropriate false discovery rate (FDR) settings or the use of a restricted OpenProt database, users will gain a more realistic view of the proteomic landscape. Overall, OpenProt is a freely available tool that will foster proteomic discoveries.

The workflow described above was applied to a MS dataset available on the PRIDE repository38,39. The original study developed a method (iMixPro), using stable isotope labeling of amino acids in cell culture (SILAC), to eliminate false positives from affinity-purification MS (AP-MS) experiments38. In brief, an AP-MS experiment consists of using beads-bound antibodies to fetch a protein of interest (bait) and its interactors (preys). The collected proteins are then digested and prepared for MS. The sample preparation method and the instrument settings are described in the original study and on the PRIDE repository (PXD004246). A challenge in such experiments is the abundance of false positives, notably from proteins binding to the beads but not the bait. Here, we used SILAC to generate different isotope ratios between true preys and false positives: 3 control samples (no bait) cultured in light medium, 1 sample expressing the bait cultured in light medium, and 1 sample expressing the bait cultured in heavy medium are processed with the beads and further mass spectrometry analysis. With such design, non-specific proteins binding to the beads will have an heavy-to-light ratio of 1:4; when true preys will have a ratio of 1:138.
We re-analyzed their AP-MS data using the OpenProt database; the baits included three endogenous proteins (PTPN14, JIP3 and IQGAP1), and two over-expressed proteins (RAF1 and RNF41). Since the experiments used SILAC, the Galaxy workflow for protein quantification was used (Supplementary Material S3, Figure 2). The workflow was run using the whole OpenProt database (OpenProt_all) or a restricted OpenProt database (OpenProt_2pep, including only proteins previously detected with a minimum of two unique peptides).
Protein identification and quantification were good and reproducible across the different used databases. As shown in Figure 3, most proteins identified in the original paper were also identified using either the OpenProt_2pep or OpenProt_all database (a detailed list is available in Supplementary Material S5). This result shows that the pipeline described here and the OpenProt databases are able to produce protein identification and quantification comparable to that of current procedures based on the UniProtKB databases40. However, the use of OpenProt databases has the unique advantage of allowing detection of novel and previously undetectable proteins, as demonstrated in this case study.
11 well-supported proteins (1 Isoform and 10 AltProts), yet currently not annotated in databases, were identified across all datasets, with confident peptides, using the OpenProt_2pep database (all protein accessions, along with the number of supporting peptides, are available in Supplementary Material S5). This database allows the use of a traditional 1% FDR as the search space increase remains moderate. These 11 proteins were not identified in the original study as they were absent from the database.
29 novel proteins (16 isoforms and 13 AltProts) were discovered across all datasets, with confident peptides, using the OpenProt_all database (all protein accessions, along with the number of supporting peptides, are available in Supplementary Material S6). As shown in Figure 3, the recommended stringent FDR did not affect the most confident protein identifications, although it did decrease the total number of identified proteins. Comparatively to the OpenProt_2pep database, a higher number of novel proteins can be confidently identified. All of these novel proteins are absent from the OpenProt_2pep database. This highlights the crucial role of the chosen database for MS-based proteomics.
One novel protein was discovered as an interactor of the RAF1 protein (IP_637643). Using the OpenProt website, one can see this protein had not been detected by MS nor ribosome profiling until now (OpenProt v1.3). The protein is 46 amino acids long and can only give two unique peptides upon tryptic digestion. The peptide detected in the RAF1 AP-MS dataset (fraction 18) had a good quality spectrum, as shown in Figure 4, and displayed a heavy-to-light ratio of 1,09. The protein is encoded in the NANOGNBP1 gene, which is a pseudogene of NANOGNB. The transcript (ENST00000448444), currently annotated as non-coding, was detected across several tissues according to the GTEx portal40. The protein contains a predicted functional domain associated with DNA binding (Gene Ontology GO:0003677)41.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/59589/59589fig1.jpg" /> 
Figure 1: Database choice for proteomics analyses chart. Analyses of MS data, notably the database choice, depend on the research objectives. Three common objectives are outlined in blue (classic proteomic pipeline), green (exhaustive proteomic search) and orange (proteomic discovery). Each objective depends on an appropriate database and pipeline. A single identification tool may be used for an exhaustive and classical proteomics pipelines. For the proteomic discovery pipeline, we strongly recommend using multiple identification engines. Recommended FDRs are indicated in red, and protein database sizes are indicated in grey boxes. <a href="https://www.jove.com/files/ftp_upload/59589/59589fig1large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/59589/59589fig2.jpg" /> 
Figure 2: Graphical representation of the Galaxy workflow used. Step-by-step representation of the proteomic analysis workflow used for re-analysis of Eyckerman et al. data38. Input files, peptide search, and protein quantification are indicated by orange boxes. Blue boxes correspond to the tools used and grey boxes correspond to the output files generated. The different search engines (MS-GF+ and X!Tandem) are indicated by different colors (respectively red and purple) as well as the arrows indicating their necessary inputs and outputs. The green box highlights the tool generating a list of protein identifications. When multiple outputs are generated, the one used for downstream steps is indicated as the closest to the arrow. This workflow is freely available in Supplementary Material S2. The X!Tandem default parameters configuration file is available in Supplementary Material S4. <a href="https://www.jove.com/files/ftp_upload/59589/59589fig2large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/59589/59589fig3.jpg" /> 
Figure 3: Comparison of interactor identification per bait using different databases. Venn diagrams of protein identifications using the most confident OpenProt database (in orange, supporting evidence of minimum 2 unique peptides, OpenProt_2pep) with a 1% FDR, or the whole OpenProt database (in blue, OpenProt_all) with a 0.001% FDR, or as reported in the original paper (in grey)38. Each diagram corresponds to identified interactors for the mentioned bait: RAF1, RNF41, PTPN14, JIP3 and IQGAP1. <a href="https://www.jove.com/files/ftp_upload/59589/59589fig3large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/59589/59589fig4.jpg" /> 
Figure 4: MS/MS spectrum of identified MDNLWAK(13C6) peptide from novel protein IP_637643. Intensity is relative (0 to 100%). Selected peaks are indicated in red, y ions annotations are in dark red and b ions annotations in green. Extracted from the TOPPview software34. Precursor Error = 2.70 ppm, PEP score = 0.12. <a href="https://www.jove.com/files/ftp_upload/59589/59589fig4large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>Term</td>
			<td>Definition</td>
			<td>Reference</td>
		</tr>
		<tr>
			<td>Alternative ORF (AltORF)</td>
			<td>non-canonical ORF currently not annotated in genome annotations, but annotated in OpenProt.</td>
			<td>15</td>
		</tr>
		<tr>
			<td>Reference ORF (RefORF)</td>
			<td>canonical ORF annotated in genome annotations and OpenProt.</td>
			<td>15</td>
		</tr>
		<tr>
			<td>Alternative protein (AltProt)</td>
			<td>novel protein coded by an AltORF, with no significant similarity with a RefProt. Accession prefix: IP_.</td>
			<td>15</td>
		</tr>
		<tr>
			<td>Reference protein (RefProt)</td>
			<td>protein currently annotated in protein sequence databases such as UniProtKB, Ensembl or NCBI RefSeq, and also in OpenProt.</td>
			<td>15</td>
		</tr>
		<tr>
			<td>Novel Isoform</td>
			<td>novel protein coded by an AltORF, with a significant similarity with a RefProt. Accession prefix: II_.</td>
			<td>15</td>
		</tr>
		<tr>
			<td>OpenProt_2pep database</td>
			<td>contains the sequence of all RefProts and novel proteins predicted by OpenProt, already detected with a minimum of 2 unique peptides.</td>
			<td>15</td>
		</tr>
		<tr>
			<td>OpenProt_1pep database</td>
			<td>contains the sequence of all RefProts and novel proteins predicted by OpenProt, already detected with a minimum of 1 unique peptide.</td>
			<td>15</td>
		</tr>
		<tr>
			<td>OpenProt_all database</td>
			<td>contains the sequence of all RefProts and novel proteins predicted by OpenProt.</td>
			<td>15</td>
		</tr>
	</tbody>
</table>
Table 1: Definition of terms used in OpenProt and throughout the protocol
Supplementary Material S1: Galaxy workflow for database handling. This will append the CRAPome and decoy sequences (reverse) to the input database. Output is a Fasta file.&#160;<a href="https://www.jove.com/files/ftp_upload/59589/Supplementary_S1_Galaxy-Workflow_Database_Handling.ga" target="_blank">Please click here to download.</a>
Supplementary Material S2: Galaxy workflow for protein identification. This will identify proteins from a mass spectrometry data file using two search engines (MS-GF+ and X!Tandem). Each parameter can be tuned as desired before running the workflow.&#160;<a href="https://www.jove.com/files/ftp_upload/59589/S2_Galaxy_Pep_Prot_Id.ga" target="_blank">Please click here to download.</a>
Supplementary Material S3: Galaxy workflow for protein quantification using stable isotope labeling (SIL). This will identify and quantify proteins from a mass spectrometry data file using two search engines (MS-GF+ and X!Tandem). Each parameter can be tuned as desired before running the workflow.&#160;<a href="https://www.jove.com/files/ftp_upload/59589/S3_Galaxy_SILAC_Quant.ga" target="_blank">Please click here to download.</a>
Supplementary Material S4: X!Tandem default parameters configuration file. This XML file is necessary for running the X!TandemAdapter tool on the Galaxy platform.&#160;<a href="https://www.jove.com/files/ftp_upload/59589/S4_XTandem_Config.xml" target="_blank">Please click here to download.</a>
Supplementary Material S5: Quantified proteins from iMixPro datasets. Data files from Eyckerman et al. 201638 were processed using OpenProt databases and quantified proteins are listed for each condition. Baits are PTPN14, JIP3, IQGAP1, RAF1 and RNF41. Gene names indicated in green correspond to proteins also identified in the original paper38. Gene names indicated in orange correspond to known interactors according to BioGrid that were not reported in the original paper. Gene names indicated in light blue correspond to novel proteins identified as interactors (the corresponding protein accession number is indicated in brackets). Gene names indicated in light grey and italics correspond to likely contaminants (keratin proteins).&#160;<a href="https://www.jove.com/files/ftp_upload/59589/S5_and_S6__DataMS.xlsx" target="_blank">Please click here to download.</a>
Supplementary Material S6: Identified novel proteins from iMixPro datasets. Data files from Eyckerman et al. 201638 were processed using OpenProt databases and novel identified proteins are listed for each condition. Baits are PTPN14, JIP3, IQGAP1, RAF1 and RNF41. Protein accession numbers are listed, starting with II_ for novel isoforms of a known protein, and with IP_ for novel proteins from an alternative ORF (AltProt).The number of supporting peptides are indicated in brackets.&#160;<a href="https://www.jove.com/files/ftp_upload/59589/S5_and_S6__DataMS.xlsx" target="_blank">Please click here to download.</a>

Watch this Scientific Journal Video about Mass Spectrometry-Based Proteomics Analyses Using the OpenProt Database to Unveil Novel Proteins Translated from Non-Canonical Open Reading Frames at JoVE.com

Mass Spectrometry-Based Proteomics Analyses Using the OpenProt Database to Unveil Novel Proteins Translated from Non-Canonical Open Reading Frames

OpenProt is a freely accessible database that enforces a polycistronic model of eukaryotic genomes. Here, we present a protocol for the use of OpenProt databases when interrogating mass spectrometry datasets. Using OpenProt database for analysis of proteomic experiments allows for discovery of novel and previously undetectable proteins.

Genome annotation is central to today's proteomic research as it draws the outlines of the proteomic landscape. Traditional models of open reading frame ...

mass-spectrometry-based-proteomics-analyses-using-openprot-database

Research

JoVE Journal

Genetics

12.5K Views.  Université de Sherbrooke. OpenProt is a freely accessible database that enforces a polycistronic model of eukaryotic genomes. Here, we present a protocol for the use of OpenProt databases when interrogating mass spectrometry datasets. Using OpenProt database for analysis of proteomic experiments allows for discovery of novel and previously undetectable proteins. 

Video: Mass Spectrometry-Based Proteomics Analyses Using the OpenProt Database to Unveil Novel Proteins Translated from Non-Canonical Open Reading Frames

Subtyping of Campylobacter jejuni ssp. doylei Isolates Using Mass Spectrometry-based PhyloProteomics (MSPP)

MALDI-TOF MS offers the possibility to differentiate some bacteria not only at the species and subspecies level but even below, at the strain level. Allelic isoforms of the detectable biomarker ions result in isolate-specific mass shifts. Mass spectrometry-based phyloproteomics (MSPP) is a novel technique that combines the mass spectrometric detectable biomarker masses in a scheme that allows deduction of phyloproteomic relations from isolate specific mass shifts compared to a genome sequenced reference strain. The deduced amino acid sequences are then used to calculate MSPP-based dendrograms.
Here we describe the workflow of MSPP by typing a Campylobacter jejuni ssp. doylei isolate collection of seven strains. All seven strains were of human origin and multilocus sequence typing (MLST) demonstrated their genetic diversity. MSPP-typing resulted in seven different MSPP sequence types, sufficiently reflecting their phylogenetic relations.
The C. jejuni ssp. doylei MSPP scheme includes 14 different biomarker ions, mostly ribosomal proteins in the mass range of 2 to 11 kDa. MSPP can in principle, be adapted to other mass spectrometric platforms with an extended mass range. Therefore, this technique has the potential to become a useful tool for strain level microbial typing.

Subtyping of Campylobacter jejuni ssp. doylei Isolates Using Mass Spectrometry-based PhyloProteomics MSPP

Mass spectrometry-based phyloproteomics (MSPP) was used to type a collection of Campylobacter jejuni ssp. doylei isolates at the strain level in comparison to multilocus sequence typing (MLST).

MALDI-TOF MS offers the possibility to differentiate some bacteria not only at the species and subspecies level but even below, at the strain level. ...

Novel RNA-Binding Proteins Isolation by the RaPID Methodology

RNA-binding proteins (RBPs) play important roles in every aspect of RNA metabolism and regulation. Their identification is a major challenge in modern biology. Only a few in vitro and in vivo methods enable the identification of RBPs associated with a particular target mRNA. However, their main limitations are the identification of RBPs in a non-cellular environment (in vitro) or the low efficiency isolation of RNA of interest (in vivo). An RNA-binding protein purification and identification (RaPID) methodology was designed to overcome these limitations in yeast and enable efficient isolation of proteins that are associated in vivo. To achieve this, the RNA of interest is tagged with MS2 loops, and co-expressed with a fusion protein of an MS2-binding protein and a streptavidin-binding protein (SBP). Cells are then subjected to crosslinking and lysed, and complexes are isolated through streptavidin beads. The proteins that co-purify with the tagged RNA can then be determined by mass spectrometry. We recently used this protocol to identify novel proteins associated with the ER-associated PMP1 mRNA. Here, we provide a detailed protocol of RaPID, and discuss some of its limitations and advantages.

RNA-protein interactions lie at the heart of many cellular processes. Here, we describe an in vivo method to isolate specific RNA and identify novel proteins that are associated with it. This could shed new light on how RNAs are regulated in the cell.

RNA-binding proteins (RBPs) play important roles in every aspect of RNA metabolism and regulation. Their identification is a major challenge in modern ...

From a Natural Product to Its Biosynthetic Gene Cluster: A Demonstration Using Polyketomycin from Streptomyces diastatochromogenes T&#252;6028

Streptomyces strains are known for their capability to produce a lot of different compounds with various bioactivities. Cultivation under different conditions often leads to the production of new compounds. Therefore, production cultures of the strains are extracted with ethyl acetate and the crude extracts are analyzed by HPLC. Furthermore, the extracts are tested for their bioactivity by different assays. For structure elucidation the compound of interest is purified by a combination of different chromatography methods.
Genome sequencing coupled with genome mining allows the identification of a natural product biosynthetic gene cluster using different computer programs. To confirm that the correct gene cluster has been identified, gene inactivation experiments have to be performed. The resulting mutants are analyzed for the production of the particular natural product. Once the correct gene cluster has been inactivated, the strain should fail to produce the compound.
The workflow is shown for the antibacterial compound polyketomycin produced by Streptomyces diastatochromogenes T&#252;6028. Around ten years ago, when genome sequencing was still very expensive, the cloning and identification of a gene cluster was a very time-consuming process. Fast genome sequencing combined with genome mining accelerates the trial of cluster identification and opens up new ways to explore biosynthesis and to generate novel natural products by genetic methods. The protocol described in this paper can be assigned to any other compound derived from a Streptomyces strain or another microorganism.

From a Natural Product to Its Biosynthetic Gene Cluster: A Demonstration Using Polyketomycin from &lt;em&gt;Streptomyces diastatochromogenes&lt;/em&gt; T&amp;#252;6028

Here we present a detailed protocol of (A) the identification of a natural product with antibiotic activity, (B) the purification of the compound, (C) the first model of its biosynthesis, (D) genome sequencing/-mining and the (E) verification of the biosynthetic gene cluster.

Streptomyces strains are known for their capability to produce a lot of different compounds with various bioactivities. Cultivation under different ...

Determination of the Optimal Chromosomal Location(s) for a DNA Element in Escherichia coli Using a Novel Transposon-mediated Approach

The optimal chromosomal position(s) of a given DNA element was/were determined by transposon-mediated random insertion followed by fitness selection. In bacteria, the impact of the genetic context on the function of a genetic element can be difficult to assess. Several mechanisms, including topological effects, transcriptional interference from neighboring genes, and/or replication-associated gene dosage, may affect the function of a given genetic element. Here, we describe a method that permits the random integration of a DNA element into the chromosome of Escherichia coli and select the most favorable locations using a simple growth competition experiment. The method takes advantage of a well-described transposon-based system of random insertion, coupled with a selection of the fittest clone(s) by growth advantage, a procedure that is easily adjustable to experimental needs. The nature of the fittest clone(s) can be determined by whole-genome sequencing on a complex multi-clonal population or by easy gene walking for the rapid identification of selected clones. Here, the non-coding DNA region DARS2, which controls the initiation of chromosome replication in E. coli, was used as an example. The function of DARS2 is known to be affected by replication-associated gene dosage; the closer DARS2 gets to the origin of DNA replication, the more active it becomes. DARS2 was randomly inserted into the chromosome of a DARS2-deleted strain. The resultant clones containing individual insertions were pooled and competed against one another for hundreds of generations. Finally, the fittest clones were characterized and found to contain DARS2 inserted in close proximity to the original DARS2 location.

Determination of the Optimal Chromosomal Locations for a DNA Element in Escherichia coli Using a Novel Transposon-mediated Approach

Here, the power of a transposon-mediated random insertion of a non-coding DNA element was used to resolve its optimal chromosomal position.

The optimal chromosomal position(s) of a given DNA element was/were determined by transposon-mediated random insertion followed by fitness selection. In ...

Adaptation of Hybridization Capture of Chromatin-associated Proteins for Proteomics to Mammalian Cells

The hybridization capture of chromatin-associated proteins for proteomics (HyCCAPP) technology was initially developed to uncover novel DNA-protein interactions in yeast. It allows analysis of a target region of interest without the need for prior knowledge about likely proteins bound to the target region. This, in theory, allows HyCCAPP to be used to analyze any genomic region of interest, and it provides sufficient flexibility to work in different cell systems. This method is not meant to study binding sites of known transcription factors, a task better suited for Chromatin Immunoprecipitation (ChIP) and ChIP-like methods. The strength of HyCCAPP lies in its ability to explore DNA regions for which there is limited or no knowledge about the proteins bound to it. It can also be a convenient method to avoid biases (present in ChIP-like methods) introduced by protein-based chromatin enrichment using antibodies. Potentially, HyCCAPP can be a powerful tool to uncover truly novel DNA-protein interactions. To date, the technology has been predominantly applied to yeast cells or to high copy repeat sequences in mammalian cells. In order to become the powerful tool we envision, HyCCAPP approaches need to be optimized to efficiently capture single-copy loci in mammalian cells. Here, we present our adaptation of the initial yeast HyCCAPP capture protocol to human cell lines, and show that single-copy chromatin regions can be efficiently isolated with this modified protocol.

This is a method to identify novel DNA-interacting proteins at specific target loci, relying on sequence-specific capture of crosslinked chromatin for subsequent proteomic analyses. No prior knowledge about potential binding proteins, nor cell modifications are required. Initially developed for yeast, the technology has now been adapted for mammalian cells.

The hybridization capture of chromatin-associated proteins for proteomics (HyCCAPP) technology was initially developed to uncover novel DNA-protein ...

A Fast and Quantitative Method for Post-translational Modification and Variant Enabled Mapping of Peptides to Genomes

Cross-talk between genes, transcripts, and proteins is the key to cellular responses; hence, analysis of molecular levels as distinct entities is slowly being extended to integrative studies to enhance the understanding of molecular dynamics within cells. Current tools for the visualization and integration of proteomics with other omics datasets are inadequate for large-scale studies. Furthermore, they only capture basic sequence identify, discarding post-translational modifications and quantitation. To address these issues, we developed PoGo to map peptides with associated post-translational modifications and quantification to reference genome annotation. In addition, the tool was developed to enable the mapping of peptides identified from customized sequence databases incorporating single amino acid variants. While PoGo is a command line tool, the graphical interface PoGoGUI enables non-bioinformatics researchers to easily map peptides to 25 species supported by Ensembl genome annotation. The generated output borrows file formats from the genomics field and, therefore, visualization is supported in most genome browsers. For large-scale studies, PoGo is supported by TrackHubGenerator to create web-accessible repositories of data mapped to genomes that also enable an easy sharing of proteogenomics data. With little effort, this tool can map millions of peptides to reference genomes within only a few minutes, outperforming other available sequence-identity based tools. This protocol demonstrates the best approaches for proteogenomics mapping through PoGo with publicly available datasets of quantitative and phosphoproteomics, as well as large-scale studies.

Here we present the proteogenomic tool PoGo and protocols for fast, quantitative, post-translational modification and variant enabled mapping of peptides identified through mass spectrometry onto reference genomes. This tool is of use to integrate and visualize proteogenomic and personal proteomic studies interfacing with orthogonal genomics data.

Cross-talk between genes, transcripts, and proteins is the key to cellular responses; hence, analysis of molecular levels as distinct entities is slowly ...

Exploring the Effects of Spaceflight on Mouse Physiology using the Open Access NASA GeneLab Platform

Performing biological experiments in space requires special accommodations and procedures to ensure that these investigations are performed effectively and efficiently. Moreover, given the infrequency of these experiments it is imperative that their impacts be maximized. The rapid advancement of omics technologies offers an opportunity to dramatically increase the volume of data produced from precious spaceflight specimens. To capitalize on this, NASA has developed the GeneLab platform to provide unrestricted access to spaceflight omics data and encourage its widespread analysis. Rodents (both rats and mice) are common model organisms used by scientists to investigate space-related biological impacts. The enclosure that house rodents during spaceflight are called Rodent Habitats (formerly Animal Enclosure Modules), and are substantially different from standard vivarium cages in their dimensions, air flow, and access to water and food. In addition, due to environmental and atmospheric conditions on the International Space Station (ISS), animals are exposed to a higher CO2 concentration. We recently reported that mice in the Rodent Habitats experience large changes in their transcriptome irrespective of whether animals were on the ground or in space. Furthermore, these changes were consistent with a hypoxic response, potentially driven by higher CO2 concentrations. Here we describe how a typical rodent experiment is performed in space, how omics data from these experiments can be accessed through the GeneLab platform, and how to identify key factors in this data. Using this process, any individual can make critical discoveries that could change the design of future space missions and activities.

The NASA GeneLab platform provides unfettered access to precious omics data from biological spaceflight experiments. We describe how a typical mouse experiment is conducted in space and how data from such experiments can be accessed and analyzed.

Performing biological experiments in space requires special accommodations and procedures to ensure that these investigations are performed effectively ...

Using In Vitro and In-cell SHAPE to Investigate Small Molecule Induced Pre-mRNA Structural Changes

In the process of drug development of RNA-targeting small molecules, elucidating the structural changes upon their interactions with target RNA sequences is desired. We herein provide a detailed in vitro and in-cell selective 2&#39;-hydroxyl acylation analyzed by primer extension (SHAPE) protocol to study the RNA structural change in the presence of an experimental drug for spinal muscular atrophy (SMA), survival of motor neuron (SMN)-C2, and in exon 7 of the pre-mRNA of the SMN2 gene. In in vitro SHAPE, an RNA sequence of 140 nucleotides containing SMN2 exon 7 is transcribed by T7 RNA polymerase, folded in the presence of SMN-C2, and subsequently modified by a mild 2&#39;-OH acylation reagent, 2-methylnicotinic acid imidazolide (NAI). This 2&#39;-OH-NAI adduct is further probed by a 32P-labeled primer extension and resolved by polyacrylamide gel electrophoresis (PAGE). Conversely, 2&#39;-OH acylation in in-cell SHAPE takes place in situ with SMN-C2 bound cellular RNA in living cells. The pre-mRNA sequence of exon 7 in the SMN2 gene, along with SHAPE-induced mutations in the primer extension, was then amplified by PCR and subject to next-generation sequencing. Comparing the two methodologies, in vitro SHAPE is a more cost-effective method and does not require computational power to visualize results. However, the in vitro SHAPE-derived RNA model sometimes deviates from the secondary structure in a cellular context, likely due to the loss of all interactions with RNA-binding proteins. In-cell SHAPE does not need a radioactive material workplace and yields a more accurate RNA secondary structure in the cellular context. Furthermore, in-cell SHAPE is usually applicable for a larger range of RNA sequences (~1,000 nucleotides) by utilizing next-generation sequencing, compared to in vitro SHAPE (~200 nucleotides) that usually relies on PAGE analysis. In case of exon 7 in SMN2 pre-mRNA, the in vitro and in-cell SHAPE derived RNA models are similar to each other.

Detailed protocols for both in vitro and in-cell selective 2&#39;-hydroxyl acylation analyzed by primer extension (SHAPE) experiments to determine the secondary structure of pre-mRNA sequences of interest in the presence of an RNA-targeting small molecule are presented in this article.

In the process of drug development of RNA-targeting small molecules, elucidating the structural changes upon their interactions with target RNA sequences ...

Sequencing of mRNA from Whole Blood using Nanopore Sequencing

Sequencing in remote locations and resource-poor settings presents unique challenges. Nanopore sequencing can be successfully used under such conditions, and was deployed to West Africa during the recent Ebola virus epidemic, highlighting this possibility. In addition to its practical advantages (low cost, ease of equipment transport and use), this technology also provides fundamental advantages over second-generation sequencing approaches, particularly the very long read length, ability to directly sequence RNA, and real-time availability of data. Raw read accuracy is lower than with other sequencing platforms, which represents the main limitation of this technology; however, this can be partially mitigated by the high read depth generated. Here, we present a field-compatible protocol for sequencing of the mRNAs encoding for Niemann-Pick C1, which is the cellular receptor for ebolaviruses. This protocol encompasses extraction of RNA from animal blood samples, followed by RT-PCR for target enrichment, barcoding, library preparation, and the sequencing run itself, and can be easily adapted for use with other DNA or RNA targets.

Nanopore sequencing is a novel technology that allows cost-effective sequencing in remote locations and resource-poor settings. Here, we present a protocol for sequencing of mRNAs from whole blood that is compatible with such conditions.

Sequencing in remote locations and resource-poor settings presents unique challenges. Nanopore sequencing can be successfully used under such conditions, ...

A Complete Pipeline for Isolating and Sequencing MicroRNAs, and Analyzing Them Using Open Source Tools

Half of all human transcripts are thought to be regulated by microRNAs. Therefore, quantifying microRNA expression can reveal underlying mechanisms in disease states and provide therapeutic targets and biomarkers. Here, we detail how to accurately quantify microRNAs. Briefly, this method describes isolating microRNAs, ligating them to adaptors suitable for high-throughput sequencing, amplifying the final products, and preparing a sample library. Then, we explain how to align the obtained sequencing reads to microRNA hairpins, and quantify, normalize, and calculate their differential expression. Versatile and robust, this combined experimental workflow and bioinformatic analysis enables users to begin with tissue extraction and finish with microRNA quantification.

Here, we describe a step-by-step strategy for isolating small RNAs, enriching for microRNAs, and preparing samples for high-throughput sequencing. We then describe how to process sequence reads and align them to microRNAs, using open source tools.

Half of all human transcripts are thought to be regulated by microRNAs. Therefore, quantifying microRNA expression can reveal underlying mechanisms in ...