1. Introduction
Lichens, fascinating symbiotic organisms composed of a photosynthetic partner (green algae or cyanobacteria) and a fungal partner, have been utilized for centuries in various fields, including medicine and industry. They serve as a vital food source for both animals and humans [
1]. Lichens produce a plethora of primary and secondary metabolites. While primary metabolites are essential for proper growth of lichens, secondary metabolites protect plants against environmental factors [
2]. Lichens contain 800–1050 secondary metabolites belonging to various groups, including aliphatic acids, anthraquinones, phenolic compounds, quinones, pulvinic acid derivatives, steroids, terpenes, and xanthones [
2,
3]. Thus, they are treasure troves of natural compounds with pharmacological potential, encompassing antimicrobial, antiproliferative, cytotoxic, and antioxidant effects [
4,
5,
6,
7,
8].
Among these compounds, gyrophoric acid (GA) stands out as a remarkable natural substance produced by several lichen species, prominently found in 31 out of 33 lichen species within the genus
Umbilicaria [
9]. It has a wide range of positive biological properties, including anticancer and antimicrobial effects, and is used in industrial sectors, such as the production of dyes [
6,
10,
11,
12]. Chemically, GA is a polyphenolic depside with three aromatic rings of orselinic acid linked by ester bonds [
13]. The aromatic rings in these phenolic compounds allow them to scavenge free radicals, giving GA its antioxidant properties. Moreover, these aromatic rings form a versatile structure that, along with hydroxyl and methyl groups, can anchor to the active sites of various enzymes and receptors. GA inhibits topoisomerase 1 activity in mammalian cells, leading to DNA damage and activation of the p53/p21 DNA damage pathway. It also inhibits the function of various other enzymes such as eukaryotic protein tyrosine phosphatase (PTP1B), zinc-dependent metalloprotease neprilysin (MME), bacterial urease, and enzymes responsible for glycation control [
9].
The first study describing its cytotoxic activity was published in 2004 [
14]. Since then, numerous studies have investigated the biological activities of GA in vitro, confirming its antiproliferative and cytotoxic effects in MM98, A431, HaCaT [
15] A2780, HeLa, MCF-7, SK-BR-3, HT-29, HCT-116 p53(+/+), HCT-116 p53(−/−), HL-60, and Jurkat cells [
16] or in 3D spheroid cells [
17], or antibacterial and antifungal effects [
17,
18], or antioxidant effects [
16,
18,
19,
20,
21,
22,
23].
Despite the wealth of knowledge regarding the biological activities of GA, no in vivo study has been conducted thus far. Therefore, the objective of our research was to assess the fundamental physicochemical and pharmacokinetic properties of GA and to carry out the first in vivo study on healthy animals.
3. Discussion
Over the past five years, more than 60 scientific publications have been published exploring the potential biological effects of GA. However, none of these studies have been conducted under in vivo conditions. Our study, therefore, represents the first investigation into the effects of GA in vivo. The objective of our research was to analyze the fundamental behavioral patterns, hematological and biochemical parameters, and blood metabolome of healthy Sprague-Dawley rats. Our findings reveal that GA is a molecule with low solubility, particularly in alkaline environments, yet it remains stable across different pH conditions. It has the ability to bind to HSA, potentially facilitating its transportation to target tissues. The half-life of GA was observed to be approximately 33 min. Administering a dose of 10 mg/kg of GA did not affect the food and fluid intake, nor did it impact body mass gain of experimental animals. However, we did observe changes in certain behavioral patterns. Specifically, in the OFT, the rats exhibited increased rearing activity. Additionally, in the EPM, apart from increased rearing activity, the rats demonstrated heightened center crossings and time spent in the open arms of the maze. Concomitantly, GA promotes proliferative activity of hippocampal progenitors, indicating its neurogenic potential.
Drug delivery into the organism involves various administration routes, highlighting the importance of assessing the stability of xenobiotic under different pH conditions. As early as 1956, GA was recognized as a highly pH-sensitive molecule [
36]. Nevertheless, our findings indicate that GA exhibits stability across a range of pH conditions.
HSA, as an abundant plasma protein, has garnered significant interest in the pharmaceutical industry due to its ability to bind a wide range of drugs, influencing their delivery and efficacy, and ultimately impacting the drugs’ pharmacokinetic and pharmacodynamics properties [
37]. Our findings indicate that spontaneous hydrophobic interactions play a significant role in the formation of the HSA–GA complex. Moreover, our results demonstrate that GA exhibits optimal affinity for HSA, which is crucial for the distribution of xenobiotic through plasma proteins and their subsequent release in targeted tissues. It is likely that GA binds to a specific binding site on HSA, with docking analysis data suggesting that GA binds in proximity to the TRP214 binding site on HSA.
The synchronous spectra obtained from studying the interaction of GA with HSA suggest that there are no significant changes in the microenvironment polarity near the Trp and Tyr residues. However, the 3D fluorescence spectra present contrasting results. This disparity may be attributed to the nature of the 3D measurement, as the prolonged measurement time could potentially induce structural changes in GA within the aqueous environment, consequently affecting the measurement results.
Because of the lack of available GA amount, we were not able to find lethal and semi-lethal concentrations of this compound, which seems to be a limitation of our study. Due to the absence of prior in vivo studies involving GA, we opted to utilize a similar dosage that has been employed for other secondary metabolites, such as atranorin [
38,
39], or usnic acid [
40]. Accordingly, GA was administered orally at a concentration of 10 mg/kg body weight on a daily basis for a duration of one month, as outlined in the Materials and Methods section. To assess potential health-related changes, we measured food and fluid intake, as well as body mass gain, which are crucial indicators. Due to a number of studies dealing with ethanol consumption, we missed out the control group treated with ethanol as a solvent, which is a limitation of our study.
Our findings indicate that GA had no discernible impact on food and fluid intake. Furthermore, there were no significant differences in body mass gain between the GA-treated group and the control group. Based on these results, we can infer that GA does not adversely affect body weight gain in both sexes.
Concomitantly with these results, we saw no substantial effects of GA on hematological parameters. The only noticeable parameters were depressed MON and GRA levels without any change of total count of WBC. It is important to note that a decrease in MON and GRA levels within the reference range is not considered clinically relevant. In healthy animals, these values can be as low as zero. Therefore, our results suggest that GA does not exert any negative impact on blood parameters, reinforcing its potential as a safe candidate in the field of herbal medicine.
From biochemical parameters, ALT was decreased. ALT is generally recognized as a specific indicator of hepatocellular damage. Elevated ALT activity typically signifies primary liver damage, which can be caused by factors such as drug-induced injury or secondary reactive hepatopathy associated with other organ or systemic diseases [
41]. On the other hand, a decrease in ALT levels may be a promising indicator of liver parenchymal regeneration, especially if another drug with potential hepatotoxicity is used in therapy together with GA. A study of Park et al. (2019) revealed that low ALT was a significant independent predictor for liver-related mortality [
42]. Additionally, two other parameters, urea and creatinine, showed decreased levels after GA treatment. Furthermore, the blood metabolome analysis revealed decreased levels of arginine, indicating the involvement of the urea cycle. Moreover, lipids from the phosphatidylcholine, lysophosphatidylcholine, and sphingomyelin groups were affected. These various mechanisms suggest ongoing processes within the liver. As GA is insoluble in majority of safe solvents, we have chosen ethanol as it has been shown to be safe at applied concentration [
43,
44]. Despite all of those facts, the liver remained histologically unchanged, suggesting the possible protective antioxidant effect of GA.
Considering the minimal alterations in hematological and biochemical parameters, we can conclude that administration of GA in the dose of 10 mg/kg to rats for a one-month period of time is safe and does not cause any changes in laboratory parameters that could be termed as an adverse effect of the therapy.
In behavioral tests, we monitored the animals to assess levels of various forms of behavior, such as locomotor and exploratory activity, comfort behavior, and anxiety. OFT is a test used to evaluate exploratory behavior and locomotor activity, as well as anxiety behavior in rodents [
45]. In the OFT, we observed an increased frequency of rearing in both experimental groups compared to the INT groups. Exploratory rearing is one of the markers to evaluate the level of exploration of a novel environment [
46]. Other parameters in this test, such as defecation, washing activity, time in the center, time in the periphery, travelled distance, and average speed, followed the INT groups analogously. Concomitantly, an increased rearing frequency after GA administration was observed in EPM. In addition to the aforementioned behaviors, we also observed an increased number of crosses through the center of the maze and an extended duration spent in the open arms of the maze. These observations suggest reduced anxiety-like behavior and increased exploratory tendencies in the tested animals. The travelled distance is an expression of locomotor activity in rodents [
47], while time spent in open arms, frequency of defecation, and frequency of rearing are expressions of the level of depression in rodents [
48]. Based on the EPM, we can evaluate the level of anxiety in rodents, which increased with decreasing time spent in the open arms of the maze by rodents and with lower number of passes through the center of the maze [
49]. All those forms of behavior give indices that GA may depress anxiety behavior of the animals and sustains the natural exploration.
Because the changes in behavior are associated with the changes in brain structures, we focused also on the permeability coefficient of blood-brain-barrier transport. As shown, the permeability coefficient was 0.037944 ± 0.009511 × 10
−3 cm/min. This value is very low as compared with for example propranolol [
50] or sucrose [
51]. As indicated before, the compounds with permeability coefficient *10
−6 cm/s > 1.5 are classified as high permeation predicted, with Pe *10
−6 cm/s < 1.5 classified as low permeation predicted [
52]. For example, physodic acid from the extract of
Hypogymnia physodes shows a permeability coefficient regarded as high (>1.5 × 10
−6 cm/s) [
53].
4. Materials and Methods
4.1. pH Properties and Stability
A 2 mM stock solution of GA (provided by RNDr. Michal Goga, PhD., UPJŠ Košice) was prepared by dissolving it in 96% ethanol and storing it in a dark environment at 8 °C. To create the working solution of GA, 5 μL of the stock solution was added to 2000 μL of distilled water. The stability of the GA solution in different pH conditions (acidic, neutral, and alkaline) was assessed by conducting repeated UV-VIS spectral measurements over a period of time. The pH measurements were conducted using a Hanna Instruments HI 2211 pH/ORP meter (Hanna Instruments, Woonsocket, RI, USA) equipped with a Hanna Instruments HI 1131 B electrode. UV-VIS spectra were acquired using a Specord S300 UV-VIS spectrometer (Analytik Jena AG, Jena, Germany) within the wavelength range of 184.5–550 nm. A quartz cuvette with a 1 cm light path was utilized for this purpose. The pH-related characteristics were investigated by observing pH-dependent alterations in the UV-VIS spectra of the GA aqueous solution, which was continuously agitated. The GA concentration in the working solution was maintained at 5 μM, and the spectra were recorded at a temperature of 22 °C. Initially, the aqueous GA solution underwent acidification through the addition of concentrated hydrochloric acid. The pH of the solution was adjusted to 3.8 initially. Subsequently, small increments of a 4 M NaOH aqueous solution were systematically introduced into the working solution within a cuvette, and both the pH and UV-VIS spectra were recorded. The titration process was concluded when the pH reached 10.1. Data obtained at wavelengths 239 nm, 308 nm, and 332 nm were subjected to analysis using GraphPad Prism 7.00 software. The pH curves were fitted using the sigmoidal standard curve model (Sigmoidal, 4PL, Warsaw, Poland) with X representing log(Concentration). To validate the results, the experimentally acquired data were compared with GA pKa values calculated using the Chemicalize.com web service, provided by © 2021 ChemAxon.
4.2. Microsomal Stability
The microsomal stability assay for insoluble compounds was performed as described by Di et al. [
54]. Briefly, GA was solubilized in DMSO (CAS no.67-68-5, ≥99.8% purity, Sigma Aldrich, St. Louis, MO, USA). Concentrated rat liver microsomes (Thermo Fisher Scientific, Waltham, MA, USA) were mixed with pre-warmed phosphate buffer (100 mM), and GA was added. The reaction was terminated by addition of 20 mM NADPH (Merck Life Science, Bratislava, Slovakia). The mixture was incubated at 37 °C for up to 120 min, terminated by the addition of acetonitrile, and centrifuged (30,000×
g/10 min). The supernatant was transferred into vials and analyzed using LC-MS/MS.
4.3. Ultra-Performance Liquid Chromatography Coupled to Tandem Mass Spectrometry
A Waters (Waters, Prague, Czechia) triple-quadrupole (Quattro premier XE) mass spectrometer with an electrospray ionization source in positive mode (ESI+) was used. The mass spectrometer was operated with the following parameters: capillary voltage, 3 kV; source temperature, 150 °C; desolvation temperature, 350 °C; cone gas, 50 L/h; and desolvation gas, 800 L/h. The source cone voltages and collision energies were manually optimized for each SRM transition. The chromatographic apparatus consisted of a Waters ACQUITY UPLC system with a binary gradient pump, autosampler, and column thermostat (Waters, Czechia). Chromatographic separation was performed on an ACQUITY UPLC BEH C18 (Thermo Fisher Scientific, USA) (2.1 × 50 mm, 1.7 μm particle size) column. The column temperature was maintained at 30 °C. Mobile phase A consisted of 0.1% formic acid in water. Mobile phase B consisted of 100% acetonitrile. The elution started at 20% B (0–0.5 min), increased to 90% B (0.5–2 min) to 90% B (2–3 min), returned to 20% B, and re-equilibrated at 3–4 min. The flow rate was 0.5 mL/min, and the injection volume was 5 μL. MassLynx software version 4.1 was used for instrument control, data acquisition, and analysis.
4.4. Binding to Human Serum Albumin
GA was dissolved in DMSO (CAS no.67-68-5, ≥99.8% purity, Sigma Aldrich, USA) to obtain a stock solution with a concentration of 0.7 µM. The stock solution was stored at −20 °C in the dark. Albumin from human serum (A1887, CAS no. 70024-90-7) and reagents for buffer solutions were obtained from Sigma Aldrich (USA).
The HSA stock solution was prepared by dissolving 20 × 10
−3 g/mL in 150 × 10
−3 M Tris-HCl-NaCl (pH 7.4, 22 °C). The concentration of the stock solution was determined by UV-Vis spectroscopy (Specord S300 UV VIS, Analytik Jena) using a quartz cuvette with an optical length of 1 cm. The extinction coefficient of HSA used in the calculation was ε280 = 35,700 M
−1cm
−1 [
55]. The HSA stock solution was stored at 4–8 °C in the dark.
Fluorescence quenching spectra were acquired using a Varian Cary Eclipse spectrofluorimeter (Varian, Inc., Palo Alto, CA, USA) within the wavelength range of 290–500 nm. A quartz cuvette with a 1 cm optical path length was employed for the measurements. The excitation slit width was set to 5 nm, and the emission slit width was adjusted to 10 nm. An excitation wavelength of 280 nm was used, with an averaging time of 0.5 s and a data interval of 1 nm. The spectra were recorded at constant human serum albumin (HSA) concentrations (c = 4.0 × 10−6 M) while varying the concentrations of Glycyrrhetinic acid (GA) in the range of 0.4–3.5 × 10−6 M. These measurements were conducted in a 10 × 10−3 M Tris-HCl buffer solution at pH 7.4 and a temperature of 22 °C. The experiments were repeated at three different temperatures: 15 °C, 25 °C, and 37 °C. Fluorescence quenching spectra were obtained using a Varian Cary Eclipse spectrofluorimeter (Varian, Inc., USA) in the wavelength range of 290–500 nm using a quartz cuvette with an optical length of 1 cm. The excitation slit width was 5 nm, and the emission slit width was set at 10 nm. An excitation wavelength of 280 nm, averaging time of 0.5 s, and 1 nm data interval were used. Spectra were measured at fixed HSA concentrations (c = 4.0 × 10−6 M) and increasing concentrations of GA (0.4–3.5 × 10−6 M) in 10 × 10−3 M Tris-HCl buffer solution (pH 7.4, 22 °C). The measurements were performed at 15, 25, and 37 °C.
Synchronous fluorescence spectra were captured at two distinct scanning intervals, ∆λ, defined as ∆λ = λem − λex. The excitation wavelength ranged from 200 to 350 nm, and fluorescence emission was recorded at ∆ = 15 nm for tyrosine (Tyr) residues and ∆ = 60 nm for tryptophan (Trp) residues. The concentrations of human serum albumin (HSA) and the ligands matched those employed in the fluorescence spectral measurements. All measurements were conducted at a temperature of 25 °C.
The GA molecule was subjected to molecular docking within the binding site of HSA, originally occupied by Thyroxine. This binding site corresponds to the configuration observed in the original crystal structure with PDB identification number 1HK1 [
56].
4.5. Experimental Design of In Vivo Experiments
Twenty-four Sprague-Dawley rats (Velaz, Czechia) were used in the experiment. The animals were maintained under standard vivarium conditions with a light regimen of 12:12 h, stable temperature (21–24 °C), and humidity (50–65%). The rats had free access to water and were fed with standard pellets (Velaz, Czechia) ad libitum. The animals were handled according to the guidelines established by Law Nos. 377 and 436/2012 of the Slovak Republic for the Care and Use of Laboratory Animals (approval numbers: Ro-2866/16-221, Ro-2219/19-221/3) and in accordance with ARRIVE guidelines. The animals were divided into four experimental groups. Each group consisted of six animals. The animals were divided according to sex into two control intact healthy groups (INT males and INT females) and two groups with daily administration of GA (GA males and GA females). GA at a dose of 10 mg/kg body weight (dissolved in 10% ethanol solution) was administered orally for a period of one month. Control rats were left undisturbed throughout the entire protocol, except for blood collection, behavioral testing, and regular cage cleaning, which was done for all animals by the same experimenter. Blood was collected from the tail vein before and after treatment for hematological analyses. Whole blood in K3EDTA tubes was analyzed using an automated veterinary hematology analyzer (MindrayBC 2800VET, Guangzhou, China). At the end of the experiment, blood was collected again (for biochemical and metabolomics analyses). Serum samples were stored at 4 °C, and biochemical analysis was performed over two consecutive days using an automated clinical chemistry analyzer (Cobas c 111, Roche, Basel, Switzerland) according to the manufacturer’s instructions (
Figure 8).
4.6. Behavioral Analyses
The conduct of an organism mirrors its mental and physical condition. To assess animal behavior, two distinct behavioral assessments were employed to observe fundamental exploratory conduct: movement (open field test (OFT)) and levels of anxiety (elevated plus-maze (EPM)).
During the open field test (OFT), various parameters were measured, including locomotor activity and exploration. These measurements encompassed the total distance covered and the average speed of the animals, as well as the frequency of rearing. Additionally, anxiety levels were assessed by recording the number of center crossings, the time spent in the center versus the periphery, the number of defecation boluses, and the frequency of washing. At the commencement of each trial, the animals were positioned in the central area of the testing field, and each trial had a duration of 6 min. Throughout the test, the animals were continuously recorded using video footage, and the subsequent analysis of the results was conducted using Smart Junior software (ver. 3.0) (Panlab, Barcelona, Spain).
During the elevated plus-maze (EPM) test, various parameters were recorded, including locomotor activity (quantified by the number of center crossings), exploratory activity (measured by rearing frequency), and the level of anxiety (evaluated through time spent in open arms, the number of defecation events, and washing frequency). To initiate the test, the animals were positioned in the central square of the maze, and each trial had a duration of 6 min. For a detailed description, see [
38].
4.7. Metabolomics
A total volume of 100 µL of blood was obtained from the tail vein at a single time point and collected into microtubes. The collected blood serum was then preserved by storing it at −80 °C. After thawing, the metabolites were analyzed using liquid chromatography-tandem mass spectrometry (LC-MS/MS) with the assistance of a Biocrates AbsoluteIDQ p180 kit (Biocrates Life Sciences AG, Innsbruck, Austria). Chromatographic separation was accomplished using an ultra-performance liquid chromatography system (ACQUITY UPLC I-Class, Waters, Czechia). This chromatography system was coupled with a triple quadrupole mass spectrometer (XEVO TQ-S, Waters, Czechia). Metabolites were extracted and quantified following the manufacturer’s guidelines, as previously [
27,
57,
58]. The values were subjected to a log2 transformation to achieve normally distributed data and to stabilize the variance, consistent with the previous description [
38].
4.8. The Blood–Brain Barrier Transport
Rat-mixed glial cell culture was prepared from the cerebral cortices of 0–2-day-old Sprague-Dawley rats as described previously [
59]. Briefly, cerebral cortices were dissected, stripped of the meninges, and mechanically dissociated by repeated pipetting, followed by passage through a 20 µm nylon mesh. Cells were plated in 6-well plates pre-coated with poly-l-lysine (10 µg/mL; Sigma-Aldrich, USA) and cultivated in DMEM containing 10% FCS and 2 mM ultraglutamine (LONZA, Basel, Switzerland) at 37 °C and 5% CO
2 in a humidified atmosphere. For further experiments, three-week-old glial cultures were used.
Primary rat brain endothelial cells were isolated as described previously [
60]. Briefly, Sprague-Dawley rats were euthanized, and their whole brains were removed. Under sterile conditions, the brainstem and cerebellum were dissected, and the midbrain, white matter, and choroid plexus were removed. The remaining cortical tissue was removed from the meninges using dry Whatman paper. Tissue was homogenized on ice in DMEM-F12 (PAA Laboratories, Cölbe, Germany). Brain homogenate was centrifuged at 800×
g/10 min at 4 °C. The supernatant was aspirated, and the pellet was re-suspended in pre-warmed digestion mix containing collagenase/dispase (Roche Diagnostics, Indianapolis, IN, USA) and DNase I (Sigma-Aldrich, USA). The tissue was incubated with the prepared digestion mix at 37 °C for 30 min with gentle shaking. Digested tissue was centrifuged at 800×
g/10 min at 4 °C, and the pellet was re-suspended in 20% bovine serum albumin (BSA, Sigma-Aldrich, USA). The tissue was centrifuged at 1500×
g/15 min at 4 °C to obtain pellets containing microvessels with a fraction of myelin and BSA at the top, which was centrifuged again. The microvessels were pooled, re-suspended in a pre-warmed digestion mix, and digested for 15 min at 37 °C. The microvessel pellet was again centrifuged at 800×
g/10 min at 4 °C and washed with a culture medium containing serum. The microvessels were cultured in DMEM-F12 medium (PAA Laboratoris, Germany) containing 15% plasma-derived serum (PDS) (First Link, Birmingham, UK), 2 mM L-glutamine (GE Healthcare, Buckinghamshire, UK), BME vitamins (Sigma–Aldrich, USA), heparin (Sigma-Aldrich, USA), and 3 µM puromycine (Sigma-Aldrich, USA). After 7 days, cells were seeded on the upper chamber (apical compartment) of Transwell insert with a 0.4 µm pore size (Becton Dickinson, Franklin Lakes, NJ, USA) pre-coated with 10 µg/cm
2 collagen type IV (Sigma-Aldrich, USA) and 5 µg/cm
2 fibronectin (Sigma-Aldrich, USA). Endothelial cells were cultivated together with mixed glial cells cultivated as an adherent monolayer in the basolateral compartment of 12-well Cortar plates for further seven days in EBM-2 medium containing 15% PDS, 2 mM L-glutamine, BME vitamins, and BulletKit SingleQuots (Lonza, Slough, UK). After seven days, the inserts were ready for the experiments.
4.9. Permeability Assay
For the permeability experiments, only inserts with TEER values higher than 300 Ω·cm2 were used. Transwell inserts were washed with Ringer-HEPES buffer (150 mM NaCl, 5.2 mM KCl, 2.2 mM CaCl2, 0.2 mM MgCl2·6H2O, 6 mM NaHCO3, 5 mM HEPES, 2.8 mM glucose; pH 7.4). KYNA-1 was dissolved in Ringer-HEPES and applied to the upper (apical) compartment at a final concentration of 1 mg mL−1. All incubations were performed at 37 °C. Samples were collected from the basolateral side after 5, 15, 30, and 60 min. The concentrations of GA in the abluminal and luminal compartments were measured using a validated LC-MS/MS method. The permeability values of the inserts (PSf, for inserts with a coating only) and the insert plus endothelium (PSt, for inserts with a coating and cells) were considered by applying the following equation:1/PSe = 1/PSt − 1/PSf. To obtain the endothelial permeability coefficient (Pe, in cm·s−1), the permeability value (PSe) corresponding to the endothelium alone was divided by the porous membrane surface area of the insert.
4.10. Statistical Analysis
Statistical analyses were performed using GraphPad Prism 8.0.1 (GraphPad Software, Inc., San Diego, CA, USA). All data were subjected to statistical analysis based on their normal distribution, utilizing the appropriate statistical tests as chosen. Results are presented as mean values with accompanying standard deviations (SD) to indicate the variation.
For quantifying metabolite concentrations and conducting quality assessments, the MetIQ software (
https://biocrates.com/webidq/ accessed on 16 June 2024) package was employed (BIOCRATES Life Sciences AG, Innsbruck, Austria). The internal standards served as a reference for metabolite concentration calculations. Univariate (
t-tests and ANOVA) and multivariate statistics (partial least squares-discrimination analysis PLS-DA), as well as variable importance in projection (VIP) plots, were performed using MetaboAnalyst 5.0 (Xia Lab @ McGill, St. Anne de Bellevue, QC, Canada). Cross-validation of the PLS-DA classification applied five components for selection as the optimal number of components. The LOOCV cross-validation method was used. The performance of measurements, such as Accuracy, R2, and Q2, were considered. The VIP score was measured using the PLS-DA method. The VIP score is a measure of the importance of a feature in the PLS-DA model. This summarizes the contribution of a feature to the model. The VIP score of a feature was calculated as the weighted sum of the PLS weights. The PLS weight is the squared correlation between PLS-DA components and the original feature.