Addressing Missing Data in GC × GC Metabolomics: Identifying Missingness Type and Evaluating the Impact of Imputation Methods on Experimental Replication

imputomics: web server and R package for missing values imputation in metabolomics data

MOTIVATION

Missing values are commonly observed in metabolomics data from mass spectrometry. Imputing them is crucial because it assures data completeness, increases the statistical power of analyses, prevents inaccurate results, and improves the quality of exploratory analysis, statistical modeling, and machine learning. Numerous Missing Value Imputation Algorithms (MVIAs) employ heuristics or statistical models to replace missing information with estimates. In the context of metabolomics data, we identified 52 MVIAs implemented across 70 R functions. Nevertheless, the usage of those 52 established methods poses challenges due to package dependency issues, lack of documentation, and their instability.

RESULTS

Our R package, 'imputomics', provides a convenient wrapper around 41 (plus random imputation as a baseline model) out of 52 MVIAs in the form of a command-line tool and a web application. In addition, we propose a novel functionality for selecting MVIAs recommended for metabolomics data with the best performance or execution time.

AVAILABILITY AND IMPLEMENTATION

'imputomics' is freely available as an R package (github.com/BioGenies/imputomics) and a Shiny web application (biogenies.info/imputomics-ws). The documentation is available at biogenies.info/imputomics.

The optimization and comparison of two high-throughput faecal headspace sampling platforms: the microchamber/thermal extractor and hi-capacity sorptive extraction probes (HiSorb)

Disease detection and monitoring using volatile organic compounds (VOCs) is becoming increasingly popular. For a variety of (gastrointestinal) diseases the microbiome should be considered. As its output is to large extent volatile, faecal volatilomics carries great potential. One technical limitation is that current faecal headspace analysis requires specialized instrumentation which is costly and typically does not work in harmony with thermal desorption units often utilized in e.g. exhaled breath studies. This lack of harmonization hinders uptake of such analyses by the Volatilomics community. Therefore, this study optimized and compared two recently harmonized faecal headspace sampling platforms:High-capacity Sorptive extraction (HiSorb) probesand theMicrochamber thermal extractor (Microchamber). Statistical design of experiment was applied to find optimal sampling conditions by maximizing reproducibility, the number of VOCs detected, and between subject variation. To foster general applicability those factors were defined using semi-targeted as well as untargeted metabolic profiles. HiSorb probes were found to result in a faster sampling procedure, higher number of detected VOCs, and higher stability. The headspace collection using the Microchamber resulted in a lower number of detected VOCs, longer sampling times and decreased stability despite a smaller number of interfering VOCs and no background signals. Based on the observed profiles, recommendations are provided on pre-processing and study design when using either one of both platforms. Both can be used to perform faecal headspace collection, but altogether HiSorb is recommended.

Statistical methods and resources for biomarker discovery using metabolomics

Metabolomics is a dynamic tool for elucidating biochemical changes in human health and disease. Metabolic profiles provide a close insight into physiological states and are highly volatile to genetic and environmental perturbations. Variation in metabolic profiles can inform mechanisms of pathology, providing potential biomarkers for diagnosis and assessment of the risk of contracting a disease. With the advancement of high-throughput technologies, large-scale metabolomics data sources have become abundant. As such, careful statistical analysis of intricate metabolomics data is essential for deriving relevant and robust results that can be deployed in real-life clinical settings. Multiple tools have been developed for both data analysis and interpretations. In this review, we survey statistical approaches and corresponding statistical tools that are available for discovery of biomarkers using metabolomics.

The influence of the co-exposure to polycyclic aromatic hydrocarbons and toxic metals on DNA damage in brazilian lactating women and their infants: A cross-sectional study using machine learning approaches

Polycyclic aromatic hydrocarbons (PAHs) and toxic metals are widely spread pollutants of public health concern. The co-contamination of these chemicals in the environment is frequent, but relatively little is known about their combined toxicities. In this context, this study aimed to evaluate the influence of the co-exposure to PAHs and toxic metals on DNA damage in Brazilian lactating women and their infants using machine learning approaches. Data were collected from an observational, cross-sectional study with 96 lactating women and 96 infants living in two cities. The exposure to these pollutants was estimated by determining urinary levels of seven mono-hydroxylated PAH metabolites and the free form of three toxic metals. 8-Hydroxydeoxyguanosine (8-OHdG) levels in the urine were used as the oxidative stress biomarker and set as the outcome. Individual sociodemographic factors were also collected using questionnaires. Sixteen machine learning algorithms were trained using 10-fold cross-validation to investigate the associations of urinary OH-PAHs and metals with 8-OHdG levels. This approach was also compared with models attained by multiple linear regression. The results showed that the urinary concentration of OH-PAHs was highly correlated between the mothers and their infants. Multiple linear regression did not show a statistically significant association between the contaminants and urinary 8OHdG levels. Machine learning models indicated that all investigated variables did not present predictive performance on 8-OHdG concentrations. In conclusion, PAHs and toxic metals were not associated with 8-OHdG levels in Brazilian lactating women and their infants. These novelty and originality results were achieved even after applying sophisticated statistical models to capture non-linear relationships. However, these findings should be interpreted cautiously because the exposure to the studied contaminants was considerably low, which may not reflect other populations at risk.

Comprehensive Two-Dimensional Gas Chromatography as a Bioanalytical Platform for Drug Discovery and Analysis

Over the last decades, comprehensive two-dimensional gas chromatography (GC×GC) has emerged as a significant separation tool for high-resolution analysis of disease-associated metabolites and pharmaceutically relevant molecules. This review highlights recent advances of GC×GC with different detection modalities for drug discovery and analysis, which ideally improve the screening and identification of disease biomarkers, as well as monitoring of therapeutic responses to treatment in complex biological matrixes. Selected recent GC×GC applications that focus on such biomarkers and metabolite profiling of the effects of drug administration are covered. In particular, the technical overview of recent GC×GC implementation with hyphenation to the key mass spectrometry (MS) technologies that provide the benefit of enhanced separation dimension analysis with MS domain differentiation is discussed. We conclude by highlighting the challenges in GC×GC for drug discovery and development with perspectives on future trends.

Big Data in Gastroenterology Research

Studying individual data types in isolation provides only limited and incomplete answers to complex biological questions and particularly falls short in revealing sufficient mechanistic and kinetic details. In contrast, multi-omics approaches to studying health and disease permit the generation and integration of multiple data types on a much larger scale, offering a comprehensive picture of biological and disease processes. Gastroenterology and hepatobiliary research are particularly well-suited to such analyses, given the unique position of the luminal gastrointestinal (GI) tract at the nexus between the gut (mucosa and luminal contents), brain, immune and endocrine systems, and GI microbiome. The generation of 'big data' from multi-omic, multi-site studies can enhance investigations into the connections between these organ systems and organisms and more broadly and accurately appraise the effects of dietary, pharmacological, and other therapeutic interventions. In this review, we describe a variety of useful omics approaches and how they can be integrated to provide a holistic depiction of the human and microbial genetic and proteomic changes underlying physiological and pathophysiological phenomena. We highlight the potential pitfalls and alternatives to help avoid the common errors in study design, execution, and analysis. We focus on the application, integration, and analysis of big data in gastroenterology and hepatobiliary research.

Volatile Metabolites in Lavage Fluid Are Correlated with Cytokine Production in a Valley Fever Murine Model

Coccidioides immitis and Coccidioides posadasii are soil-dwelling fungi of arid regions in North and South America that are responsible for Valley fever (coccidioidomycosis). Forty percent of patients with Valley fever exhibit symptoms ranging from mild, self-limiting respiratory infections to severe, life-threatening pneumonia that requires treatment. Misdiagnosis as bacterial pneumonia commonly occurs in symptomatic Valley fever cases, resulting in inappropriate treatment with antibiotics, increased medical costs, and delay in diagnosis. In this proof-of-concept study, we explored the feasibility of developing breath-based diagnostics for Valley fever using a murine lung infection model. To investigate potential volatile biomarkers of Valley fever that arise from host−pathogen interactions, we infected C57BL/6J mice with C. immitis RS (n = 6), C. posadasii Silveira (n = 6), or phosphate-buffered saline (n = 4) via intranasal inoculation. We measured fungal dissemination and collected bronchoalveolar lavage fluid (BALF) for cytokine profiling and for untargeted volatile metabolomics via solid-phase microextraction (SPME) and two-dimensional gas chromatography coupled with time-of-flight mass spectrometry (GC×GC-TOFMS). We identified 36 volatile organic compounds (VOCs) that were significantly correlated (p < 0.05) with cytokine abundance. These 36 VOCs clustered mice by their cytokine production and were also able to separate mice with moderate-to-high cytokine production by infection strain. The data presented here show that Coccidioides and/or the host produce volatile metabolites that may yield biomarkers for a Valley fever breath test that can detect coccidioidal infection and provide clinically relevant information on primary pulmonary disease severity.