2 Functional Analysis
2 Functional Analysis
0
A Web-based Tool for streamlined
metabolomics data analysis
2022.07.12
2. Functional Analysis
The Functional Analysis module of MetaboAnalyst has undergone several major updates since it’s introduction in Version 4. First,
it includes a modified Gene Set Enrichment Analysis method, which considers the overall ranks of uploaded peaks and is capable
of detecting more subtle and consistent changes than the original mummichog algorithm (Li et al. 2013). Second, it supports the
inclusion of retention time when performing functional analysis to increase the confidence and robustness of putative compound
annotation. Finally, MetaboAnalyst 5.0 has included an interactive heatmap visualization of a user’s peak intensity table to help
users perform functional interpretation of manually identified patterns of interest.
Other Highlights:
- Users can upload either a peak intensity table (generic or MZMine) or peak list.
- Added support for pathway analysis of 26 organisms including human, mouse, zebrafish, C. elegans, among other species.
- Added ~9, 000 metabolite sets (e.g. Disease-associated sets, chemical classes) to be used for functional interpretation.
2.0 Knowledge & Background
● Mass spectrometry based untargeted metabolomics traditionally require metabolites to be identified before any biological meaning
can be drawn from the data. Metabolite identification is a challenging and low throughput process, therefore becomes the
bottleneck of the filed. Li et al. report here a novel approach to predict biological activity directly from mass spectrometry data
without a priori identification of metabolites by unifying network analysis and metabolite prediction under the same computational
framework. (version 1)
● The algorithm has been further enhanced to version 2 by considering the retention time information for more accuracy by
introducing empirical compounds. Empirical Compounds are intermediaries between m/z features and compounds. The steps for
how they are formed are as follows:
First, all m/z features are matched to potential compounds considering different adducts. Then, per compound, all matching m/z features are split into
Empirical Compounds based on whether they match within an expected retention time window. The retention time window (in seconds) is calculated
as the maximum retention time * 0.02. This results in the initial Empirical Compounds list.
Next, Empirical Compounds are merged if they have the same m/z, matched form/ion, and retention time. This results in the merged Empirical
Compounds list.
Then, if primary ions are enforced, only Empirical Compounds containing at least 1 primary ion are kept. Primary ions considered are 'M+H[1+]',
'M+Na[1+]', 'M-H2O+H[1+]', 'M-H[-]', 'M-2H[2-]', 'M-H2O-H[-]', 'M+H [1+]', 'M+Na [1+]', 'M-H2O+H [1+]', 'M-H [1-]', 'M-2H [2-]', and 'M-H2O-H [1-]'. This
results in the final Empirical Compounds list.
Finally, pathway libraries are converted from "Compound" space to "Empirical Compound" space. This is done by converting all compounds in each
pathway to all Empirical Compound matches. Then the mummichog/GSEA algorithms work as before to calculate pathway enrichment.
2.1 Start Functional Analysis
Click here to
start
2.2 Starting from a list
3. Click submit to
continue
2.2.2 Data Integrity Check
1. Check Data
Integrity Result to
make sure correct
2. Click Proceed to
continue
2.2.3 Set Parameters
Set corresponding
parameters/Library
1. Switch to uploading
peak intensity table tab
3. Click submit to
continue
2.3.2 Peak uploading – Preprocessing
1. Select Heatmaps
radio to start !
2.3.5 Heatmap based pattern specific analysis - result
This section maybe too complicated to easily understand/follow for beginners, why not watch a video first?
Overview of
the peak
across the
spectrum.
Default is Pattern based
clustered Enrichment
based on p analysis panel.
value.
Focus view of
the specific peak
pattern from the
whole spectrum. Dynamic Display panel,
Default is the top used to show the
Sample
50 peaks. peak/sample information
Names dynamically.
View panel.
2.3.7 Heatmap peak clustering
3. The selected
1. Select a clustering specific spectral
method from Cluster peaks peaks’ pattern
menu, e.g. Ward’s method appears in the
focus view panel.
4. New pattern
appears immediately
at the top focus view
panel. Hold your left
click mouse and move
over from this focus
view to confirm the
area.
6. Click To
Focus view.
3. Left click and move
your mouse over a
new pattern from the 5. Newly selected area will be
overview panel. stitched with the previously 7. The stitched peaks’
selected pattern at the pattern will be presented
bottom focus view. at the Focus view.
2.3.8 Enrichment Analysis -1
TIP1 : Operation Mode could
show the hits in different way.
Annotate will annotate directly
1. Select on the heatmap, while the
database and ‘Extract’ will extract the ions
submit to do hits and hide other non-hits.
the enrichment
analysis
If you have any questions please read through the FAQs or contact us at
Zhiqiang.pang[at]xialab.ca or Jeff.xia[at]xialab.ca