- 1Applied Metabolomics Research Group, IMIM—Institut Hospital del Mar d’Investigacions Mèdiques, Barcelona, Spain
- 2Signal and Information Processing for Sensing Systems, Institute for Bioengineering of Catalonia (IBEC), The Barcelona Institute of Science and Technology, Barcelona, Spain
- 3Department of Electronics and Biomedical Engineering, Faculty of Physics, University of Barcelona, Barcelona, Spain
- 4Biomarkers and Nutrimetabolomics Laboratory, Department of Nutrition, Food Science and Gastronomy, Faculty of Pharmacy and Food Sciences, Campus Torribera, University of Barcelona, Sant Coloma de Gramanet, Spain
- 5Food Innovation Network (XIA), Santa Coloma de Gramanet, Spain
- 6Institut de Recerca en Nutrició i Seguretat Alimentària (INSA-UB), Santa Coloma de Gramanet, Spain
- 7Centro de Investigación Biomédica en Red de Fragilidad y Envejecimiento Saludable (CIBERFES), Instituto de Salud Carlos III, Madrid, Spain
Introduction: There is evidence that sample treatment of blood-based biosamples may affect integral signals in nuclear magnetic resonance-based metabolomics. The presence of macromolecules in plasma/serum samples makes investigating low-molecular-weight metabolites challenging. It is particularly relevant in the targeted approach, in which absolute concentrations of selected metabolites are often quantified based on the area of integral signals. Since there are a few treatments of plasma/serum samples for quantitative analysis without a universally accepted method, this topic remains of interest for future research.
Methods: In this work, targeted metabolomic profiling of 43 metabolites was performed on pooled plasma to compare four methodologies consisting of Carr-Purcell-Meiboom-Gill (CPMG) editing, ultrafiltration, protein precipitation with methanol, and glycerophospholipid solid-phase extraction (g-SPE) for phospholipid removal; prior to NMR metabolomics analysis. The effect of the sample treatments on the metabolite concentrations was evaluated using a permutation test of multiclass and pairwise Fisher scores.
Results: Results showed that methanol precipitation and ultrafiltration had a higher number of metabolites with coefficient of variation (CV) values above 20%. G-SPE and CPMG editing demonstrated better precision for most of the metabolites analyzed. However, differential quantification performance between procedures were metabolite-dependent. For example, pairwise comparisons showed that methanol precipitation and CPMG editing were suitable for quantifying citrate, while g-SPE showed better results for 2-hydroxybutyrate and tryptophan.
Discussion: There are alterations in the absolute concentration of various metabolites that are dependent on the procedure. Considering these alterations is essential before proceeding with the quantification of treatment-sensitive metabolites in biological samples for improving biomarker discovery and biological interpretations. The study demonstrated that g-SPE and CPMG editing are effective methods for removing proteins and phospholipids from plasma samples for quantitative NMR analysis of metabolites. However, careful consideration should be given to the specific metabolites of interest and their susceptibility to the sample treatment procedures. These findings contribute to the development of optimized sample preparation protocols for metabolomics studies using NMR spectroscopy.
1 Introduction
Metabolomics refers to the investigation of low-molecular-weight compounds present in biological samples (Wishart, 2019). Liquid chromatography-tandem mass spectrometry (LC-MS) and nuclear magnetic resonance (NMR) are popular techniques to explore alterations in the metabolome (Guijas et al., 2018). Metabolomic studies applied to human research generally employ biological samples such as urine, blood, saliva, and feces, among others. Often, extraction of metabolites from different matrices, such as blood, is necessary (Vignoli et al., 2019).
Blood-based samples—plasma and serum—are key in human metabolomics investigations (Zhang et al., 2012). Blood is the primary carrier of small compounds in the human body. This matrix carries nutrients, hormones, dissolved gases, and breakdown products, and it is involved in regulating body temperature, pressure, pH stabilization, and defense system, among other functions (Martini and Ober, 2006). Plasma/serum components are dysregulated by the presence of dysfunctions in the organism and/or pathological conditions (Psychogios et al., 2011). This makes these biofluids appropriate samples for the study of mechanisms of action, disease characterization, and the discovery of biomarkers of a specific condition through comprehensive fingerprinting (Scalbert et al., 2014; Yan and Xu, 2018; Attard et al., 2019; Aderemi et al., 2021).
Plasma/serum samples are also rich in other nonpolar molecules such as cholesterol and triglycerides, that need to be transported in the biofluid associated with various lipoprotein particles. Lipoproteins have, in general, heavier molecular weight than polar metabolites. They are divided based on size, density, and their relative content of triglycerides, cholesterol, and protein into high-density lipoproteins (HDL), low-density lipoproteins (LDL), intermediate-density lipoproteins (IDL), very-low-density lipoproteins (VLDL), and chylomicrons (Cox and García-Palmieri, 1990). The plasma/serum lipid fraction is a subject of systematic study in another omics technology, lipidomics. The lipoprotein profile can be measured when lipids and lipid composition change as a result of pathological states or disorders (Li et al., 2017). One of the most recent and applicable examples is the COVID-19 disease (Bruzzone et al., 2020; Kimhofer et al., 2020; Meoni et al., 2021). However, the presence of these macromolecules implies a great challenge to exploring low-molecular-weight metabolites (Ulaszewska et al., 2019; Ulmer et al., 2021). For this latter purpose, prior extraction of the low-molecular-weight metabolites from this matrix is necessary to be able to detect and, if necessary, quantify such metabolites.
The International Organization for Standardization (ISO) has already published specifications for pre-examination processes for metabolomics in urine, venous blood serum, and plasma (ISO 23118, ISO, 2021); however, these recommendations are yet to be universally adopted. In a multicenter study performed by Ghini et al. (2022), they examined how different operating procedures, such as different collection tubes, processing time, and storage conditions, affected the level of metabolites in serum and plasma samples (Ghini et al., 2022). In addition to these differences in the stage of sample pre-processing, sample treatment also affects the final metabolite levels. The reliability of the sample treatment utilized for the extraction requires that the chemical nature and relative concentration of metabolites found in the plasma/serum samples stay unaltered by the technique. However, in NMR-based metabolomics, the consistent quantification of plasma/serum metabolites remains a difficult concern (Ghini et al., 2019; Gowda et al., 2021), finding differences among sample treatments (Nagana Gowda and Raftery, 2014). Proteins and phospholipids show a large number of mostly broad NMR peaks, which overlap with the resonances of low-molecular-weight metabolites. They introduce an unstable spectral baseline, which makes identification and quantification of these metabolites difficult, and in addition, proteins may interact with the peak intensity of specific metabolites (Tiziani et al., 2008; Nagana Gowda and Raftery, 2014).
One of the approaches addressed in NMR may be partially overcome by a spin-spin relaxation edition, using the Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence (Wang et al., 2003). Nevertheless, this approach has important drawbacks, such as the intensity loss of peak signals, the waste of relevant information about the chemical nature of phenolic conjugates, and the introduction of artifacts. Several metabolites such as tyrosine, histidine, and lactate, may bind to plasma/serum proteins, and consequently, the absolute concentration of several metabolites is undervalued and unreliable without cleaving the bonds to proteins (Bell et al., 1988; Nicholson and Gartland, 1989; Chatham and Forder, 1999). Similarly, trimethylsilylpropanoic acid (TSP), the quantification standard commonly utilized in aqueous solutions such as plasma/serum, displays the tendency to bind to protein (Wallmeier et al., 2017). Hence, the utilization of TSP as quantification standard in plasma/serum samples acquired by the CPMG sequence presents several challenges. Firstly, chemical shift variations can occur in TSP when using CPMG sequences due to pulse imperfections and relaxation differences during refocusing periods. These variations may deviate from the expected chemical shift of TSP at around 0 ppm, potentially introducing errors in quantification. Moreover, elimination of macromolecules that are bound to TSP clearly causes changes in TSP intensity and peak shape. Such interference compromises the accuracy and reliability of quantification using TSP as an internal standard (Gowda et al., 2021). Suitable experimental designs and appropriate data processing techniques may be required to overcome these challenges and improve the reliability in CPMG sequences. Exploring alternative internal standards may enhance the accuracy of quantification. An option mostly preferred for quantitative NMR are maleic, fumaric and formic acids, but to date has been poorly utilized in metabolic studies (Bliziotis et al., 2020; Gowda et al., 2021). On the other hand diffusion-weighted NMR pulse sequence has been also used to detect signals from only the macromolecular components, e.g., lipoproteins (Liu et al., 1996), and along with computational algorithms, it may be useful in quantitative metabolomics (de Graaf and Behar, 2003). Nevertheless, diffusion sensitization easily erases signals from fast-diffusing metabolites and slow-diffusing lipoproteins, requiring distinct diffusion coefficients for each chemical structure (de Graaf et al., 2015) and has not yet been tested for metabolites at concentrations below 10 μM. Several sample treatments have emerged to remove proteins and phospholipids from serum/plasma.
A more comprehensive range of metabolites can be quantified by ultrafiltration compared to CPMG by removing macromolecules with filters (Psychogios et al., 2011). However, due to its time-consuming filtration step that requires expensive filters, this gold standard has several limitations, including laboriousness and cost (Wallmeier et al., 2017). Moreover, likewise as for CPMG editing, the application of ultrafiltration may underestimate certain metabolites bound to macromolecules. Alternatively, precipitation with organic solvents have been primarily applied for quantifying plasma/serum metabolites before NMR analysis (Nagana Gowda and Raftery, 2014; Nagana Gowda and Raftery, 2017). A less popular technique is the addition of electrically charged silica nanoparticles, which can be utilized combined with ultrafiltration and solvent precipitation for the aggregation and co-precipitation of proteins (Zhang et al., 2016). Although this treatment exhibited spectral improvements at the qualitative level, the technique performance does not allow quantitative analysis by itself.
In MS, glycerophospholipid-SPE (g-SPE) has been described, and compared to ultrafiltration and solvent extraction with different gradients, finding that SPE-mediated phospholipid removal was the treatment with the best coverage of non-lipid metabolites, extraction reproducibility, and minimization of matrix effects (Tulipani et al., 2013; González-Domínguez et al., 2020). Although the SPE method has led to considerable improvements in MS-based metabolomics (Liu et al., 2019), this methodology has not yet been described for treating plasma/serum before quantitative NMR analysis. Considering that sample treatments in quantitative NMR plasma analysis may involve metabolite-dependent particularities, this topic remains of interest, especially for targeting a particular set of metabolites (Nagana Gowda and Raftery, 2017). This work aimed to evaluate the performance of four procedures to remove plasma proteins and phospholipids, CPMG, ultrafiltration, protein precipitation, and g-SPE in pooled human plasma. To this end, each sample treatment combined with 1H-NMR analysis enabled the quantitation of up to 40 metabolites, including amino acids, carboxylic acids, and short-chain fatty acids, among other compounds.
2 Materials and methods
2.1 Subjects
Venous blood samples from five non-smoking fasting healthy volunteers with 37.9 ± 10.3 years of age (mean ± SD) were collected into heparin-containing vials at the Hospital Clinic of Barcelona (Spain). Exclusion criteria included serious illness, supplement intake, medication, and pregnancy. Plasma samples were acquired after the removal of cells by centrifugation at 1,600 g for 15 min at room temperature. Pooled plasma aliquots were used to avoid biological variability in the comparative analysis among sample preparation procedures, isolating variability derived from the techniques. Aliquots were stored at −80°C for analysis and processing.
The Bioethical Committee of the University of Barcelona approved the research protocol, and all the participants provided written informed consent. The study was inscribed as ISRCTN17200423 in the International Standard Randomized Controlled Trial Number registry.
2.2 Sample treatment
The different sample treatments of plasma samples were conducted independently.
Protein precipitation with methanol. Four 200 μL of pooled plasma were thawed, vortexed, mixed with methanol in a 2:1 solvent-to-plasma ratio (v/v), and incubated at −20°C for 20 min. The samples were centrifuged at 13,400 g for 30 min to pellet macromolecules and proteins. Supernatants were decanted into new vials and dried under a nitrogen stream. Dried samples were mixed with 100 μL of phosphate buffer in D2O, which contained 2.32 mM of 3-(trimethylsilyl)-proprionate-2,2,3,3-d4 (TSP), and the pH was adjusted to 7.0.
Glycerophospholipid SPE. SPE procedure using Ostro® 96 plate (Ostro plates, Waters) for protein precipitation and phospholipid removal plates was adapted to NMR-based metabolomics spectroscopy based on a previously published procedure for LC-MS (Tulipani et al., 2013). Four 200 μL aliquots of pooled plasma were thawed, vortexed, and mixed with acidic solvent, followed by g-SPE with pressure valves. Samples were pipetted into the plate wells, followed by the forceful addition of 1% formic acid in acetonitrile 3:1 solvent-to-plasma ratio (v/v). After drying, samples were mixed with 100 μL of the phosphate buffer in D2O (pH 7.0) containing 2.32 mM of TSP.
Ultrafiltration. Ultrafiltration was applied based on a published methodology with modifications (Nagana Gowda and Raftery, 2014). Centrifugal filters (3 kDa cutoff; Amicon Microcon, YM-3; Sigma-Aldrich) were initially washed with distilled water. Additionally, filter tubes were centrifuged three times with 500 μL of distilled water, each time at 13,400 g for 20 min, to remove residual glycerol from filters. Four 200 μL of pooled plasma were then transferred to filter tubes and centrifuged as previously described. The filtrates were mixed with 100 μL of the phosphate buffer.
Untreated. Four 200 μL of pooled plasma were then transferred to filter tubes and centrifuged as previously described. The filtrates were mixed with 100 μL of the phosphate buffer, utilizing 0.3 μM of formic acid as internal standard.
All the solutions were individually made up to 600 μL with D2O and transferred to 5 mm NMR tubes.
2.3 NMR spectra acquisition
All NMR experiments were performed at 298 K on a Bruker Avance Neo 400 MHz spectrometer equipped with a cryoprobe. One-dimensional proton spectra were acquired using a pulse sequence (1D NOESY) with presaturation of the water resonance. Additionally, untreated plasma samples were acquired using a CPMG pulse sequence. For both sequences, an exponential window function was applied to the free induction decay (FID) with a line-broadening factor of 0.3 Hz prior to the Fourier transformation. For each spectrum, FIDs were collected into 32 K data points (128 scans) with a spectral width of 14 ppm, an acquisition time of 3 s, a relaxation delay of 5 s, and a mixing time of 10 ms. Parameters were adjusted to allow metabolite quantification by Chenomx NMR Suite 10 profiler (Chenomx Edmonton, Canada). All the NMR spectra were phased, baseline corrected, and calibrated (TSP, 0.0 ppm) using TopSpin software (version 3.6, Bruker, BioSpin, Germany).
2.4 NMR metabolic profiling and quantitation
Identification and quantitation of detected metabolites were performed using Chenomx NMR Suite 10 profiler. Chenomx is a robust, reliable, and widely spread software for NMR-based metabolite quantification (Jung et al., 2016; Crook and Powers, 2020). This software enables metabolite deconvolution in complex samples and can determine concentrations in overlapped spectral regions. The software fits the peaks with a set of model spectra characterizing the chemical environments of each metabolite. Finally, the software calculates the concentrations of each metabolite using the peak areas and the calibration curves generated from the model spectra.
The fit of spectral signals was performed with a standard metabolite library for 400 MHz 1H NMR spectra. For treated plasma samples, peak fitting with reference to the internal TSP signal allowed quantification of absolute concentrations for all identified metabolites. For untreated samples variations in previous formic acid were standardized (blank) and the formic acid signal (δ8.44, 35.05 Hz) was utilized for untreated plasma samples. Further contributions to the proton peak assignment were provided by comparing the chemical shifts with those available in the Human Metabolome Database (http://www.hmdb.ca).
2.5 Data analysis
The dataset was imported to R software version 4.2 for statistical analyses. The coefficient of variation (CV) was calculated to determine the replicability for each metabolite and those with a CV > 20% in more than one group were subsequently removed from the multiclass analysis. The differences in metabolite concentration associated with the different treatments were evaluated using the Fisher Score (Tang et al., 2014). In order to provide stable estimates of the statistical properties, permutation tests with random sampling with replacement. This method is often used when the sample size is compromised to increase the effective sample power (Dwivedi et al., 2017; Anderson and Verbeeck, 2023). We performed 5,000 multiclass permutations on Fisher Scores to distinguish their sensitivity to sample treatment. Original Fisher Scores were compared with Fisher Scores randomly distributed, and 95% interval confidence and p-values were determined from permutations. Pairwise assessments of significant metabolites were applied to test binary permutation test (5,000 binary permutations) of Fisher Scores. The Benjamini–Hochberg procedure was carried out on all analyses to control the false discovery rate (FDR) (Benjamini and Hochberg, 1995). An FDR-corrected p-value < 0.05 was considered statistically significant.
3 Result and discussion
This study focused on investigating the influence of different sample treatments on the concentrations of a panel of 43 metabolites found in human plasma (Supplemental Figure S1). Compounds, concentrations (mg/dL), and CV (%) are detailed in Table 1 and Figure 1. Based on the present results, methanol precipitation exhibited a number of six compounds with a CV > 20%, while ultrafiltration presented eight compounds with a CV above this threshold. The metabolites meeting this criteria in at least three of the methodologies corresponded with acetate, isobutyrate and N,N-dimethylglycine, and therefore, did not meet the benchmark for acceptable precision in analytical methods (Bioanalytical Method Validation Guidance for Industry FDA, 2023; European Medicine Agency, 2011; ICHH Guideline, 2019). Although acetate is a volatile short-chain fatty acid, it is generally considered stable for quantification under standard storage and handling conditions (Psychogios et al., 2011). However, acetate can be subject to degradation or changes in concentration over time if not properly handled or stored (Bernini et al., 2011). Factors such as temperature, pH, exposure to air or light, and the presence of enzymes or microorganisms can affect the stability of acetate (Brunius et al., 2017). As with acetate, isobutyrate, methanol, and compounds with low-molecular weight (<90 g/mol) such as 2-hydroxybutyrate and ethanol, and might be affected by the evaporation process of the sample. Remarkably, the average concentration of 2-hydroxybutyrate after ultrafiltration and CPMG was three times lower compared to g-SPE. The measure of this metabolite was uniquely precise in the g-SPE (0.240 ± 0.019 mg/dL, CV of 8.1%) and CPMG (0.090 ± 0.017 mg/dL, CV of 19.3) methodologies. On the other hand, N,N-dimethylglycine did not show a clear and distinguishable spectral signal after ultrafiltration, while in methanol precipitation, the presence of a baseline made its quantification difficult, and therefore, an increased error was expected. Similarly, 3-methyl-2-oxovalerate, isobutyrate and methionine neither showed a clear spectral signal in CPMG sequence, finding lipoprotein signals destabilizing the baseline. Citrate, an important tricarboxylic acid cycle intermediate, also displayed a high CV after g-SPE and ultrafiltration. This may be due to the long sample processing time, especially when applying ultrafiltration. Differential sample processing times can impact the enzymatic reactions occurring within the samples. In particular, plasma citrate, lactate and acetoacetate levels are known to decrease over time after sample collection, likely due to enzymatic degradation (Ghini et al., 2019). On the other hand, prolonged processing times can result in increased red blood cells lysis, thereby promoting the release of citrate synthase and subsequent conversion of glucose to citrate (Larsen et al., 2012). Brunius et al., also reported citrate as one of the features related to kinetic drift of NMR plasma metabolomics (Brunius et al., 2017). Citrate displayed a CV of 10.6% and 8.0% after the methanol precipitation and CPMG methodologies, respectively, which should be the techniques of choice for the quantification of this plasma metabolite. Despite g-SPE being a fast procedure, citrate was not generally detected. In line with these findings, Tiziani and co-workers found a lower concentration of citrate when a mix of methanol and chloroform was used for precipitating the plasmatic proteins (Tiziani et al., 2008), and Gowda et al. showed that citrate was also at higher concentration after ultrafiltration, compared with methanol precipitation (Nagana Gowda and Raftery, 2014).
TABLE 1. Compounds, absolute concentrations (mg/dL) and coefficient of variation (CV, %) in plasma pools (mean ± SD) using 1H-NMR spectroscopy. Metabolites are sorted alphabetically.
FIGURE 1. Concentrations of metabolites with a top-10 fisher score in any of the comparisons (extension in Supplemental Figure S1). PR, methanol precipitation; SPE, glycerophospholipid-solid phase extraction; UF, ultrafiltration.
Fisher scores of measured metabolites are presented in Figure 2. The permutation test of multiclass Fisher scores showed that the concentrations of 18 metabolites were modified among the sample treatments with an FDR-corrected p-value < 0.05, shown in Table 1. Significant metabolites were then tested on a pairwise basis, shown in Table 2. The results showed that the efficiency of analyte extraction depended on the metabolite and the methodology utilized. Pairwise analysis denoted that concentrations of tryptophan, ornithine and acetoacetate were different among nearly all the procedures. While compared to CPMG, the techniques of SPE and ultrafiltration appeared to display the greatest differences, methanol precipitation and ultrafiltration seemed to exhibit more similar results.
FIGURE 2. Fisher scores for multiclass and pairwise comparisons. Metabolites were sorted according to their median fisher criterion across the four comparisons. CPMG, Carr-Purcell-Meiboom-Gill; PR, methanol precipitation; SPE, glycerophospholipid-solid phase extraction; UF, ultrafiltration.
Overall, branched-chain amino acid related valine, 3-methyl-2-oxovalerate and 2-aminoisobutyrate, aromatic amino acids such as phenylalanine, and tyrosine, and other amino acids such as alanine, glutamine, lysine, threonine and ornithine, along with creatinine, were found at lower concentrations in CPMG compared with the other sample treatments. On the other hand, acetone and glycine were found at higher concentration in untreated samples followed by CPMG experiments compared to the three pre-analytical treatments. One possible explanation for the lower concentrations in CPMG experiments is the lack of additional steps for concentrating and isolating metabolites with weak or overlapping signals (Weljie et al., 2006). Ultrafiltration, for example, employs a selective membrane to remove high molecular weight proteins and retain smaller metabolites. This process effectively concentrates the target analytes, resulting in higher concentrations. Similarly, SPE utilizes specific sorbents to extract metabolites, enabling their enrichment and subsequent elution for analysis (Tulipani et al., 2015). Methanol precipitation involves the addition of methanol to the plasma samples, causing protein precipitation and facilitating the recovery of the supernatant containing metabolites. These concentration steps can contribute to higher metabolite levels compared to CPMG experiments, where no such concentration steps are involved (Weljie et al., 2006). Another issue in the use of CPMG pulse sequences involves a series of refocusing pulses to counteract the signal decay caused by transverse relaxation (T2) processes. This approach can be affected by metabolite-macromolecule interactions, leading to reduced signal intensities and consequently lower observed concentrations. Metabolites in plasma samples have the potential to interact with macromolecules such as proteins, lipoproteins, and other biomolecules (Barrilero et al., 2017). It is worth noting that the binding of metabolites to macromolecules can vary depending on factors such as metabolite structure, pH, ionic strength, and the composition of the sample matrix (Kirwan et al., 2018). Different metabolites exhibit varying affinities for macromolecules, leading to differential effects on their signal intensities and concentrations in untreated samples using CPMG experiments and treated with ultrafiltration. This binding can occur through electrostatic interactions, hydrogen bonding, or hydrophobic interactions. As a result, the signals originating from bound metabolites decay more rapidly during the relaxation period compared to the signals from freely mobile metabolites (Bliziotis et al., 2020). The reduction of signal intensities for bound metabolites during the time relaxation period can result in underestimation of their concentrations in CPMG experiments (Nagana Gowda et al., 2018). Since the concentration determination is based on the signal intensities, the reduced signals of bound metabolites contribute to the overall lower concentrations observed in plasma samples prepared using CPMG compared to methods that involve additional processing steps like SPE or methanol precipitation (Nagana Gowda and Raftery, 2014). This was particularly observed for certain aromatic amino acids such as tryptophan (0.039 ± 0.046 mg/dL) and phenylalanine (0.414 ± 0.072 mg/dL) in CPMG experiments. Likewise, a total absence of tryptophan signals was found after ultrafiltration while in methanol precipitation, tryptophan concentration was 0.551 ± 0.066 mg/dL and even higher after g-SPE (0.780 ± 0.011 mg/dL). Since most of the plasma tryptophan is bound to albumin in the blood (McMenamy and Oncley, 1958), other metabolites bound to protein, such as creatinine (Varshney et al., 2011; Wallmeier et al., 2017), lactate (Bell et al., 1988; Chatham and Forder, 1999), tyrosine and histidine (Nicholson and Gartland, 1989), may require a protein degradation, such as methanol precipitation or SPE, to be correctly quantified in NMR metabolomics analysis. Nevertheless, in our results, the high levels of lactate in both ultrafiltration and methanol precipitation compared to those of g-SPE and CPMG are remarkable whereas histidine was not detected after g-SPE but was found to be at regular concentration in CPMG. Creatinine (0.317 ± 0.036 mg/dL) and tyrosine (0.420 ± 0.064 mg/dL) were found at lower concentrations uniquely in untreated samples processed using CPMG sequence. A study by Gowda and colleagues exhibited higher concentration of citrate after the ultrafiltration of 300 µL of serum compared to methanol precipitation. In contrast with the present study, the authors were able to quantify tryptophan after the ultrafiltration (0.474 ± 0.004 mg/dL). These results might be due to the free fraction of the metabolite in a larger amount of serum utilized. Nevertheless, the authors also reported that tryptophan concentration was much lower in ultrafiltration compared to methanol precipitation (1.291 ± 0.049 mg/dL) (Nagana Gowda and Raftery, 2014). Compared to ultrafiltration, the ketones acetoacetate and 3-hydroxybutyrate decreased after g-SPE, while levels of phenylalanine, tryptophan, and threonine were higher in ultrafiltration. Considering that ultrafiltration is the longest methodology, lipolytic action through microbial degradation could explain the lower levels of ketone bodies in ultrafiltration at the expense of increased lactate. In a study on sample stability, Bernini et al. reported increased lactate, while glucose decreased within 2 h at room temperature after blood collection (Bernini et al., 2011). However, despite the methanol precipitation may deactivate bacterial degradation rapidly, levels of lactate after ultrafiltration and methanol precipitation were similar. In this regard, Pinto and colleagues observed the occurrence of enzymatic lipolytic action without microbial growth (Pinto et al., 2014). In our study, all samples were collected and immediately stored at −80°C, and the only time discrepancy occurred during sample treatment.
Plasma treatments specifically alter the concentrations of certain metabolites. For instance, the concentration of methanol, formate and glycerol was influenced by the sample treatments, which have modified concentrations of these metabolites inherent to the used procedures (Wishart, 2008; Ghini et al., 2019). A high concentration of methanol is related to a residue of the methanol used for the protein precipitation instead of a higher extraction of the endogenous metabolome. Although the mix with chloroform should improve the extraction of lipoproteins and lipids, the final ratio of methanol-plasma was set to 2:1 (v/v) (Nagana Gowda et al., 2015). Protein precipitation by using a 1:1 methanol-to-sample ratio retains a high level of residual proteins that complicates the identification/quantification (Snytnikova et al., 2019). Higher methanol-sample ratios may increase the spectral quality but also evaporation time and solvent residue, affecting sample stability (Crook and Powers, 2020). Likewise, a prominent singlet at δ8.49 corresponding to formate was present in the aromatic region of the spectra and it is related to a residue of the formate used during g-SPE treatment (Tulipani et al., 2013; González-Domínguez et al., 2020). Glycerol concentration was very elevated after ultrafiltration due to insufficient washing of filters (Ghini et al., 2019). These metabolites were then excluded from further analyses.
In the present study, the use of methanol precipitation retained ethanol and serine concentrations more effectively compared to g-SPE, while myo-inositol extraction was higher in the methanol precipitation and g-SPE groups. Other authors found a better recovery of most metabolites when g-SPE was compared to precipitation with methanol or ultrafiltration (Tulipani et al., 2013). The authors showed that the combination of solvent extraction and SPE-mediated removal of phospholipids, prior analysis by MS, was the most suitable sample preparation for detecting subtle quantitative changes for the majority of the remaining metabolites. The extraction of polar compounds such as acetylcholine, acetyl-l-carnitine, leucine, isoleucine, and phenylalanine was higher compared to ultrafiltration and solvent precipitation. With the exception of phenylalanine, none of these metabolites were statistically different across the methodologies used in the present study. Further optimization of the procedure may enhance the global extraction efficiency of this sample treatment. Interestingly, ethanol levels after ultrafiltration and CPMG were extremely unstable, and levels of this metabolite should be interpreted with caution when using these procedures. Many NMR-based metabolomics studies have used ultrafiltration due to its efficiency for protein removal with optimal metabolite extraction (Bathe et al., 2011; Psychogios et al., 2011; Farshidfar et al., 2012; Nagana Gowda and Raftery, 2014). Tiziani and colleagues reported that relative to the protein precipitation methods, ultrafiltration performs better in retaining metabolite concentrations (Tiziani et al., 2008). Nevertheless, this statement was discussed in other studies (Nagana Gowda et al., 2015). In the present study, several considerations have been encountered. After ultrafiltration, glycerol concentration was around 1.5 orders of magnitude higher than in methanol precipitation and g-SPE, despite that centrifugal filters were washed as published elsewhere (Nagana Gowda and Raftery, 2014). Further filter washing would critically increase sample handling and time preparation, and this may affect the stability of plasma/serum metabolites (Brunius et al., 2017).
In the present NMR metabolomics study, four methodologies for removing macromolecules in plasma samples, CPMG, ultrafiltration, methanol precipitation, and adapted g-SPE were evaluated. The results indicate that this step may be critical for the reliable quantitation of many metabolites, while the claimed reproducibility of NMR is also investigated. There are improvements in quantification performance between treatments that are specific to certain metabolites, finding differences that are metabolite-dependent. Each methodology also presents several considerations. For instance, the cost, the need to wash centrifugal filters several times critically increases the time procedure, and the loss of signal of protein-associated metabolites in plasma/serum samples are some of the limitations of ultrafiltration. On the other hand, in light of the present results, untreated plasma samples followed by CPMG experiments offer an attractive alternative in terms of processing time, cost and reliability of results, provided that limitations in the replicability of some metabolites and underestimation in the concentrations of several metabolites are considered. To mitigate the impact of metabolite-macromolecule binding, additional sample processing steps such as ultrafiltration, SPE, or methanol precipitation can be employed. These techniques aim to separate metabolites from macromolecules, thereby minimizing their interactions and allowing for more accurate quantification of metabolite concentrations. Although methanol precipitation is a popular sample treatment in plasma/serum samples, the extraction efficiency of methanol precipitation can depend on the solvent-to-sample ratio, and in our case, we have observed a baseline that made quantification difficult, whereas g-SPE provides spectra without macromolecule residues, and in terms of repeatability, it was the technique of choice. Yet, this treatment is also costly and presented various metabolites with higher concentrations while underestimating others. To conclude, this work highlights that the different methodologies affect extraction performance in a metabolite-dependent manner in quantitative metabolomics, and this allows the selection of the most appropriate sample treatment for each case.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving human participants were reviewed and approved by Bioethical Committee of the University of Barcelona. The patients/participants provided their written informed consent to participate in this study.
Author contributions
FM-G conducted the metabolomics analysis, and its interpretation, and wrote this manuscript. SO and SM conducted statistical analyses. OP conducted quantification analysis and edited this manuscript. CA-L and RL reviewed this manuscript and supervised this study. All authors contributed to the article and approved the submitted version.
Funding
This work was accomplished as part of the Food4Brain project and the coordinated project TargetML granted by the Spanish MINECO (PID2020-114921RB-C21 and PID2021-126543OB-C21, respectively) and with the support of the Fundació La Marató de TV3 (project 202123). The work also received funding from CIBERFES, funded by the Instituto de Salud Carlos III and co-funded by the European Regional Development Fund’s “A way to make Europe” and the Generalitat de Catalunya’s Agency AGAUR (2017SGR1546 and 2017SGR1721). F.M-G was supported by Grant FJC2018-035791-I funded by MCIN/AEI/10.13039/501100011033. C.A-L. is grateful for the ICREA Academia Award 2018. Additional financial support was provided by the Institut de Bioenginyeria de Catalunya (IBEC). IBEC is a member of the CERCA Programme/Generalitat de Catalunya.
Acknowledgments
CA-L is grateful for the ICREA Academia Award 2018. We would like to acknowledge the Departament d’Universitats, Recerca i Societat de la Informació de la Generalitat de Catalunya (expedient 2021 SGR 01393). The authors are also grateful for the Fundació La Marató de TV3 (project 202123).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2023.1125582/full#supplementary-material
References
Aderemi, A. V., Ayeleso, A. O., Oyedapo, O. O., and Mukwevho, E. (2021). Metabolomics: A scoping review of its role as a tool for disease biomarker discovery in selected non-communicable diseases. Metab 11, 418. doi:10.3390/METABO11070418
Anderson, W. N., and Verbeeck, J. (2023). Exact permutation and bootstrap distribution of generalized pairwise comparisons statistics. Math 11, 1502. doi:10.3390/MATH11061502
Attard, J. A., Dunn, W. B., Mergental, H., Mirza, D. F., Afford, S. C., and Perera, M. T. P. R. (2019). Systematic review: Clinical metabolomics to forecast outcomes in liver transplantation surgery. A J. Integr. Biol. 23, 463–476. doi:10.1089/OMI.2019.0086
Barrilero, R., Ramírez, N., Vallvé, J. C., Taverner, D., Fuertes, R., Amigó, N., et al. (2017). Unravelling and quantifying the “nMR-invisible” metabolites interacting with human serum albumin by binding competition and T2 relaxation-based decomposition analysis. J. Proteome Res. 16, 1847–1856. doi:10.1021/ACS.JPROTEOME.6B00814
Bathe, O. F., Shaykhutdinov, R., Kopciuk, K., Weljie, A. M., McKay, A., Sutherland, F. R., et al. (2011). Feasibility of identifying pancreatic cancer based on serum metabolomics. Cancer Epidemiol. Biomarkers Prev. 20, 140–147. doi:10.1158/1055-9965.EPI-10-0712
Bell, J. D., Brown, J. C. C., Kubal, G., and Sadler, P. J. (1988). NMR-invisible lactate in blood plasma. FEBS Lett. 235, 81–86. doi:10.1016/0014-5793(88)81238-9
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 57, 289–300. doi:10.2307/2346101
Bernini, P., Bertini, I., Luchinat, C., Nincheri, P., Staderini, S., and Turano, P. (2011). Standard operating procedures for pre-analytical handling of blood and urine for metabolomic studies and biobanks. J. Biomol. NMR 49, 231–243. doi:10.1007/S10858-011-9489-1/SCHEMES/3
Bliziotis, N. G., Engelke, U. F. H., Aspers, R. L. E. G., Engel, J., Deinum, J., Timmers, H. J. L. M., et al. (2020). A comparison of high-throughput plasma NMR protocols for comparative untargeted metabolomics. Metabolomics 16, 64. doi:10.1007/S11306-020-01686-Y
Brunius, C., Pedersen, A., Malmodin, D., Karlsson, B. G., Andersson, L. I., Tybring, G., et al. (2017). Prediction and modeling of pre-analytical sampling errors as a strategy to improve plasma NMR metabolomics data. Bioinformatics 33, 3567–3574. doi:10.1093/BIOINFORMATICS/BTX442
Bruzzone, C., Bizkarguenaga, M., Gil-Redondo, R., Diercks, T., Arana, E., García de Vicuña, A., et al. (2020). SARS-CoV-2 infection dysregulates the metabolomic and lipidomic profiles of serum. iScience 23, 101645. doi:10.1016/j.isci.2020.101645
Chatham, J. C., and Forder, J. R. (1999). Lactic acid and protein interactions: Implications for the NMR visibility of lactate in biological systems. Biochim. Biophys. Acta - Gen. Subj. 1426, 177–184. doi:10.1016/S0304-4165(98)00154-8
Cox, R. A., and García-Palmieri, M. R. (1990). Cholesterol, triglycerides, and associated lipoproteins. Boston: Butterworths.
Crook, A. A., and Powers, R. (2020). Quantitative NMR-based biomedical metabolomics: Current status and applications. Mol 25, 5128. doi:10.3390/MOLECULES25215128
de Graaf, R. A., and Behar, K. L. (2003). Quantitative 1H NMR spectroscopy of blood plasma metabolites. Anal. Chem. 75, 2100–2104. doi:10.1021/ac020782+
de Graaf, R. A., Prinsen, H., Giannini, C., Caprio, S., and Herzog, R. I. (2015). Quantification of 1H NMR spectra from human plasma. Metabolomics 11, 1702–1707. doi:10.1007/s11306-015-0828-1
Dwivedi, A. K., Mallawaarachchi, I., and Alvarado, L. A. (2017). Analysis of small sample size studies using nonparametric bootstrap test with pooled resampling method. Stat. Med. 36, 2187–2205. doi:10.1002/SIM.7263
Farshidfar, F., Weljie, A. M., Kopciuk, K., Buie, W. D., Maclean, A., Dixon, E., et al. (2012). Serum metabolomic profile as a means to distinguish stage of colorectal cancer. Genome Med. 4, 42. doi:10.1186/gm341
FDA (2023). Bioanalytical method validation guidance for Industry | FDA. Available at: https://www.fda.gov/regulatory-information/search-fda-guidance-documents/bioanalytical-method-validation-guidance-industry (Accessed April 4, 2023).
Ghini, V., Abuja, P. M., Polasek, O., Kozera, L., Laiho, P., Anton, G., et al. (2022). Impact of the pre-examination phase on multicenter metabolomic studies. N. Biotechnol. 68, 37–47. doi:10.1016/J.NBT.2022.01.006
Ghini, V., Quaglio, D., Luchinat, C., and Turano, P. (2019). NMR for sample quality assessment in metabolomics. N. Biotechnol. 52, 25–34. doi:10.1016/J.NBT.2019.04.004
González-Domínguez, R., Jáuregui, O., Queipo-Ortuño, M. I., and Andrés-Lacueva, C. (2020). Characterization of the human exposome by a comprehensive and quantitative large-scale multianalyte metabolomics platform. Anal. Chem. 92, 13767–13775. doi:10.1021/ACS.ANALCHEM.0C02008
Gowda, G. A. N., Raftery, D., and Hong, N. N. (2021). Evaluation of fumaric acid and maleic acid as internal standards for nmr analysis of protein precipitated plasma, serum, and whole blood. Anal. Chem. 93, 3233–3240. doi:10.1021/ACS.ANALCHEM.0C04766
Guijas, C., Montenegro-Burke, J. R., Warth, B., Spilker, M. E., and Siuzdak, G. (2018). Metabolomics activity screening for identifying metabolites that modulate phenotype. Nat. Biotechnol. 36, 316–320. doi:10.1038/nbt.4101
ICHH Guideline (2019). Bioanalytical method validation M10. Available at: https://www.wrib.org.
ISO (2021). Molecular in vitro diagnostic examinations - specifications for pre-examination processes in metabolomics in urine, venous blood serum and plasma. Berlin, Germany: Springer.
Jung, Y. S., Hyeon, J. S., and Hwang, G. S. (2016). Software-assisted serum metabolite quantification using NMR. Anal. Chim. Acta 934, 194–202. doi:10.1016/J.ACA.2016.04.054
Kimhofer, T., Lodge, S., Whiley, L., Gray, N., Loo, R. L., Lawler, N. G., et al. (2020). Integrative modeling of quantitative plasma lipoprotein, metabolic, and amino acid data reveals a multiorgan pathological signature of SARS-CoV-2 infection. J. Proteome Res. 19, 4442–4454. doi:10.1021/ACS.JPROTEOME.0C00519
Kirwan, J. A., Brennan, L., Broadhurst, D., Fiehn, O., Cascante, M., Dunn, W. B., et al. (2018). Preanalytical processing and biobanking procedures of biological samples for metabolomics research: A white paper, community perspective (for “precision medicine and pharmacometabolomics task group”—the metabolomics society initiative). Clin. Chem. 64, 1158–1182. doi:10.1373/CLINCHEM.2018.287045
Larsen, S., Nielsen, J., Hansen, C. N., Nielsen, L. B., Wibrand, F., Stride, N., et al. (2012). Biomarkers of mitochondrial content in skeletal muscle of healthy young human subjects. J. Physiol. 590, 3349–3360. doi:10.1113/JPHYSIOL.2012.230185
Li, J., Vosegaard, T., and Guo, Z. (2017). Applications of nuclear magnetic resonance in lipid analyses: An emerging powerful tool for lipidomics studies. Prog. Lipid Res. 68, 37–56. doi:10.1016/J.PLIPRES.2017.09.003
Liu, M., Nicholson, J. K., and Lindon, J. C. (1996). High-Resolution diffusion and relaxation edited one- and two-dimensional 1H NMR spectroscopy of biological fluids. Anal. Chem. 68, 3370–3376. doi:10.1021/AC960426P
Liu, X., Zhou, L., Shi, X., and Xu, G. (2019). New advances in analytical methods for mass spectrometry-based large-scale metabolomics study. Trac. Trends Anal. Chem. 121, 115665. doi:10.1016/J.TRAC.2019.115665
Martini, F., and Ober, W. C. (2006). Fundamentals of anatomy & physiology. 6th ed. San Francisco: Benjamin Cummings.
McMenamy, R. H., and Oncley, J. L. (1958). The specific binding of L-tryptophan to serum albumin. J. Biol. Chem. 233, 1436–1447. doi:10.1016/s0021-9258(18)49353-2
Meoni, G., Ghini, V., Maggi, L., Vignoli, A., Mazzoni, A., Salvati, L., et al. (2021). Metabolomic/lipidomic profiling of COVID-19 and individual response to tocilizumab. PLOS Pathog. 17, e1009243. doi:10.1371/JOURNAL.PPAT.1009243
Nagana Gowda, G. A., Djukovic, D., Bettcher, L. F., Gu, H., and Raftery, D. (2018). NMR-guided mass spectrometry for absolute quantitation of human blood metabolites. Anal. Chem. 90, 2001–2009. doi:10.1021/ACS.ANALCHEM.7B04089
Nagana Gowda, G. A., Gowda, Y. N., and Raftery, D. (2015). Expanding the limits of human blood metabolite quantitation using NMR spectroscopy. Anal. Chem. 87, 706–715. doi:10.1021/ac503651e
Nagana Gowda, G. A., and Raftery, D. (2014). Quantitating metabolites in protein precipitated serum using NMR spectroscopy. Anal. Chem. 86, 5433–5440. doi:10.1021/ac5005103
Nagana Gowda, G. A., and Raftery, D. (2017). Recent advances in NMR-based metabolomics. Anal. Chem. 89, 490–510. doi:10.1021/acs.analchem.6b04420
Nicholson, J. K., and Gartland, K. P. R. (1989). 1H NMR studies on protein binding of histidine, tyrosine and phenylalanine in blood plasma. NMR Biomed. 2, 77–82. doi:10.1002/nbm.1940020207
Pinto, J., Domingues, M. R. M., Galhano, E., Pita, C., Do Céu Almeida, M., Carreira, I. M., et al. (2014). Human plasma stability during handling and storage: Impact on NMR metabolomics. Analyst 139, 1168–1177. doi:10.1039/C3AN02188B
Psychogios, N., Hau, D. D., Peng, J., Guo, A. C., Mandal, R., Bouatra, S., et al. (2011). The human serum metabolome. PLoS One 6, e16957. doi:10.1371/journal.pone.0016957
Scalbert, A., Brennan, L., Manach, C., Andres-Lacueva, C., Dragsted, L. O., Draper, J., et al. (2014). The food metabolome: A window over dietary exposure. Am. J. Clin. Nutr. 99, 1286–1308. doi:10.3945/ajcn.113.076133
Snytnikova, O. A., Khlichkina, A. A., Sagdeev, R. Z., and Tsentalovich, Y. P. (2019). Evaluation of sample preparation protocols for quantitative NMR-based metabolomics. Metabolomics 15, 1–9. doi:10.1007/S11306-019-1545-Y/TABLES/3
Tang, J., Alelyani, S., and Liu, H. (2014). Feature selection for classification: A review. Data Classif. Algorithms Appl. 2014, 37–64. doi:10.1201/B17320
Tiziani, S., Emwas, A.-H., Lodi, A., Ludwig, C., Bunce, C. M., Viant, M. R., et al. (2008). Optimized metabolite extraction from blood serum for 1H nuclear magnetic resonance spectroscopy. Anal. Biochem. 377, 16–23. doi:10.1016/j.ab.2008.01.037
Tulipani, S., Llorach, R., Urpi-Sarda, M., and Andres-Lacueva, C. (2013). Comparative analysis of sample preparation methods to handle the complexity of the blood fluid metabolome: When less is more. Anal. Chem. 85, 341–348. doi:10.1021/ac302919t
Tulipani, S., Mora-Cubillos, X., Jáuregui, O., Llorach, R., García-Fuentes, E., Tinahones, F. J., et al. (2015). New and vintage solutions to enhance the plasma metabolome coverage by LC-ESI-MS untargeted metabolomics: The not-so-simple process of method performance evaluation. Anal. Chem. 87, 2639–2647. doi:10.1021/ac503031d
Ulaszewska, M. M., Weinert, C. H., Trimigno, A., Portmann, R., Andres Lacueva, C., Badertscher, R., et al. (2019). Nutrimetabolomics: An integrative action for metabolomic analyses in human nutritional studies. Mol. Nutr. Food Res. 63, 1800384. doi:10.1002/mnfr.201800384
Ulmer, C. Z., Maus, A., Hines, J., and Singh, R. (2021). Challenges in translating clinical metabolomics data sets from the bench to the bedside. Clin. Chem. 67, 1581–1583. doi:10.1093/CLINCHEM/HVAB210
Varshney, A., Rehan, M., Subbarao, N., Rabbani, G., and Khan, R. H. (2011). Elimination of endogenous toxin, creatinine from blood plasma depends on albumin conformation: Site specific uremic toxicity & impaired drug binding. PLoS One 6, e17230. doi:10.1371/JOURNAL.PONE.0017230
Vignoli, A., Ghini, V., Meoni, G., Licari, C., Takis, P. G., Tenori, L., et al. (2019). High-Throughput metabolomics by 1D NMR. Angew. Chem. Int. Ed. 58, 968–994. doi:10.1002/anie.201804736
Wallmeier, J., Samol, C., Ellmann, L., Zacharias, H. U., Vogl, F. C., Garcia, M., et al. (2017). Quantification of metabolites by NMR spectroscopy in the presence of protein. J. Proteome Res. 16, 1784–1796. doi:10.1021/ACS.JPROTEOME.7B00057
Wang, Y., Bollard, M. E., Keun, H., Antti, H., Beckonert, O., Ebbels, T. M., et al. (2003). Spectral editing and pattern recognition methods applied to high-resolution magic-angle spinning 1H nuclear magnetic resonance spectroscopy of liver tissues. Anal. Biochem. 323, 26–32.
Weljie, A. M., Newton, J., Mercier, P., Carlson, E., and Slupsky, C. M. (2006). Targeted profiling: Quantitative analysis of 1H NMR metabolomics data. Anal. Chem. 78, 4430–4442. doi:10.1021/AC060209G
Wishart, D. S. (2019). Metabolomics for investigating physiological and pathophysiological processes. Physiol. Rev. 99, 1819–1875. doi:10.1152/physrev.00035.2018
Wishart, D. S. (2008). Metabolomics: Applications to food science and nutrition research. Trends Food Sci. Technol. 19, 482–493. doi:10.1016/j.tifs.2008.03.003
Yan, M., and Xu, G. (2018). Current and future perspectives of functional metabolomics in disease studies–A review. Anal. Chim. Acta 1037, 41–54. doi:10.1016/J.ACA.2018.04.006
Zhang, A., Sun, H., Wang, P., Han, Y., and Wang, X. (2012). Recent and potential developments of biofluid analyses in metabolomics. J. Proteomics 75, 1079–1088. doi:10.1016/J.JPROT.2011.10.027
Keywords: metabolomics, nuclear magnetic resonance, plasma, quantitative analysis, quantification, sample treatment
Citation: Madrid-Gambin F, Oller S, Marco S, Pozo ÓJ, Andres-Lacueva C and Llorach R (2023) Quantitative plasma profiling by 1H NMR-based metabolomics: impact of sample treatment. Front. Mol. Biosci. 10:1125582. doi: 10.3389/fmolb.2023.1125582
Received: 16 December 2022; Accepted: 22 May 2023;
Published: 02 June 2023.
Edited by:
Justin J. J. Van Der Hooft, Wageningen University and Research, NetherlandsReviewed by:
Veronica Ghini, University of Florence, ItalyElena Chekmeneva, Imperial College London, United Kingdom
Copyright © 2023 Madrid-Gambin, Oller, Marco, Pozo, Andres-Lacueva and Llorach. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Francisco Madrid-Gambin, Zm1hZHJpZEBpbWltLmVz; Rafael Llorach, cmFmYWxsb3JhY2hAdWIuZWR1