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.
posted
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.
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 Scienti c) for swabbing and sample extraction
CAUTION: ethanol is highly ammable and irritating to eyes and skin.
Water of LC-MS (Optima) grade (ThermoFisher Scienti c).
ESI-L Low Concentration Tuning Mix (Agilent Technologies, CA, USA) for external calibration
of the MaXis II QTOF mass spectrometer
Hexakis(1H,1H,3H-tetra uoropropoxy)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 Scienti c). CAUTION: Acetonitrile is highly
ammable.
Formic acid of LC-MS grade, Optima grade (ThermoFisher Scienti c). 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 Scienti c, MA,
USA) reference 278753 or an equivalent generic 96 wells microtiter plate, 400 µL, polystyrene.
Reusable lid for a 96 wells plate (ThermoFisher Scienti c), 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 Scienti c, 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.
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.
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 certi ed shipping service.
Metabolite extraction of skin samples
temperature and longer extraction time. However, this can lead to the degradation of some
molecules.
this step prevent contamination of other wells by covering them with aluminium lm.
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.
swabs and extraction solvent onto each plate.
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 speci c 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.
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.
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.
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.
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%.
samples order should be randomized.
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.
(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.
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).
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 speci cs 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 fi les are loaded and batch-processed with the following settings for each step in the batch:
Mass detection
Noise level 1000
Minimum peak height 3000
m/z tolerance 0.1 m/z or 20 ppm
m/z tolerance 0.1 m/z or 20 ppm
Minimum peak in a row 3
based on speci cs 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.
the nodes appeared in each other's respective top 10 most similar nodes.
level 2/3 according to the proposed minimum standards in metabolomics[7].
the status page.
Open the Cytoscape software and import the downloaded file: File ! Import ! Network ! File. . . and then select the .graphml fi le 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.
(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 fi le 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.
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.
"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)
Under the "Compounds" tab (left hand side), the currently active job is selected.
displayed below for each formula or could be viewed individually in the corresponding tabs.
experimental spectrum is not present in selected source databases.
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 speci c 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 speci c 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 speci c to the described work ow 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 pro le 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
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.
posted
You are reading this latest protocol version