nf-core/quantms

URI/path to an SDRF file (.sdrf.tsv) OR OpenMS-style experimental design with paths to spectra files (.tsv)

The output directory where the results will be saved.

Email address for completion summary.

MultiQC report title. Printed as page header, used for filename if not otherwise specified.

Root folder in which the spectrum files specified in the SDRF/design are searched

Overwrite the file type/extension of the filename as specified in the SDRF/design

Proteomics data acquisition method

The fasta protein database used during database search. Note: For DIA data, it must not contain decoys.

Generate and append decoys to the given protein database

Pre- or suffix of decoy proteins in their accession

Location of the decoy marker string in the fasta accession. Before (prefix) or after (suffix)

Choose the method to produce decoys from input target database.

Maximum nr. of attempts to lower the amino acid sequence identity between target and decoy for the shuffle algorithm

Target-decoy amino acid sequence identity threshold for the shuffle algorithm. if the sequence identity is above this threshold, shuffling is repeated. In case of repeated failure, individual amino acids are ‘mutated’ to produce a difference amino acid sequence.

Activate OpenMS-internal peak picking

Perform peakpicking in memory

Which MS levels to pick as comma separated list. Leave empty for auto-detection.

Convert bruker .d files to mzML

Force initial re-indexing of input mzML files. Also fixes some common mistakes in slightly incomplete/outdated mzMLs. (Default: true for safety)

A comma separated list of search engines to use (and combine). Valid: comet, msgf, sage

Number of sage processes to be spawned.

The enzyme to be used for in-silico digestion, in ‘OpenMS format’

Specify the amount of termini matching the enzyme cutting rules for a peptide to be considered. Valid values are fully (default), semi, or none

Specify the maximum number of allowed missed enzyme cleavages in a peptide. The parameter is not applied if unspecific cleavage is specified as enzyme.

Precursor mass tolerance used for database search. For High-Resolution instruments a precursor mass tolerance value of 5 ppm is recommended (i.e. 5). See also --precursor_mass_tolerance_unit.

Precursor mass tolerance unit used for database search. Possible values are ‘ppm’ (default) and ‘Da’.

Fragment mass tolerance used for database search. The default of 0.03 Da is for high-resolution instruments.

Fragment mass tolerance unit used for database search. Possible values are ‘ppm’ (default) and ‘Da’.

A comma-separated list of fixed modifications with their Unimod name to be searched during database search

A comma-separated list of variable modifications with their Unimod name to be searched during database search

Comma-separated range of integers with allowed isotope peak errors for precursor tolerance (like MS-GF+ parameter ‘-ti’). E.g. -1,3

Type of instrument that generated the data. ‘low_res’ or ‘high_res’ (default; refers to LCQ and LTQ instruments)

MSGF only: Labeling or enrichment protocol used, if any. Default: automatic

Minimum precursor ion charge. Omit the ’+’

Maximum precursor ion charge. Omit the ’+’

Minimum peptide length to consider (works with MSGF and in newer Comet versions)

Maximum peptide length to consider (works with MSGF and in newer Comet versions)

Specify the maximum number of top peptide candidates per spectrum to be reported by the search engine. Default: 1

Maximum number of modifications per peptide. If this value is large, the search may take very long.

The minimum precursor m/z for the in silico library generation or library-free search

The maximum precursor m/z for the in silico library generation or library-free search

The minimum fragment m/z for the in silico library generation or library-free search

The maximum fragment m/z for the in silico library generation or library-free search

Turn the mechanism on.

Which variable modifications to use for scoring their localization.

What to do when peptides are found that do not follow a unified set of rules (since search engines sometimes differ in their interpretation of them).

Should isoleucine and leucine be treated interchangeably when mapping search engine hits to the database? Default: true

How to calculate posterior probabilities for PSMs:

• ‘percolator’ = Re-score based on PSM-feature-based SVM and transform distance to hyperplane for posteriors
• ‘fit_distributions’ = Fit positive and negative distributions to scores (similar to PeptideProphet)

FDR cutoff on PSM level (or peptide level; see Percolator options) per run before going into feature finding, map alignment and inference. This can be seen as a pre-filter. See

Calculate FDR on PSM (‘psm-level-fdrs’) or peptide level (‘peptide-level-fdrs’)?

The FDR cutoff to be used during training of the SVM.

The FDR cutoff to be used during testing of the SVM.

Only train an SVM on a subset of PSMs, and use the resulting score vector to evaluate the other PSMs. Recommended when analyzing huge numbers (>1 million) of PSMs. When set to 0, all PSMs are used for training as normal. This is a runtime vs. quality tradeoff. Default: 300,000

Use additional features whose values are learnt by correct entries. See help text. Default: 0 = none

How to handle outliers during fitting:

• ignore_iqr_outliers (default): ignore outliers outside of 3*IQR from Q1/Q3 for fitting
• set_iqr_to_closest_valid: set IQR-based outliers to the last valid value for fitting
• ignore_extreme_percentiles: ignore everything outside 99th and 1st percentile (also removes equal values like potential censored max values in XTandem)
• none: do nothing

How to combine the probabilities from the single search engines: best, combine using a sequence similarity-matrix (PEPMatrix), combine using shared ion count of peptides (PEPIons). See help for further info.

Only use the top N hits per search engine and spectrum for combination. Default: 0 = all

A threshold for the ratio of occurrence/similarity scores of a peptide in other runs, to be reported. See help.

Operate only on MSn scans where any of its precursors features a certain activation method. Set to empty to disable.

Allowed shift (left to right) in Th from the expected position

Minimum intensity of the precursor to be extracted

Minimum fraction of the total intensity. 0.0:1.0

Minimum intensity of the individual reporter ions to be extracted.

Maximum allowed deviation (in ppm) between theoretical and observed isotopic peaks of the precursor peak

Enable isotope correction (highly recommended)

Enable normalization of the channel intensities

The reference channel, e.g. for calculating ratios.

The inference method to use. ‘aggregation’ (default) or ‘bayesian’.

[Ignored in Bayesian] How to aggregate scores of peptides matching to the same protein

[Ignored in Bayesian] Also use shared peptides during score aggregation to protein level

[Ignored in Bayesian] Minimum number of peptides needed for a protein identification

Consider only the top X PSMs per spectrum to find the best PSM per peptide. 0 considers all.

The experiment-wide protein (group)-level FDR cutoff. Default: 0.01

Use picked protein FDRs

The experiment-wide PSM-level FDR cutoff. Default: 0.01

Specify the labelling method that was used. Will be ignored if SDRF was given but is mandatory otherwise

Calculate protein abundance from this number of proteotypic peptides (most abundant first; ‘0’ for all, Default 3)

Averaging method used to compute protein abundances from peptide abundances.

Distinguish between fraction and charge states of a peptide. (default: ‘false’)

Add the log2 ratios of the abundance values to the output.

Scale peptide abundances so that medians of all samples are equal.(Default false)

Use the same peptides for protein quantification across all samples.(Default false)

Include results for proteins with fewer proteotypic peptide than indicated by top.

Quantify proteins based on:

• ‘unique_peptides’ = use peptides mapping to single proteins or a group of indistinguishable proteins (according to the set of experimentally identified peptides)
• ‘strictly_unique_peptides’ (only LFQ) = use peptides mapping to a unique single protein only
• ‘shared_peptides’ = use shared peptides, too, but only greedily for its best group (by inference score and nr. of peptides)

Export the results in mzTab format.

Choose between feature-based quantification based on integrated MS1 signals (‘feature_intensity’; default) or spectral counting of PSMs (‘spectral_counting’). WARNING: ‘spectral_counting’ is not compatible with our MSstats step yet. MSstats will therefore be disabled automatically with that choice.

Recalibrates masses based on precursor mass deviations to correct for instrument biases. (default: ‘false’)

Only looks for quantifiable features at locations with an identified spectrum. Set to false to include unidentified features so they can be linked and matched to identified ones (= match between runs). (default: ‘true’)

The minimum probability (e.g.: 0.25) an identified (=id targeted) feature must have to be kept for alignment and linking (0=no filter).

The minimum probability (e.g.: 0.75) an unidentified feature must have to be kept for alignment and linking (0=no filter).

The minimum intensity for a feature to be considered for quantification. (default: ‘10000’)

The order in which maps are aligned. Star = all vs. the reference with most IDs (default). TreeGuided = an alignment tree is calculated first based on similarity measures of the IDs in the maps.

Also quantify decoys? (Usually only needed for Triqler post-processing output with --add_triqler_output, where it is auto-enabled)

Choosing the MS2 mass accuracy setting automatically

Choosing scan_window setting automatically

Set the scan window radius to a specific value

Only peaks with correlation sum exceeding min_corr will be considered

Peaks with correlation sum below corr_diff from maximum will not be considered

A single score will be used until RT alignment to save memory

Controls the protein inference mode

Instructs DIA-NN to add the organism identifier to the gene names

The spectral library to use for DIA-NN

Enable cross-run normalization between runs by diann.

Skip MSstats/MSstatsTMT for statistical post-processing?

Experimental: Instead of all pairwise contrasts (default), uses the given condition name/number (corresponding to your experimental design) as a reference and creates pairwise contrasts against it.

Experimental: Allows full control over contrasts by specifying a set of contrasts in a semicolon separated list of R-compatible contrasts with the condition names/numbers as variables (e.g. 1-2;1-3;2-3). Overwrites --ref_condition.

The threshold value for differential expressed proteins in MSstats plots based on adjusted p-value

Also create an output in Triqler’s format for an alternative manual post-processing with that tool

Which features to use for quantification per protein: ‘top3’ or ‘highQuality’ which removes outliers only

which summary method to use: ‘TMP’ (Tukey’s median polish) or ‘linear’ (linear mixed model)

Omit proteins with only one quantified feature?

Keep features with only one or two measurements across runs?

Use unique peptide for each protein

Remove the features that have 1 or 2 measurements within each run

select the feature with the largest summmation or maximal value

summarization methods to protein-level can be perfomed

Reference channel based normalization between MS runs on protein level data?

Remove ‘Norm’ channels from protein level data

Reference channel based normalization between MS runs on protein level data

Export MSstats profile QC plots including all proteins

Enable generation of pmultiqc report? default: ‘false’

Skip idXML files (do not generate search engine scores) in pmultiqc report? default: ‘true’

Skip protein/peptide table plots with pmultiqc for large dataset.

Custom MultiQC yaml file containing HTML including a methods description.