-
Notifications
You must be signed in to change notification settings - Fork 7
/
README.Rmd
107 lines (64 loc) · 6.24 KB
/
README.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
---
title: "Targeted Metabolomics Processing Pipeline"
output:
github_document
---
> Microbial Metabolomics Research Center
[](https://github.com/nhn/tui.editor/issues?q=is%3Aissue+is%3Aopen+label%3A%22help+wanted%22)
## Table of Contents
- [Description](#-Description)
- [Packages and External Processes](#-Packages-and-External-Processes)
- [Usage](#-Usage)
- [User-Defined QC Paramters](#User-Defined-QC-Parameters)
- [Visualization](#-Visualization)
------------------------------------------------------------------------
### Description {#description}
**Important:** *This repository is a template and training tool! Each dataset is different; while this pipeline handles many of the scenarios that arise within the data wrangling process, the user will need to adjust the parameters for their specific input data. This also applies to the example data from the Ingalls Drive. The example csvs are a handful of possible combinations of run type/instrument/software, but there is not an example sheet for every case.*
The Targeted Pipeline is a combined process for importing, tidying, Quality Controlling, BMISing (see below for more details on this step), and quantifying (where possible) targeted metabolomics data produced by the Ingalls Lab.
It is run by a **"Control.Rmd"** markdown script that moves through the pipeline step by step, accounting for two mass spectrometer instrument types ([TQS](https://www.waters.com/waters/en_US/Xevo-TQ-S/nav.htm?cid=10160596&locale=en_US) or [QE](https://www.thermofisher.com/order/catalog/product/IQLAAEGAAPFALGMBDK)), run type (HILIC/Reverse phase), ion mode (positive/negative), and metabolomics processing software ([MS-DIAL](http://prime.psc.riken.jp/compms/msdial/main.html) or [Skyline](https://skyline.ms/wiki/home/software/Skyline/page.view?name=tutorial_hi_res_metabolomics)).
The Control.Rmd script contains detailed instructions, and *the user must read carefully through each section description while progressing through the pipeline*.
The pipeline contains four major sections:
1. Import and cleaning/rearranging of data.
2. Quality control using user-defined parameters.
3. Applying a Best-Matched Internal Standard (B-MIS).
4. Quantifying peak area to umol/vial when possible.
### 📦 Packages and External Processes {#packages-and-external-processes}
Below is a list of packages and external sofwares that this pipeline utilizes. Please ensure you have the package(s) installed and have access to the tools listed below.
```{r packages, echo=FALSE, message=FALSE, warnings=FALSE, results='asis'}
packages.table <- "
|Name | Description
|:-------------|:-------------|
|[`tidyverse`](https://www.tidyverse.org/) | Data science package in R. |
|[`data.table`](https://cran.r-project.org/web/packages/data.table/vignettes/datatable-intro.html) | Provides an enhanced version of data.frame. |
|[`anytime`](https://cran.r-project.org/web/packages/anytime/index.html) | General purpose date converter. This package helps the user track analysis reruns and produce unique identifiers for files. |
|[`rlist`](https://cran.r-project.org/web/packages/rlist/index.html) | Non-tabular data manipulation toolbox. |
|[`BMIS`](https://github.com/IngallsLabUW/B-MIS-normalization) | R script for batch-specific normalization of raw peak areas. Matches measured metabolites with isotope-labeled internal standards that behave similarly during the analysis, nicknamed Best-Matched Internal Standard or B-MIS. This pipeline contains a modified script taken from the original repository. Published paper here: <https://pubs.acs.org/doi/10.1021/acs.analchem.7b04400> |
|[`Example Data`](https://drive.google.com/drive/folders/1E-B4dyDTkwOfupHJK0uL8BgKj586OU31) | Pretidied data from the Ingalls Lab for learning purposes. If you do not have access to this drive, please contact the owners of this repository.|
"
cat(packages.table)
```
### Usage {#usage}
The first few chunks of the Control.Rmd script will create empty folders for different types of data (raw, processed, extra, intermediate). Some of this data will be produced by the script (processed, intermediate) and some needs to be moved into those folders before proceeding. Because this is a template script, the data has been modified for easy input, which will NEVER be the case in a real run.
All the necessary example data can be downloaded here: <https://drive.google.com/drive/folders/1E-B4dyDTkwOfupHJK0uL8BgKj586OU31>
If you do not have access to this drive, please contact the owners of this repository.
> #### User-Defined QC Parameters {#user-defined-QC-paramters}
These parameters will change depending on your data and level of stringency. Below are the "default" values.
**area.min**: The minimum area of an integrated peak to be considered a real peak.
HILIC - 1000, Reverse Phase - 5000.
**RT.flex**: Flexibility of the deviation of the sample retention time from the appropriate standard. HILIC +/- 0.4 min, Reverse Phase +/- 0.2 min.
**blk.thresh**: Cutoff ratio value of sample area to a blank reference sample. HILIC +/-0.3, Reverse Phase +/- 0.2.
**SN.min**: Cutoff ratio value for Signal to Noise. HILIC +/- 4, Reverse Phase +/- 4.
#### Additional QC parameters for Skyline
**height.min**: Minimum height for a peak. 1000, HILIC and Reverse Phase.
**height.max**: Maximum height for a peak. 1.0e8, HILIC and Reverse Phase.
#### Additional QC parameter for TQS
**IR.flex**: The TQS produces quantitative and secondary trace values. Find the Ion Ratio by dividing the area of the quantitative trace by the area of the secondary trace for standards, and use these ratios to create a reference table. Deviations from this reference table that are larger than the user-defined values will be flagged. 0.3, HILIC and Reverse Phase.
## Visualization {#visualization}
------------------------------------------------------------------------
## 🔧 Pull Requests
Pull requests are welcome. For major changes, please open an issue first.
## 💬 Contributors
- [Regina Lionheart](https://github.com/R-Lionheart)
- [Katherine Heal](https://github.com/kheal)
- [Angie Boysen](https://github.com/Angie-B)