Deciphering the Underlying Metabolomic and Lipidomic Patterns Linked to Thermal Acclimation in Saccharomyces cerevisiae

Multivariate curve resolution-based data fusion approaches applied in 1H NMR metabolomic analysis of healthy cohorts

BACKGROUND

Metabolomics plays a critical role in deciphering metabolic alterations within individuals, demanding the use of sophisticated analytical methodologies to navigate its intricate complexity. While many studies focus on single biofluid types, simultaneous analysis of multiple matrices enhances understanding of complex biological mechanisms. Consequently, the development of data fusion methods enabling multiblock analysis becomes essential for comprehensive insights into metabolic dynamics.

RESULTS

This study introduces a novel guideline for jointly analyzing diverse metabolomic datasets (serum, urine, metadata) with a focus on metabolic differences between groups within a healthy cohort. The guideline presents two fusion strategies, 'Low-Level data fusion' (LLDF) and 'Mid-Level data fusion' (MLDF), employing a sequential application of Multivariate Curve Resolution with Alternating Least Squares (MCR-ALS), linking the outcomes of successive analyses. MCR-ALS is a versatile method for analyzing mixed data, adaptable at various stages of data processing-encompassing resonance integration, data compression, and exploratory analysis. The LLDF and MLDF strategies were applied to 1H NMR spectral data extracted from urine and serum samples, coupled with biochemical metadata sourced from 145 healthy volunteers.

SIGNIFICANCE

Both methodologies effectively integrated and analysed multiblock datasets, unveiling the inherent data structure and variables associated with discernible factors among healthy cohorts. While both approaches successfully detected sex-related differences, the MLDF strategy uniquely revealed components linked to age. By applying this analysis, we aim to enhance the interpretation of intricate biological mechanisms and uncover variations that may not be easily discernible through individual data analysis.

Integrated metabolomics and lipidomics analyses suggest the temperature-dependent lipid desaturation promotes aflatoxin biosynthesis in Aspergillus flavus

Temperature is one of the main factors affecting aflatoxin (AF) biosynthesis in Aspergillus flavus. Previous studies showed that AF biosynthesis is elevated in A. flavus at temperatures between 28°C-30°C, while it is inhibited at temperatures above 30°C. However, little is known about the metabolic mechanism underlying temperature-regulated AF biosynthesis. In this study, we integrated metabolomic and lipidomic analyses to investigate the endogenous metabolism of A. flavus across 6 days of mycelia growth at 28°C (optimal AF production) and 37°C (no AF production). Results showed that both metabolite and lipid profiles were significantly altered at different temperatures. In particular, metabolites involved in carbohydrate and amino acid metabolism were up-regulated at 37°C on the second day but down-regulated from days three to six. Moreover, lipidomics and targeted fatty acids analyses of mycelia samples revealed a distinct pattern of lipid species and free fatty acids desaturation. High degrees of polyunsaturation of most lipid species at 28°C were positively correlated with AF production. These results provide new insights into the underlying metabolic changes in A. flavus under temperature stress.

Quantification strategies for two-dimensional liquid chromatography datasets using regions of interest and multivariate curve resolution approaches

In this work, three chemometrics-based approaches are compared for quantification purposes when using two-dimensional liquid chromatography (LC×LC-MS), taking as a study case the quantification of amino acids in commercial drug mixtures. Although the approaches have been already used for one-dimensional gas or liquid chromatography, the main novelty of this work is the demonstration of their applicability to LC×LC-MS datasets. Besides, steps such as peak alignment and modelling, commonly applied in this type of data analysis, are not required with the approaches proposed here. In a first step, regions of interest (ROI) strategy is used for the spectral compression of the LC×LC-MS datasets. Then the first strategy consists of building a calibration curve from the areas obtained in this ROI compression step. Alternatively, the ROI intensity matrices can be used as input for a second analysis step employing the multivariate curve resolution alternating least squares (MCR-ALS) method. The main benefit of MCR-ALS is the resolution of elution and spectral profiles for each of the analytes in the mixture, even in the case of strong coelutions and high signal overlapping. Classical MCR-ALS based calibration curve from the peak areas resolved only applying non-negativity constraints (second strategy) is compared to the results obtained when an area correlation constraint is imposed during the ALS optimization (third strategy). All in all, similar quantification results were achieved by the three approaches but, especially in prediction studies, the more accurate quantification is obtained when the calibration curve is built from the peak areas obtained with MCR-ALS when the area correlation constraint is imposed.

Comparison of Multivariate ANOVA-Based Approaches for the Determination of Relevant Variables in Experimentally Designed Metabolomic Studies

The use of chemometric methods based on the analysis of variances (ANOVA) allows evaluation of the statistical significance of the experimental factors used in a study. However, classical multivariate ANOVA (MANOVA) has a number of requirements that make it impractical for dealing with metabolomics data. For this reason, in recent years, different options have appeared that overcome these limitations. In this work, we evaluate the performance of three of these multivariate ANOVA-based methods (ANOVA simultaneous component analysis—ASCA, regularized MANOVA–rMANOVA, and Group-wise ANOVA-simultaneous component analysis—GASCA) in the framework of metabolomics studies. Our main goals are to compare these various ANOVA-based approaches and evaluate their performance on experimentally designed metabolomic studies to find the significant factors and identify the most relevant variables (potential markers) from the obtained results. Two experimental data sets were generated employing liquid chromatography coupled to mass spectrometry (LC-MS) with different complexity in the design to evaluate the performance of the statistical approaches. Results show that the three considered ANOVA-based methods have a similar performance in detecting statistically significant factors. However, relevant variables pointed by GASCA seem to be more reliable as there is a strong similarity with those variables detected by the widely used partial least squares discriminant analysis (PLS-DA) method.

Adverse Effects of Arsenic Uptake in Rice Metabolome and Lipidome Revealed by Untargeted Liquid Chromatography Coupled to Mass Spectrometry (LC-MS) and Regions of Interest Multivariate Curve Resolution

Rice crops are especially vulnerable to arsenic exposure compared to other cereal crops because flooding growing conditions facilitates its uptake. Besides, there are still many unknown questions about arsenic’s mode of action in rice. Here, we apply two untargeted approaches using liquid chromatography coupled to mass spectrometry (LC-MS) to unravel the effects on rice lipidome and metabolome in the early stages of growth. The exposure is evaluated through two different treatments, watering with arsenic-contaminated water and soil containing arsenic. The combination of regions of interest (ROI) and multivariate curve resolution (MCR) strategies in the ROIMCR data analyses workflow is proposed and complemented with other multivariate analyses such as partial least square discriminant analysis (PLS-DA) for the identification of potential markers of arsenic exposure and toxicity effects. The results of this study showed that rice metabolome (and lipidome) in root tissues seemed to be more affected by the watering and soil treatment. In contrast, aerial tissues alterations were accentuated by the arsenic dose, rather than with the watering and soil treatment itself. Up to a hundred lipids and 40 metabolites were significantly altered due to arsenic exposure. Major metabolic alterations were found in glycerophospholipids, glycerolipids, and amino acid-related pathways.

Yeast metabolic innovations emerged via expanded metabolic network and gene positive selection

Yeasts are known to have versatile metabolic traits, while how these metabolic traits have evolved has not been elucidated systematically. We performed integrative evolution analysis to investigate how genomic evolution determines trait generation by reconstructing genome-scale metabolic models (GEMs) for 332 yeasts. These GEMs could comprehensively characterize trait diversity and predict enzyme functionality, thereby signifying that sequence-level evolution has shaped reaction networks towards new metabolic functions. Strikingly, using GEMs, we can mechanistically map different evolutionary events, e.g. horizontal gene transfer and gene duplication, onto relevant subpathways to explain metabolic plasticity. This demonstrates that gene family expansion and enzyme promiscuity are prominent mechanisms for metabolic trait gains, while GEM simulations reveal that additional factors, such as gene loss from distant pathways, contribute to trait losses. Furthermore, our analysis could pinpoint to specific genes and pathways that have been under positive selection and relevant for the formulation of complex metabolic traits, i.e. thermotolerance and the Crabtree effect. Our findings illustrate how multidimensional evolution in both metabolic network structure and individual enzymes drives phenotypic variations.

The dynamics and role of sphingolipids in eukaryotic organisms upon thermal adaptation

All living beings have an optimal temperature for growth and survival. With the advancement of global warming, the search for understanding adaptive processes to climate changes has gained prominence. In this context, all living beings monitor the external temperature and develop adaptive responses to thermal variations. These responses ultimately change the functioning of the cell and affect the most diverse structures and processes. One of the first structures to detect thermal variations is the plasma membrane, whose constitution allows triggering of intracellular signals that assist in the response to temperature stress. Although studies on this topic have been conducted, the underlying mechanisms of recognizing thermal changes and modifying cellular functioning to adapt to this condition are not fully understood. Recently, many reports have indicated the participation of sphingolipids (SLs), major components of the plasma membrane, in the regulation of the thermal stress response. SLs can structurally reinforce the membrane or/and send signals intracellularly to control numerous cellular processes, such as apoptosis, cytoskeleton polarization, cell cycle arresting and fungal virulence. In this review, we discuss how SLs synthesis changes during both heat and cold stresses, focusing on fungi, plants, animals and human cells. The role of lysophospholipids is also discussed.

ROIMCR: a powerful analysis strategy for LC-MS metabolomic datasets

BACKGROUND

The analysis of LC-MS metabolomic datasets appears to be a challenging task in a wide range of disciplines since it demands the highly extensive processing of a vast amount of data. Different LC-MS data analysis packages have been developed in the last few years to facilitate this analysis. However, most of these strategies involve chromatographic alignment and peak shaping and often associate each "feature" (i.e., chromatographic peak) with a unique m/z measurement. Thus, the development of an alternative data analysis strategy that is applicable to most types of MS datasets and properly addresses these issues is still a challenge in the metabolomics field.

RESULTS

Here, we present an alternative approach called ROIMCR to: i) filter and compress massive LC-MS datasets while transforming their original structure into a data matrix of features without losing relevant information through the search of regions of interest (ROIs) in the m/z domain and ii) resolve compressed data to identify their contributing pure components without previous alignment or peak shaping by applying a Multivariate Curve Resolution-Alternating Least Squares (MCR-ALS) analysis. In this study, the basics of the ROIMCR method are presented in detail and a detailed description of its implementation is also provided. Data were analyzed using the MATLAB (The MathWorks, Inc., www.mathworks.com ) programming and computing environment. The application of the ROIMCR methodology is described in detail, with an example of LC-MS data generated in a lipidomic study and with other examples of recent applications.

CONCLUSIONS

The methodology presented here combines the benefits of data filtering and compression based on the searching of ROI features, without the loss of spectral accuracy. The method has the benefits of the application of the powerful MCR-ALS data resolution method without the necessity of performing chromatographic peak alignment or modelling. The presented method is a powerful alternative to other existing data analysis approaches that do not use the MCR-ALS method to resolve LC-MS data. The ROIMCR method also represents an improved strategy compared to the direct applications of the MCR-ALS method that use less-powerful data compression strategies such as binning and windowing. Overall, the strategy presented here confirms the usefulness of the ROIMCR chemometrics method for analyzing LC-MS untargeted metabolomics data.