Data generation and analysis with SIRIUS 4 on two biological case studies

doi:10.1038/protex.2018.133

Method Article

Data generation and analysis with SIRIUS 4 on two biological case studies

https://doi.org/10.1038/protex.2018.133

This work is licensed under a CC BY 4.0 License

This protocol has been posted on Protocol Exchange, an open repository of community-contributed protocols sponsored by Nature Portfolio. These protocols are posted directly on the Protocol Exchange by authors and are made freely available to the scientific community for use and comment.

Version 1

posted 20 Sep, 2019

You are reading this latest protocol version

The procedures for the metabolomics analysis performed in the present work has been described previously for skin[1] and fecal material[2]. The detailed step by step instructions to reproduce the biological case studies presented in this work are given below.

Chemistry

Biochemistry

Biological techniques

mass spectrometry

sirius

in-silico ms-ms

This protocol summarizes the workflow on how to get from sample to the results described in the paper Sirius 4

Ethanol of LC-MS (Optima) grade (ThermoFisher Scientic) for swabbing and sample extraction

CAUTION: ethanol is highly ammable and irritating to eyes and skin.

Water of LC-MS (Optima) grade (ThermoFisher Scientic).

ESI-L Low Concentration Tuning Mix (Agilent Technologies, CA, USA) for external calibration

of the MaXis II QTOF mass spectrometer

Hexakis(1H,1H,3H-tetrauoropropoxy)phosphazene (Synquest Laboratories, FL, USA), m/z

922.009798, for internal calibration (lockmass) of the MaXis II QTOF mass spectrometer.

CAUTION: This compound is irritating to eyes and skin.

Acetonitrile LC-MS (Optima) grade (ThermoFisher Scientic). CAUTION: Acetonitrile is highly

ammable.

Formic acid of LC-MS grade, Optima grade (ThermoFisher Scientic). CAUTION: Formic acid

is highly corrosive.

(Skin sampling): Cotton swabs, model Puritan® (Spectrum, CA, USA), reference 150-533353

or 150-2971696.

(Faeces sampling): BBL culture swabs (Becton, Dickinson and Company, Sparks, MD).

Deep-wells microtiter plate, 2 mL, polypropylene, model NUNC (ThermoFisher Scientic, MA,

USA) reference 278753 or an equivalent generic 96 wells microtiter plate, 400 µL, polystyrene.

Reusable lid for a 96 wells plate (ThermoFisher Scientic), reference 351190 or equivalent.

Stainless steel tweezers.

Single or multi channel 20-200 µL micropipette.

Plate mat for storage, model Corning Storage Mat III (ThermoFisher Scientic, MA, USA).

Benchtop vacuum concentrator compatible with 96-well plates evaporation (Centrivap Labconco,

MO, USA).

Reverse phase C18 LC column 1.7 µm particle size, 50 2.1 mm (Phenomenex, CA, USA) or

equivalent.

Ultra high-performance liquid chromatography (UHPLC) system coupled to a tandem mass

spectrometer with an ESI source; for the biological studies described here we have used UltiMate

3000 UHPLC (Dionex, Idstein, Germany) coupled to MaXis II QTOF (Bruker Daltonics,

Bremen, Germany).

Sampling of skin

Prepare the swabs by disposing of the ethanol:water 1:1 (v/v) from the swab container and by replacing it with fresh ethanol:water 1:1 (v/v).
Prepare a 96-well plate (volume 2 mL) by adding 500 µL of the solvent mix into each well, and placing it on an ice bed.
For each sampling spot:
Take one swab, and swab vigorously the spot of interest with the cotton bud side of the swab for at least five seconds. CRITICAL STEP: The swabbing procedure must be as consistent

as possible across spots.

Place the swab into the corresponding well of the 96-well plate and cut it with a pair of scissors in order to leave only the cotton bud in the well. Dispose the swab wood stick.

CRITICAL STEP: During the placing of the swab and the cutting of the cotton bud, use

aluminium fold lm to cover the rest of the plate to prevent cross-contamination of other

wells.

Close the 96 well plate with a polypropylene reusable cover. CAUTION: All samples must be prepared using the same extraction protocol, using the same

reagents, and ideally by the same operator.

Sampling of stool

The sampling is carried out by individuals as described in [2]. Samples are

collected using swabs and returned by mail. Samples collected outside of the United States are

shipped using domestic post within each country to an aggregation site and stored at -80 at the

aggregation site until shipment to the United States. Shipment to the United States is done on dry

ice using a certied shipping service.

Metabolite extraction of skin samples

Leave the 96-well plate with soaked cotton buds at a temperature of 4 for 24 to 48 hours. CRITICAL STEP: Note that the extraction yield can be increased with higher extraction

temperature and longer extraction time. However, this can lead to the degradation of some

molecules.

Remove each cotton bud with a pair of tweezers, and retrieve as much solvent as possible by pressing the bud on the side of the well. Dispose of the cotton bud. CRITICAL STEP: During

this step prevent contamination of other wells by covering them with aluminium lm.

Evaporate the 96-well plate solvent content using the benchtop vacuum centrifuge concentrator. CAUTION: The use of centrifugation device requires specic training, in particular to ensure

that the rotor is properly equilibrated before spinning.

Redissolve the plate content in 120 µL of 1:1 (v/v) acetonitrile:water at ambient temperature.
Cover the plate with a lid. Agitate gently the plate and use ultrasonic bath for 5 min to ensure the appropriate dissolution of the sample.
Centrifugate the 96-well plate at 2000 rpm for 15 minutes at ambient temperature.
Transfer 100 µL of the plate content with multi-channel pipette into a new 96-well plate (well volume 400 µL). CRITICAL STEP: Make sure to reproduce the order of the wells and avoid cross

contamination of the plate by covering it with aluminium foil, and avoid undissolved material

upon sample transfer.

(Optional, but recommended). Pool 10 µL from each of twelve representative samples to create a pooled QC sample.
Store the plate with extracts at -20 until further analysis. Note that for longer storage, storing the samples at -80 is recommended, followed by sonication of extracts to redissolve

any precipitated material.

Metabolite extraction of stool samples

Remove the swab tubes scheduled for analysis from the -80 freezer and place on dry ice for the duration of sample processing.
Place each swab onto a Thermo Fisher Scientific (Waltham, MA) 2-ml deep-well 96-well plate set on top of dry ice coolant. CRITICAL STEP: Make sure to reproduce the order of the wells and

avoid cross contamination of the plate by covering it with aluminium foil, and avoid undissolved

material upon sample transfer.

Snap the top part of each swab stick o and discard.
Add 200 µl of HPLC-grade 90% (vol/vol) ethanol-water solvent to each well using a multichannel pipette immediately after lling all wells. It is recommended to include four blanks of unused

swabs and extraction solvent onto each plate.

Seal with a 96-well plate lid, sonicate for 10 min
Place into the refrigerator at 2 to extract samples overnight.
Remove and discard swabs.
Place the plates into a lyophilizer, and dry down the entire sample.
Add 200 µl 90% (vol/vol) ethanol-water for resuspension.
Resealed the plates and centrifuge at 2,000 rpm for 10 min. CAUTION: The use of centrifugation device requires specic training, in particular to ensure that the rotor is properly equilibrated

before spinning.

Transfer 100-µl aliquots of sample onto a Falcon 96-well MS plate using a multichannel pipette, seal plate immediately with sealing lm.
Centrifuge sealed plates at 2,000 rpm for 10 min and analyze immediately or store at 2 until analysis. CAUTION: The use of centrifugation device requires specic training, in particular to

ensure that the rotor is properly equilibrated before spinning.

Mass spectrometry analysis

Prepare the LC mobile phase A of 100:0.1 (v/v) water:formic acid.
Prepare the LC mobile phase B of 100:0.1 (v/v) acetonitrile:formic acid.
Prime the C18 HPLC column according to the manufacturer's guidelines.
Set up the LC linear gradient method as follows: 0-0.5 min 2% B, 0.5-2 min 2-20% B, 2-8 min 20-98% B, 8-9 min 98-98% B, and 9-10 min 2% B with a constant of flow rate of 500 µL/min.
Set up the LC parameters as following: column temperature (40), loop factor wash (3) (Dionex UltiMate 3000, or as appropriate for other systems), injection volume (10 µL for skin samples,

5 µL for fecal samples). Prime the injection system multiple time.

(Optional) Clean-up the ESI source of the MS instrument according to the manufacturer's guidelines. Calibrate instruments using the calibration solution according to manufacturer's

guideline. For the MaXis II QTOF, add the solution containing internal lockmass (m/z 922.009)

in the ESI source in the positive ion mode, following the the manufacturer's guideline.

Set up the MS method as follows: for the MaXis II QTOF MS, acquire spectra in positive ion mode in the mass range of m/z 802,000 with the following settings: capillary voltage, 4,500

V; ion source temperature, 180 ; dry gas ow, 9 L/min; spectra rate acquisition, 3 spectra/s.

Perform MS/MS fragmentation of the seven most intense selected ions per spectrum (TOP7)

using ramped collision-induced dissociation (CID) energy, ranging from 16 to 48 eV, to get

diverse fragmentation patterns. Set MS/MS exclusion after 3 spectra to be released after 30 s.

Set an MS/MS exclusion list for the mass range of m/z 921.5924.5 to exclude the lockmass, and

if needed include other adducts of the lockmass if their intensities are in the range of MS/MS

threshold.

Run a blank sample and verify that no signicant contamination from source is impacting the system by assessing abundances of observed ions and ensuring that they do not exceed

thresholds appropriate for the instrument (e.g. below few thousand counts for the MaXis II

QTOF instrument). If high abundance contaminants are present, consider additional source

cleaning, replacement of mobile phases, glassware, LC column or tubing, as needed.

(Optional, but recommended): Start by running the QC sample, followed by a blank, and repeat it at least two times. Check for satisfactory reproducibility by comparing multiple QC

samples (coefficient of variation, CV, below 0.1-0.3). Check for carryover by inspecting the blank

that followed a QC sample. If needed, adjust the injection volume to prevent chromatographic

saturation (peak broadening), and mass spectrometer detector saturation (peak flattening). If

needed, optimize the LC gradient to ensure appropriate separation of representative features,

and the number of needle/injector washes to reduce carry-over below 1%.

Prepare the LC-MS sequence for the samples injections. Include the analysis of the blank and the QC sample every 12 samples or less, depending on total number of samples. If possible, the

samples order should be randomized.

(Optional): Include a clean-up sample every 12 samples to limit column saturation. This is recommended if sample contains components that are strongly retained by the column and are

not eluted completely by the mobile phase. Typically, this clean-up run is composed of multiple

mobile phase composition variation according to the manufacturer's guidelines.

Run the LC-MS sequence, and monitor regularly that (i) elution prole of every QC sample is reproducible (5-10 sec window maximum for UHPLC, or a coefficient of variation of 0.1-0.3), and

(ii) that the calibration is stable (< 10 ppm for MaXis II QTOF). If needed, stop the LC-MS

sequence, and calibrate the instrument according to manufacturer guidelines. For the MaXis II

QTOF, lock mass calibrant must be added every 12 hours. If the LC-MS sequence was stopped,

run a QC sample and a blank prior continuing the rest of your analysis.

After nishing the LC-MS experiments, open representative les, from both the beginning and the end of the sequence, with the MS data viewing software (Data Analysis for MaXis II QTOF).

Overlap features across samples, such as the internal standard, and assess the stability of the

chromatography (retention shift in seconds), and the MS calibration (m/z window in ppm).

Export all LC-MS/MS les into the centroided .mzXML format. For the MaXis II QTOF MS, use CompassXport to apply calibration, and convert the data from .d format to .mzXML format

in centroid mode (select Filters: Peak Picking, MS-Levels 1-2).

Data analysis

Creating feature tables. Multiple software packages for MS feature extraction exist and can be used;

the present protocol is given for the use of the open-source OpenMS 2.0[3] software utilized for feature

detection in the fecal sampling biological study example and MZmine2[4] used in the skin sampling

biological study. The recommended settings were found to be appropriate for the data obtained in

our experiments; however, experienced users may nd changes in parameters necessary, especially

if dierent instrumentation is used.

Creating feature tables with OpenMS.

Prior to feature extraction, the collected HPLC-MS raw data les are converted from Bruker's .d to .mzXML format
Peaks across dierent chromatograms deconvolved and aligned (feature detection). The alignment window is set at 0.5 min, the noise threshold is set at 1,000 counts, the

chromatographic peak full width at half-maximum (FWHM) value is set at 20, and the mass

error is set at 30 ppm.

An output table with detected MS features across all samples and their corresponding abundances is obtained.
(Recommended): All of the peaks that are present in any of the blanks with a signal-to-noise ratio (S/N) below 10:1 are removed from the final feature table. The threshold can be adjusted

based on specics of each experiment.

Creating feature tables with MZmine2.

Prior to feature extraction, the collected HPLC-MS raw data files are converted from Bruker's .d to .mzXML format.
The experimental files are loaded and batch-processed with the following settings for each step in the batch:

Mass detection

Noise level 1000
- Chromatogram builder
- Minimum time span 0.01 min
- Minimum peak height 3000
- m/z tolerance 0.1 m/z or 20 ppm _Chromatogram deconvolution - Baseline cutoff
Minimum peak height 3000
- Peak duration range (min) 0.01 - 3.00
- Baseline level 300 Deisotopisation - Isotopic peak grouper
m/z tolerance 0.1 m/z or 20 ppm
- RT tolerance 0.1 min
- Maximum charge 4 Peak alignment - Join aligner
m/z tolerance 0.1 m/z or 20 ppm
- Weight for m/z 75
- Weight for RT 25
- RT tolerance (min) 0.1 Peak ltering - Peak list raw lter
Minimum peak in a row 3
- Minimum peak in an isotope pattern 2
1. An output table with detected MS features across all samples and their corresponding abundances is obtained.

(Recommended): All of the peaks that are present in any of the blanks with a signal-to-noise ratio (S/N) below 10:1 are removed from the nal feature table. The threshold can be adjusted

based on specics of each experiment.

Molecular networking

Raw mass spectrometry data files in the .mzXML format are uploaded MassIVE mass spectrometry database (https://massive.ucsd.edu/).
Molecular networking[5] is performed to identify spectra shared between dierent sample types and to identify known molecules in the data set.
MS/MS spectra are window ltered by choosing only the top 6 peaks in the 50-Da window throughout the spectrum.
The MS spectra were are clustered with a parent mass tolerance of 0.02 Da and an MS/MS fragment ion tolerance of 0.02; consensus spectra that contained fewer than 4 spectra are

discarded.

A network is created where with edges ltered to have a cosine score above 0.65 and more than 5 matched peaks. The edges between two nodes are kept in the network if and only if each of

the nodes appeared in each other's respective top 10 most similar nodes.

Set required library matches to have a score above 0.7 and at least 6 matched peaks, and search the spectra in the network against GNPS spectral libraries[6]. All resulting annotations are at

level 2/3 according to the proposed minimum standards in metabolomics[7].

Molecular networks are visualized and using Cytoscape software[8]: download the required files from the generated GNPS job by clicking on the Download GraphML for Cytoscape link on

the status page.

Open the Cytoscape software and import the downloaded file: File ! Import ! Network ! File. . . and then select the .graphml file in the downloaded folder.
Select appropriate labels for the network nodes and edges under the Style tab.
The structure of molecules interest that are found within the cluster with existing annotations may be postulated using annotation propagation from adjacent annotated nodes in the cluster

as described in reference [6] by assessing dierences in parent mass and fragmentation patterns.

(Optional, if possible): For level 1 identications, authentic standards for postulated compounds need to be procured and analyzed under identical experimental conditions. The retention times

(RTs) and MS/MS spectra need to be compared to those experimentally obtained for the

putatively annotated compounds.

Determining the molecular structure with SIRIUS 4

Extract the MS/MS spectrum of the compound of interest; this can be achieved in several ways: Extract the MS/MS using MSConvert software (ProteoWizard package)9:
Launch MSConvert GUI
Click "Browse" for the "File" field to select the experimental .mzXML file that contains MS/MS spectrum of interest.
Click "Add" to include the file into the analysis window.
Click "Browse" for the "Output Directory" field to select the output directory.
Select "Output format": mgf
Select "Binary encoding precision": 32 bit.
Uncheck "Use zlib compression" box
Under "Filters", select "Subset" and enter the scan number of the MS/MS spectrum of interest.
Click "Add".
Click "Start".
The .mgf file containing the single MS/MS will be saved in the selected output folder. Extract the MS/MS peak list directly from GNPS job:
When viewing the in-browser network cluster, select the node of interest by left clicking on it. The top right window will display the MS/MS of the consensus spectrum for this

node.

Click on the "Download MS/MS Peaks" icon in the top right corner of the MS/MS spectrum window. The browser will download the .txt file. CAUTION: The downloaded

MS/MS peaks will belong to the consensus spectrum obtained from clustering multiple

experimental spectra and thus may dier from MS/MS patterns in individual files.

Caution must be used to ensure that the appropriate spectrum is used for further analysis.

Launch the GUI of the SIRIUS software.
Import the MS/MS spectrum (spectra) for the analysis:

Click "Import" icon
Click "+" icon
Select the file with MS/MS of interest (enter the parent m/z if not recognized from the file)
Alternatively: Click "Batch Import" icon
Select the file with MS/MS of interest (Note: for batch import, multiple MS/MS spectra can be loaded at once from a single file. The .txt files cannot be used for batch import).

For single spectrum import, select the ionization mode ([M+H]+ for the present case studies).
Click "Compute All"
In the pop-up window, select:

"Elements beside CHNOP": include sulfur for the present biological case studies; for the range of expected number of atoms use up to "+2" instead of "auto"
Instrument: "Q-TOF" for the present biological case studies
"ppm": "20" for the present biological case studies
"candidates": default of 10
In the drop down menu, select "formulas from biological databases" (in the present case studies) or other option if appropriate (Note: selecting "all possible molecular formulas"

would result in longer computing time)

Click the "CSI:FingerID" button
- Under "Search in" select "biological database"
- Select possible adducts ([M+H]+ is selected by default)
1. Click "Compute".
2. Once the "Jobs" icon is static, it indicates that the computation is completed and the results can be explored.

Under the "Compounds" tab (left hand side), the currently active job is selected.

In the "Sirius Overview" tab (main window), browse the results by clicking the suggested formulas listed on top of the window. The explained fragments and calculated fragmentation tree are

displayed below for each formula or could be viewed individually in the corresponding tabs.

Select "CSI:FingerID Details" tab for possible structures for each formula. Specic libraries for possible candidate molecules could be selected. If no entries are displayed, the molecule with the

experimental spectrum is not present in selected source databases.

In both of the biological case studies, the correct structures were found to have the highest scores for both the molecular formulas and predicted structures; however, oftentimes that may

not be the case and researcher's discretion and expertise are essential for selecting the potential

candidate molecules for further consideration and interpretation.

The sampling time for skin swabbing and faeces sampling is generally under one minute per sample.

The total time necessary for the samples preparation and data analysis depends on the number of

samples. Typical timing is approximately one hour for placing swabs into 96-well plate, 6 hours to

12 hours for metabolite extraction, 24 hours for LC-MS analysis, 2 hours for MS feature extraction

(strongly depends on the computational power of the system that runs the software), 6 hours

for preparing and running GNPS job. The time for network-based annotation propagation can vary

widely depending on the analysis circumstances. The time for a single SIRIUS/CSI:FingerID job for

a single spectrum processing is typically few seconds on a personal computer. A detailed breakdown

of timing for specic experimental steps is given in publication of Sirius 4.

The users should follow general good laboratory practices for the analysis to avoid issues such

as contamination. For detailed troubleshooting guidance specic to the LC-MS/MS analysis, the

user should refer to the manuals of the instrumentation that is used in the experiment. Additional

troubleshooting guidelines specic to the described workow are given in [1]. The troubleshooting

tips for the GNPS-based molecular networking and annotation could be accessed at: https://

ccms-ucsd.github.io/GNPSDocumentation/troubleshooting/

The anticipated results should be comparable to those described in the given examples for the two

biological case studies.

Protsyuk, I. et al. 3D molecular cartography using LC-MS facilitated by Optimus and 'ili software. Nature protocols 13, 134154 (2018).
McDonald, D. et al. American Gut: an Open Platform for Citizen Science Microbiome Research. mSystems 3 (2018).
Röst, H. L. et al. OpenMS: a exible open-source software platform for mass spectrometry data analysis. Nat Methods 13, 741748 (2016).
Pluskal, T., Castillo, S., Villar-Briones, A. & Oresic, M. MZmine 2: Modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular prole data. BMC

Bioinf 11, 395 (2010).

Watrous, J. et al. Mass spectral molecular networking of living microbial colonies. Proc Natl Acad Sci U S A 109, E1743E1752 (2012).
Wang, M. et al. Sharing and community curation of mass spectrometry data with Global Natural Products Social Molecular Networking. Nat Biotechnol 34, 828837 (2016).
Sumner, L. W. et al. Proposed minimum reporting standards for chemical analysis. Metabolomics 3, 211221 (2007).
Shannon, P. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 13, 24982504 (2003).
Adusumilli, R. & Mallick, P. Data Conversion with ProteoWizard msConvert. Methods in molecular biology (Clifton, N.J.) 1550, 339368 (2017).

No conflicting financial interests

PDF

Version 1

posted 20 Sep, 2019

You are reading this latest protocol version

Data generation and analysis with SIRIUS 4 on two biological case studies

Status:

Version 1

Abstract

Introduction

Reagents

Equipment

Procedure

Timing

Troubleshooting

Anticipated Results

References

Additional Declarations

Associated Publications

Status:

Version 1

Privacy Policy

Terms of Service