MWASTools: an R/bioconductor package for metabolome-wide association studies

Apparent Saturation of Branched-Chain Amino Acid Catabolism After High Dietary Milk Protein Intake in Healthy Adults

CONTEXT

Milk protein contains high concentrations of branched-chain amino acids (BCAA) that play a critical role in anabolism and are implicated in the onset of obesity and chronic disease. Characterizing BCAA catabolism in the postprandial phase could elucidate the impact of protein intake on obesity risk established in the "early protein hypothesis."

OBJECTIVE

To examine the acute effects of protein content of young child formulas as test meals on BCAA catabolism, observing postprandial plasma concentrations of BCAA in relation to their degradation products.

METHODS

The TOMI Add-On Study is a randomized, double-blind crossover study in which 27 healthy adults consumed 2 isocaloric young child formulas with alternating higher (HP) and lower (LP) protein and fat content as test meals during separate interventions, while 9 blood samples were obtained over 5 hours. BCAA, branched-chain α-keto acids (BCKA), and acylcarnitines were analyzed using a fully targeted HPLC-ESI-MS/MS approach.

RESULTS

Mean concentrations of BCAA, BCKA, and acylcarnitines were significantly higher after HP than LP over the 5 postprandial hours, except for the BCKA α-ketoisovalerate (KIVA). The latter metabolite showed higher postprandial concentrations after LP. With increasing mean concentrations of BCAA, concentrations of corresponding BCKA, acylcarnitines, and urea increased until a breakpoint was reached, after which concentrations of degradation products decreased (for all metabolites except valine and KIVA and Carn C4:0-iso).

CONCLUSION

BCAA catabolism is markedly influenced by protein content of the test meal. We present novel evidence for the apparent saturation of the BCAA degradation pathway in the acute postprandial phase up to 5 hours after consumption.

Bucket Fuser: Statistical Signal Extraction for 1D 1H NMR Metabolomic Data

Untargeted metabolomics is a promising tool for identifying novel disease biomarkers and unraveling underlying pathomechanisms. Nuclear magnetic resonance (NMR) spectroscopy is particularly suited for large-scale untargeted metabolomics studies due to its high reproducibility and cost effectiveness. Here, one-dimensional (1D) 1H NMR experiments offer good sensitivity at reasonable measurement times. Their subsequent data analysis requires sophisticated data preprocessing steps, including the extraction of NMR features corresponding to specific metabolites. We developed a novel 1D NMR feature extraction procedure, called Bucket Fuser (BF), which is based on a regularized regression framework with fused group LASSO terms. The performance of the BF procedure was demonstrated using three independent NMR datasets and was benchmarked against existing state-of-the-art NMR feature extraction methods. BF dynamically constructs NMR metabolite features, the widths of which can be adjusted via a regularization parameter. BF consistently improved metabolite signal extraction, as demonstrated by our correlation analyses with absolutely quantified metabolites. It also yielded a higher proportion of statistically significant metabolite features in our differential metabolite analyses. The BF algorithm is computationally efficient and it can deal with small sample sizes. In summary, the Bucket Fuser algorithm, which is available as a supplementary python code, facilitates the fast and dynamic extraction of 1D NMR signals for the improved detection of metabolic biomarkers.

Density and Diversity Differences of Contemporary and Subfossil Cladocera Assemblages: A Case Study in an Oxbow Lake

Cladocerans are biological indicators of environmental changes. Their remains provide information on past changes in lake environments. We studied the correspondence between contemporary Cladocera assemblages and their subfossil remains from an oxbow lake. We sought to demonstrate that there were differences among the various sites of an oxbow lake with different utilization based on contemporary and subfossil Cladocera assemblages and physical–chemical variables. The oxbow lake’s two sides are used as fishing sites, where angling is the main activity. The middle site of the lake is under nature protection with high macrovegetation coverage. Contemporary and subfossil Cladocera assemblages were sampled from 21 sampling sites along the oxbow lake. Our research showed that the subfossil Cladocera assemblages had higher species richness and densities (36 taxa) than the contemporary species (29 taxa). We found one species of the Polyphemidae family only in the contemporary assemblage. Among the sites, Cladocera assemblages differed in their species composition and density. The highest densities were found in the second fishing site due to the appearance of the small-sized Bosmids. The relationship between Cladocerans and the physical–chemical variables showed that some variables, such as chlorophyll-a, biological oxygen demand, dissolved oxygen, copper, phosphide, and organic matter content, significantly affected Cladocera composition. We found that the subfossil Cladocera assemblage was significantly more diverse and abundant than the contemporary one, indicating that an integrated sampling may be sufficient to provide better results on the total species composition of the water body.

Music of metagenomics-a review of its applications, analysis pipeline, and associated tools

This humble effort highlights the intricate details of metagenomics in a simple, poetic, and rhythmic way. The paper enforces the significance of the research area, provides details about major analytical methods, examines the taxonomy and assembly of genomes, emphasizes some tools, and concludes by celebrating the richness of the ecosystem populated by the "metagenome."

Software Application Profile: exposomeShiny—a toolbox for exposome data analysis

Motivation

Studying the role of the exposome in human health and its impact on different omic layers requires advanced statistical methods. Many of these methods are implemented in different R and Bioconductor packages, but their use may require strong expertise in R, in writing pipelines and in using new R classes which may not be familiar to non-advanced users. ExposomeShiny provides a bridge between researchers and most of the state-of-the-art exposome analysis methodologies, without the need of advanced programming skills.

Implementation

ExposomeShiny is a standalone web application implemented in R. It is available as source files and can be installed in any server or computer avoiding problems with data confidentiality. It is executed in RStudio which opens a browser window with the web application.

General features

The presented implementation allows the conduct of: (i) data pre-processing: normalization and missing imputation (including limit of detection); (ii) descriptive analysis; (iii) exposome principal component analysis (PCA) and hierarchical clustering; (iv) exposome-wide association studies (ExWAS) and variable selection ExWAS; (v) omic data integration by single association and multi-omic analyses; and (vi) post-exposome data analyses to gain biological insight for the exposures, genes or using the Comparative Toxicogenomics Database (CTD) and pathway analysis.

Availability

The exposomeShiny source code is freely available on Github at [https://github.com/isglobal-brge/exposomeShiny], Git tag v1.4. The software is also available as a Docker image [https://hub.docker.com/r/brgelab/exposome-shiny], tag v1.4. A user guide with information about the analysis methodologies as well as information on how to use exposomeShiny is freely hosted at [https://isglobal-brge.github.io/exposome_bookdown/].

A computationally efficient Bayesian seemingly unrelated regressions model for high-dimensional quantitative trait loci discovery

Our work is motivated by the search for metabolite quantitative trait loci (QTL) in a cohort of more than 5000 people. There are 158 metabolites measured by NMR spectroscopy in the 31-year follow-up of the Northern Finland Birth Cohort 1966 (NFBC66). These metabolites, as with many multivariate phenotypes produced by high-throughput biomarker technology, exhibit strong correlation structures. Existing approaches for combining such data with genetic variants for multivariate QTL analysis generally ignore phenotypic correlations or make restrictive assumptions about the associations between phenotypes and genetic loci. We present a computationally efficient Bayesian seemingly unrelated regressions model for high-dimensional data, with cell-sparse variable selection and sparse graphical structure for covariance selection. Cell sparsity allows different phenotype responses to be associated with different genetic predictors and the graphical structure is used to represent the conditional dependencies between phenotype variables. To achieve feasible computation of the large model space, we exploit a factorisation of the covariance matrix. Applying the model to the NFBC66 data with 9000 directly genotyped single nucleotide polymorphisms, we are able to simultaneously estimate genotype-phenotype associations and the residual dependence structure among the metabolites. The R package BayesSUR with full documentation is available at https://cran.r-project.org/web/packages/BayesSUR/.

Chronic Kidney Disease Cohort Studies: A Guide to Metabolome Analyses

Kidney diseases still pose one of the biggest challenges for global health, and their heterogeneity and often high comorbidity load seriously hinders the unraveling of their underlying pathomechanisms and the delivery of optimal patient care. Metabolomics, the quantitative study of small organic compounds, called metabolites, in a biological specimen, is gaining more and more importance in nephrology research. Conducting a metabolomics study in human kidney disease cohorts, however, requires thorough knowledge about the key workflow steps: study planning, sample collection, metabolomics data acquisition and preprocessing, statistical/bioinformatics data analysis, and results interpretation within a biomedical context. This review provides a guide for future metabolomics studies in human kidney disease cohorts. We will offer an overview of important a priori considerations for metabolomics cohort studies, available analytical as well as statistical/bioinformatics data analysis techniques, and subsequent interpretation of metabolic findings. We will further point out potential research questions for metabolomics studies in the context of kidney diseases and summarize the main results and data availability of important studies already conducted in this field.

NMR-Based Metabolomics in Cancer Research

Nuclear magnetic resonance (NMR) spectroscopy offers reproducible quantitative analysis and structural identification of metabolites in various complex biological samples, such as biofluids (plasma, serum, and urine), cells, tissue extracts, and even intact organs. Therefore, NMR-based metabolomics, a mainstream metabolomic platform, has been extensively applied in many research fields, including pharmacology, toxicology, pathophysiology, nutritional intervention, disease diagnosis/prognosis, and microbiology. In particular, NMR-based metabolomics has been successfully used for cancer research to investigate cancer metabolism and identify biomarker and therapeutic targets. This chapter highlights the innovations and challenges of NMR-based metabolomics platform and its applications in cancer research.

NMR spectroscopy in wine authentication: An official control perspective

Wine authentication is vital in identifying malpractice and fraud, and various physical and chemical analytical techniques have been employed for this purpose. Besides wet chemistry, these include chromatography, isotopic ratio mass spectrometry, optical spectroscopy, and nuclear magnetic resonance (NMR) spectroscopy, which have been applied in recent years in combination with chemometric approaches. For many years, ² H NMR spectroscopy was the method of choice and achieved official recognition in the detection of sugar addition to grape products. Recently, ¹ H NMR spectroscopy, a simpler and faster method (in terms of sample preparation), has gathered more and more attention in wine analysis, even if it still lacks official recognition. This technique makes targeted quantitative determination of wine ingredients and nontargeted detection of the metabolomic fingerprint of a wine sample possible. This review summarizes the possibilities and limitations of ¹ H NMR spectroscopy in analytical wine authentication, by reviewing its applications as reported in the literature. Examples of commercial and open-source solutions combining NMR spectroscopy and chemometrics are also examined herein, together with its opportunities of becoming an official method.

pJRES Binning Algorithm (JBA): a new method to facilitate the recovery of metabolic information from pJRES 1H NMR spectra

Motivation

Data processing is a key bottleneck for ¹H NMR-based metabolic profiling of complex biological mixtures, such as biofluids. These spectra typically contain several thousands of signals, corresponding to possibly few hundreds of metabolites. A number of binning-based methods have been proposed to reduce the dimensionality of 1 D ¹H NMR datasets, including statistical recoupling of variables (SRV). Here, we introduce a new binning method, named JBA (“pJRES Binning Algorithm”), which aims to extend the applicability of SRV to pJRES spectra.

Results

The performance of JBA is comprehensively evaluated using 617 plasma ¹H NMR spectra from the FGENTCARD cohort. The results presented here show that JBA exhibits higher sensitivity than SRV to detect peaks from low-abundance metabolites. In addition, JBA allows a more efficient removal of spectral variables corresponding to pure electronic noise, and this has a positive impact on multivariate model building

Availability and implementation

The algorithm is implemented using the MWASTools R/Bioconductor package.

Supplementary information

Supplementary data are available at Bioinformatics online.

Analysis of NMR Metabolomics Data

In this chapter, we summarize data preprocessing and data analysis strategies used for analysis of NMR data for metabolomics studies. Metabolomics consists of the analysis of the low molecular weight compounds in cells, tissues, or biological fluids, and has been used to reveal biomarkers for early disease detection and diagnosis, to monitor interventions, and to provide information on pathway perturbations to inform mechanisms and identifying targets. Metabolic profiling (also termed metabotyping) involves the analysis of hundreds to thousands of molecules using mainly state-of-the-art mass spectrometry (MS) and nuclear magnetic resonance (NMR) spectroscopy technologies. While NMR is less sensitive than mass spectrometry, NMR does provide a wealth of complex and information rich metabolite data. NMR data together with the use of conventional statistics, modeling methods, and bioinformatics tools reveals biomarker and mechanistic information. A typical NMR spectrum, with up to 64k data points, of a complex biological fluid or an extract of cells and tissues consists of thousands of sharp signals that are mainly derived from small molecules. In addition, a number of advanced NMR spectroscopic methods are available for extracting information on high molecular weight compounds such as lipids or lipoproteins. There are numerous data preprocessing, data reduction, and analysis methods developed and evolving in the field of NMR metabolomics. Our goal is to provide an extensive summary of NMR data preprocessing and analysis strategies by providing examples and open source and commercially available analysis software and bioinformatics tools.

The metaRbolomics Toolbox in Bioconductor and beyond

Metabolomics aims to measure and characterise the complex composition of metabolites in a biological system. Metabolomics studies involve sophisticated analytical techniques such as mass spectrometry and nuclear magnetic resonance spectroscopy, and generate large amounts of high-dimensional and complex experimental data. Open source processing and analysis tools are of major interest in light of innovative, open and reproducible science. The scientific community has developed a wide range of open source software, providing freely available advanced processing and analysis approaches. The programming and statistics environment R has emerged as one of the most popular environments to process and analyse Metabolomics datasets. A major benefit of such an environment is the possibility of connecting different tools into more complex workflows. Combining reusable data processing R scripts with the experimental data thus allows for open, reproducible research. This review provides an extensive overview of existing packages in R for different steps in a typical computational metabolomics workflow, including data processing, biostatistics, metabolite annotation and identification, and biochemical network and pathway analysis. Multifunctional workflows, possible user interfaces and integration into workflow management systems are also reviewed. In total, this review summarises more than two hundred metabolomics specific packages primarily available on CRAN, Bioconductor and GitHub.

Metabolomics-Based Dietary Biomarkers in Nutritional Epidemiology-Current Status and Future Opportunities

The application of metabolomics in nutrition epidemiology holds great promise and there is a high expectation that it will play a leading role in deciphering the interactions between diet and health. However, while significant progress has been made in the identification of putative biomarkers, more work is needed to address the use of the biomarkers in dietary assessment. The aim of this review is to critically evaluate progress in these areas and to identify challenges that need to be addressed going forward. The notable applications of dietary biomarkers in nutritional epidemiology include 1) determination of food intake based on biomarkers levels and calibration equations from feeding studies, 2) classification of individuals into dietary patterns based on the urinary metabolic profile, and 3) application of metabolome wide-association studies. Further work is needed to address some specific challenges to enable biomarkers to reach their full potential.

Tools and resources for metabolomics research community: A 2017-2018 update

The scale at which MS- and NMR-based platforms generate metabolomics datasets for both research, core, and clinical facilities to address challenges in the various sciences-ranging from biomedical to agricultural-is underappreciated. Thus, metabolomics efforts spanning microbe, environment, plant, animal, and human systems have led to continual and concomitant growth of in silico resources for analysis and interpretation of these datasets. These software tools, resources, and databases drive the field forward to help keep pace with the amount of data being generated and the sophisticated and diverse analytical platforms that are being used to generate these metabolomics datasets. To address challenges in data preprocessing, metabolite annotation, statistical interrogation, visualization, interpretation, and integration, the metabolomics and informatics research community comes up with hundreds of tools every year. The purpose of the present review is to provide a brief and useful summary of more than 95 metabolomics tools, software, and databases that were either developed or significantly improved during 2017-2018. We hope to see this review help readers, developers, and researchers to obtain informed access to these thorough lists of resources for further improvisation, implementation, and application in due course of time.