Q: What is GlyTrait?
A: GlyTrait is a tool to calculate derived traits for N-glycomic data.
Q: Sounds cool! So..., what are derived traits again?
A: Well, derived traits are artificial variables that summarize certain aspects of the glycome. For example, the proportion of core-fucosylated glycans, the average number of sialic acids per glycan, or the proportion of bisected glycans within bi-antennary complex glycans, etc. Derived traits are more biologically relavant and have been used a lot in the glycomics community.
Q: So, GlyTrait does the dirty work for me, saving my time and energy for more interesting analysis?
A: You bet!
Let there be no code!
-- my colleagues
Click the badge above and try the GlyTrait web app! The web app is self-documented, but in case you need help, please refer to the Usage section below.
python >= 3.10
If python hasn't been installed, download it from its website, or use Anaconda if you like.
pipx is a tool to help you install and run end-user applications written in Python.
It's roughly similar to macOS's brew, JavaScript's npx, and Linux's apt.
Install pipx following its Document.
pipx install glytraitconda create -n glytrait python=3.10and activate the environment.
conda activate glytraitpip install glytrait-
Download the example files to a new directory:
glytrait_example. -
Open a terminal and navigate to the directory (you need to replace
path/to/glytrait_examplewith the actual path):
cd path/to/glytrait_example- Run the following command:
glytrait abundance.csv --glycan-file=structures.csvThat's it! If everything goes well, a folder named "abundance_glytrait" will be created in the same directory with the abundnce.csv file. Inside the directory are four files:
derived_traits.csv: all derived traits calculated by GlyTrait.derived_traits_filtered.csv: derived traits after post-filtering (see Post-filtering).glycan_abundance_processed.csv: the glycan abundance after preprocessing (see Preprocessing).meta_properties.csv: the meta-properties of all glycans.
The detailed format of the input file will be introduced in the Input file format section.
This section is intended to give an overview of the CLI interface. Feel free to skip it right now.
As a glance, GlyTrait supports the following options:
| Option | Description |
|---|---|
| --help | Show the help message and exit. |
| --glycan-file | The glycan structure or composition file. |
| --mp-file | The meta-property file. |
| -m, --mode | The mode. "S" or "structure" for structure mode, "C" or "composition" for composition mode. Default: structure. |
| -o, --output | The output path. Default: the same directory with the input file. |
| -r, --filter-glycan-ratio | The proportion of missing values for a glycan to be ruled out. Default: 0.5. |
| -i, --inpute-method | The imputation method. "min", "mean", "median", "zero", or "lod". Default: "min". |
| -c, --corr-threshold | The correlation threshold for collinearity filtering. Default to 1.0. |
| -l, --sia-linkage | Flag to include the sialic acid linkage traits. |
| --no-filter | Flag to turn off post-filtering of derived traits. |
| -g, --group | The group file. |
| -f, --formula-file | The custom formula file to use. |
| -b, --builtin-formulas | The directory path to save a copy of the built-in formulas. |
The following sections will introduce these options in detail.
GlyTrait has two modes: the "structure" mode and the "composition" mode. In the "structure" mode, GlyTrait will calculate derived traits based on the topology properties of glycan structures. In the "composition" mode, GlyTrait will make educated guesses on the structure properties based on the glycan composition.
Note that the "composition" mode has uncertainties to some extent. Specifically:
- Estimating the number of Gal based on composition is not possible for hybrid glycans, so GlyTrait will calculate the number of Gal assuming there are no hybrid glycans. kily, hybrid glycans are usually in low abundance, so the algorithm is a good approximation for most cases.
- Estimating the number of branches is not possible based on composition, so GlyTrait will roughly classify glycans into two categories: low-branching and high-branching. Glycans with N > 4 (including bisecting diantenary glycans) are considered as high-branching, while others as low-branching.
- Telling hybrid glycans from mono-antenary complex glycans is not possible based on composition, so GlyTrait will not classify glycans into complex, hybrid and high-mannose.
Due to the ambiguities above, we recommend using the "structure" mode if possible.
You can specify the mode by the "-m" or the "--mode" option:
glytrait abundance.csv --glycan-file=composition.csv -m compositionOr in short:
glytrait abundance.csv --glycan-file=composition.csv -m CThe default mode is the "structure" mode, as in the quick start example.
Thus, using glytrait -m structure or glytrait -m S is equivalent to glytrait alone.
If you might use both modes in a project,
we recommend using the "-m" option to specify the mode explicitly.
At least two files are needed for GlyTrait to work:
A csv file with samples as rows and glycan IDs as columns. An example file would be like:
| Sample | Glycan1 | Glycan2 | Glycan3 |
|---|---|---|---|
| Sample1 | 0.0417 | 0.0503 | 0.0354 |
| Sample2 | 0.0233 | 0.0533 | 0.0593 |
| Sample3 | 0.0123 | 0.0133 | 0.0194 |
The header of the first column should be "Sample", and the header of the other columns should be glycan IDs. Glycan IDs can be any string, e.g., the composition strings ("H3N4").
Both glycan IDs and samples should be unique.
A csv file with two columns: "GlycanID" and "Structure" (or "Composition"). An example file would be like:
| GlycanID | Structure |
|---|---|
| Glycan1 | RES... |
| Glycan2 | RES... |
| Glycan3 | RES... |
The "GlycanID" column should contain all glycan IDs in the abundance file. The "Structure" column should contain the structure strings of the glycans. For now, only the GlycoCT format is supported. In the "composition" mode, the second column should be "Composition" instead of "Structure", and the composition strings should be used instead of the structure strings. Condensed format ("H3N4F1S1") is supported.
The names of the abundance file and the structure (composition) file is not important, i.e., it doesn't have to be "abundance.csv" and "structures.csv". You could use any name you want, as long as the two files are passed in order, e.g.:
glytrait experiment_1_15.csv --glycan-file=serum_structures.csv
# ~~~~~~~~~~~~~~~~~~~ ~~~~~~~~~~~~~~~~~~~~
# the abundance file the structure fileYou might have noticed before that GlyTrait saves the output file to the same directory as the abundance file with a "_glytrait" suffix. You can specify the output file path by using the "-o" or "--output-file" option:
glytrait abundance.csv --glycan-file=structure.csv -o outputGlyTrait will carry out a preprocessing step before calculating derived traits. The following steps will be done:
- Remove glycans with missing values in more than a certain proportion of samples.
- Impute missing values.
- Perform Total Abundance Normalization.
In the glycan-filtering step, the proportion threshold could be specified by the "-r" or the "--filter-glycan-ratio" option. The default value is 1, which means no glycan will be removed. You can change this value to 0.5 by:
glytrait abundance.csv --glycan-file=structure.csv -r 0.5The imputation method could be specified by the "-i" or the "--impute-method" option. The default method is "zero", which means missing values will be imputed by 0. Other supported methods are "mean", "median", "zero", "lod". You can change the imputation method to "min" by:
glytrait abundance.csv --glycan-file=structure.csv -i minA full list of supported imputation methods are:
- "min": impute missing values by the minimum value of a glycan within all samples.
- "mean": impute missing values by the mean value of a glycan within all samples.
- "median": impute missing values by the median value of a glycan within all samples.
- "zero": impute missing values by 0.
- "lod": impute missing values by the limit of detection (LOD) of the equipment. The LOD of a glycan is defined as the minimum value of the glycan within all samples divided by 5.
Sialic acids can have different linkages for N-glycans (e.g., α2,3 and α2,6). Different sialic acid linkage has different biological functions. GlyTrait supports calculating derived traits regarding these linkages. To use this feature, you need to have siaic acid linkage information.
In the structure mode, the "Structure" column or the structure file should contain the linkage information. Only linkage information about sialic acids is needed. This can be easily done using GlycoWorkbench.
In the composition mode, the "Composition" column must contain the linkage information. GlyTrait uses a common notation for sialic acid with different linkages: "E" for a2,6-linked sialic acids, and "L" a2,3-linked sialic acids. For example, "H5N4F1E1L1" contains two sialic acids, one is a2,6-linked, and the other is a2,3-linked.
You can use the "-l" or "--sia-linkage" option to include sialic-acid-linkage traits:
glytrait abundance.csv --glycan-file=structure.csv -lNote that if you use this option, all glycans with sialic acids should have linkage information. That is to say, all structure strings should have structure information in the structure mode and no "S" in composition strings in the composition mode.
Not all derived traits are informative. For example, some traits might have the same value for all samples. Some traits might be highly correlated with others.
GlyTrait carries out a two-step post-filtering process to remove these uninformative traits. First, traits with the same value for all samples will be removed. Second, highly correlated traits will be pruned, keeping only the traits considering more glycans.
GlyTrait filters out highly correlated traits, using a "trait family tree" filtering method. Briefly, for two correlated traits, the "parent" trait, which normally considers more glycans, will be kept. For example, for the two high-correlated traits: A2FG and A2G, the latter will be kept because it is more general, and more robust for considering more glycans. Thanks to the dynamic "trait family tree" generated by GlyTrait, user-defined traits will also be considered in this filtering process.
By default, GlyTrait only filters trarits with Pearson correlation coefficient of 1, i.e., traits with perfect collinearity. This threshold can be changed by the "-c" or "--corr-threshold" option:
glytrait abundance.csv --glycan-file=structure.csv -c 0.9Setting the threshold to -1 will turn off the colinearity filtering. To turn off postfiltering all together, use the "--no-filtering" option:
glytrait abundance.csv --glycan-file=structure.csv --no-filteringComing soon...
(The functionallity has been implemented, but the documentation is not ready yet. After publishing the GlyTrait paper, we will update the documentation.)
Coming soon...
(The functionallity has been implemented, but the documentation is not ready yet. After publishing the GlyTrait paper, we will update the documentation.)
To reproduce the results in our paper "GlyTrait: A versatile bioinformatics tool for Glycomics Analysis", please download the Supplementary Data files. There are three such files, corresponding to different datasets mentioned in the paper. Each Supplementary Data file has several sheets with the prefix of "Input -". These are the input data for GlyTrait. Please reformat the data according to the Input file format section, and run GlyTrait with the same options as in the paper.