mrbin - Metabolomics Data Analysis Functions

This package is a collection of functions for processing and analyzing metabolomics data.

The namesake function mrbin() uses spectral binning to convert 1D or 2D Nuclear Magnetic Resonance (NMR) data into a matrix of values suitable for further data analysis and performs basic processing steps in a reproducible way. Negative values, a common issue in NMR data, are replaced by positive values. All used parameters are stored in a readable text file and can be restored from that file to enable exact reproduction of the data at a later time.

The atnv algorithm for replacing negative values in NMR data sets can be employed using atnv().

NMR plotting functions are found in mrplot().

Artificial Neural Network features can be analyzed using Feature Impact Assessment (FIA) using the function fia().

Installation

To install mrbin, please install the latest version of R first. Then install the package as follows:

To install the latest stable version from CRAN:

install.packages("mrbin")
library("mrbin")

You can find more details and examples in the vignette file:

vignette("mrbin")

mrbin: Magnetic Resonance Binning, Integration and Normalization

To use this package, you will need your NMR data in the Bruker file format accessible on your computer. Please make sure your data is Fourier transformed, phase corrected, baseline corrected, and correctly referenced. The data has to be stored in folders according to standard Bruker folders, that means foldername/1/pdata/1 etc. Experiment numbers and processing numbers can be freely chosen.

This package has been tested for 1D NOESY and 2D 1H-13C HSQC spectra.

Before starting mrbin, take a look at your NMR data, for example in Bruker Topspin, and decide on the following:

• Bin area: Area where signals are observed in your data set
• Bin width: Should match roughly the width of a singlet peak in your data set. Given in ppm.
• Bin height (only 2D): Should match roughly the height of a singlet peak in your data set. Given in ppm.
• Solvent area: Area to exclude to remove solvent artifacts
• Additional areas to be removed: Any other area containing artifacts, such as streaks surrounding strong peaks.

mrbin will also show you preview plots for these parameters during the run.

You can start mrbin using the following code:

results<-mrbin()

This will start a series of questions that will guide you through the parameters to be used, including preview plots.

mrbin() returns an (invisible) mrbin object containing the following variables:

• bins: A matrix containing bin data for all samples, Depending on the option you chose, the data will be cleaned up and scaled.
• parameters: A list containing all parameters used to create the bin matrix.
• metadata: A list containing metadata, if provided.
• transformations: A character vector containing information on the data transformation and scaling that has been performed, for example reference scaling, PQN, atnv, log transform, etc.
• changeLog: A data.frame containing information on documented changes that were made to the data, including time stamps.
• changeValues: A list containing control values, enabling verifying changes by checkmrbin(mrbinResults)

Several files may be written to the chosen directory:

• A .txt file containing all parameters and potential warning messages from the mrbin run. This file can be reloaded to R using recreatemrbin("filename"). This will enable reusing parameters used in a previous run and can help increase reproducibility.
• A .Rdata file containing the generated mrbin data object.
• A .R text file containing R code to recreate and or load the data set.
• A .pdf file containing quality control plots of the raw binned data
• A .pdf file containing quality control plots of the processed binned data
• A .pdf file containing preview plots of the chosen signal-to-noise levels of selected spectra

Submitting Parameters at the Command Line

Parameters can be submitted at the command line, using the following syntax:

results<-mrbin(parameters=list(dimension="1D",
             binwidth1D=.01,
             referenceScaling="Yes",
             removeSolvent="Yes",
             removeAreas="No",
             sumBins="No",
             noiseRemoval="Yes",
             PQNScaling="Yes",
             fixNegatives="Yes",
             logTrafo="Yes",
             NMRfolders=c("C:/NMR/Sample01/1/pdata/1",
                          "C:/NMR/Sample19/1/pdata/1",
                          "C:/NMR/Sample61/1/pdata/1")
     ))

Users will then be asked to adjust the parameters as above. Additional parameters are available, please see vignette("mrbin") for details.

atnv: Affine Transformation of Negative Values

Please find details on the atnv algorithm in vignette("mrbin").

fia: Feature Impact Assessment

Please find details on the fia algorithm in vignette("mrbin").

Known Issues

Firewall Warnings

The function mrbin uses the socket approach for parallel computing. This requires establishing network connections to the local cluster, which might trigger the firewall. It is safe to unblock these connections.

Pop-Up Windows

mrbin is set up to ask for user input through pop-up windows. This requires graphics support, otherwise the user input will be asked through command line menus, which is less user friendly but still offers the full functionality.

Apple/Mac Computers And RStudio

In some cases, running mrbin from within RStudio on Apple computers will not generate pop-up windows. To enable pop-up windows, it might be helpful to install the newest version of xquartz from https://www.xquartz.org.

Built With

• [roxygen2]
• [devtools]

Authors

• Matthias S. Klein - KleinOmicsLab

Citation

If you are using mrbin in a publication, please cite the following manuscript: Klein, M.S. (2021): Affine Transformation of Negative Values for NMR Metabolomics Using the mrbin R Package. J. Proteome Res. 20(2):1397-1404, DOI: 10.1021/acs.jproteome.0c00684

License

This project is licensed under GPL-3.0.