Insect Protein Content Analysis in Handcrafted Fitness Bars by NIR Spectroscopy. Gaussian Process Regression and Data Fusion for Performance Enhancement of Miniaturized Cost-Effective Consumer-Grade Sensors

Future food supply will become increasingly dependent on edible material extracted from insects. The growing popularity of artisanal food products enhanced by insect proteins creates particular needs for establishing effective methods for quality control. This study focuses on developing rapid and efficient on-site quantitative analysis of protein content in handcrafted insect bars by miniaturized near-infrared (NIR) spectrometers. Benchtop (Büchi NIRFlex N-500) and three miniaturized (MicroNIR 1700 ES, Tellspec Enterprise Sensor and SCiO Sensor) in hyphenation to partial least squares regression (PLSR) and Gaussian process regression (GPR) calibration methods and data fusion concept were evaluated via test-set validation in performance of protein content analysis. These NIR spectrometers markedly differ by technical principles, operational characteristics and cost-effectiveness. In the non-destructive analysis of intact bars, the root mean square error of cross prediction (RMSEP) values were 0.611% (benchtop) and 0.545–0.659% (miniaturized) with PLSR, and 0.506% (benchtop) and 0.482–0.580% (miniaturized) with GPR calibration, while the analyzed total protein content was 19.3–23.0%. For milled samples, with PLSR the RMSEP values improved to 0.210% for benchtop spectrometer but remained in the inferior range of 0.525–0.571% for the miniaturized ones. GPR calibration improved the predictive performance of the miniaturized spectrometers, with RMSEP values of 0.230% (MicroNIR 1700 ES), 0.326% (Tellspec) and 0.338% (SCiO). Furthermore, Tellspec and SCiO sensors are consumer-oriented devices, and their combined use for enhanced performance remains a viable economical choice. With GPR calibration and test-set validation performed for fused (Tellspec + SCiO) data, the RMSEP values were improved to 0.517% (in the analysis of intact samples) and 0.295% (for milled samples).     

Determination of Ni2+ using an equilibrium ion exchange technique: Important chemical factors and applicability to environmental samples

The main purpose of this study was to investigate the relationship between some coffee roasting variables (weight loss, density and moisture) with near infrared (NIR) spectra of original green (i.e. raw) and differently roasted coffee samples, in order to test the availability of non-destructive NIR technique to predict coffee roasting degree. Separate calibration and validation models, based on partial least square (PLS) regression, correlating NIR spectral data of 168 representatives and suitable green and roasted coffee samples with each roasting variable, were developed. Using PLS regression, a prediction of the three modelled roasting responses was performed. High accuracy results were obtained, whose root mean square errors of the residuals in prediction (RMSEP) ranged from 0.02 to 1.23%. Obtained data allowed to construct robust and reliable models for the prediction of roasting variables of unknown roasted coffee samples, considering that measured vs. predicted values showed high correlation coefficients (r from 0.92 to 0.98). Results provided by calibration models proposed were comparable in terms of accuracy to the conventional analyses, revealing a promising feasibility of NIR methodology for on-line or routine applications to predict and/or control coffee roasting degree via NIR spectra.     

Improving Prediction of Peroxide Value of Edible Oils Using Regularized Regression Models

We present four unique prediction techniques, combined with multiple data pre-processing methods, utilizing a wide range of both oil types and oil peroxide values (PV) as well as incorporating natural aging for peroxide creation. Samples were PV assayed using a standard starch titration method, AOCS Method Cd 8-53, and used as a verified reference method for PV determination. Near-infrared (NIR) spectra were collected from each sample in two unique optical pathlengths (OPLs), 2 and 24 mm, then fused into a third distinct set. All three sets were used in partial least squares (PLS) regression, ridge regression, LASSO regression, and elastic net regression model calculation. While no individual regression model was established as the best, global models for each regression type and pre-processing method show good agreement between all regression types when performed in their optimal scenarios. Furthermore, small spectral window size boxcar averaging shows prediction accuracy improvements for edible oil PVs. Best-performing models for each regression type are: PLS regression, 25 point boxcar window fused OPL spectral information RMSEP = 2.50; ridge regression, 5 point boxcar window, 24 mm OPL, RMSEP = 2.20; LASSO raw spectral information, 24 mm OPL, RMSEP = 1.80; and elastic net, 10 point boxcar window, 24 mm OPL, RMSEP = 1.91. The results show promising advancements in the development of a full global model for PV determination of edible oils.     

低场核磁共振结合化学计量学方法快速检测掺假核桃油

以掺假核桃油样品为低场核磁共振检测对象,利用主成分分析法(PCA)和偏最小二乘回归法(PLSR)分析处理Carr-Purcell-Meiboom-Gill(CPMG)序列的核磁共振弛豫数据,旨在探求一种能快速检测核桃油品质的新方法。对几种常见掺假形式(掺入大豆油、玉米油、葵花油)的核桃油样品和纯核桃油样品进行检测和评价。实验结果表明:纯核桃油和掺入不同种类食用油的掺假核桃油在主成分得分图上可以得到很好的区分,且掺假样品随掺假比例在图中呈规律性分布;采用PLSR法对CPMG数据和实际掺假率进行回归,可实现对核桃油掺假水平的准确定量测定。方法快速、无损、准确,在食用油制品的品质控制及评价方面具有很大的应用潜力。     

Quantitative Detection of Chromium Pollution in Biochar Based on Matrix Effect Classification Regression Model

Returning biochar to farmland has become one of the nationally promoted technologies for soil remediation and improvement in China. Rapid detection of heavy metals in biochar derived from varied materials can provide a guarantee for contaminated soil, avoiding secondary pollution. This work aims first to apply laser-induced breakdown spectroscopy (LIBS) for the quantitative detection of Cr in biochar. Learning from the principles of traditional matrix effect correction methods, calibration samples were divided into 1–3 classifications by an unsupervised hierarchical clustering method based on the main elemental LIBS data in biochar. The prediction samples were then divided into diverse classifications of calibration samples by a supervised K-nearest neighbor (KNN) algorithm. By comparing the effects of multiple partial least squares regression (PLSR) models, the results show that larger numbered classifications have a lower averaged relative standard deviations of cross-validation (ARSDCV) value, signifying a better calibration performance. Therefore, the 3 classification regression model was employed in this study, which had a better prediction performance with a lower averaged relative standard deviations of prediction (ARSDP) value of 8.13%, in comparison with our previous research and related literature results. The LIBS technology combined with matrix effect classification regression model can weaken the influence of the complex matrix effect of biochar and achieve accurate quantification of contaminated metal Cr in biochar.     

A microwave-based extraction method for the determination of sugar and polyols: Application to the characterization of regular and peaberry coffees

The brewing properties of coffee products are defined by the chemical composition in the bean, including sugars and polyols. Some factors, such as coffee species and roasting, may affect the level of these compounds in the bean. A new analytical microwave-assisted extraction (MAE) method has been developed to extract sugars and polyols from the coffee bean. The studied extraction conditions for the MAE were temperature (30–80 °C), solvent composition (0–50% ethanol in water), and solvent-to-sample ratio (10:1–30:1 mL solvent per g sample). A Box-Behnken design was applied to study the effect of extraction variables, and subsequently, the influential variables were optimized by response surface methodology (RSM). In addition to the main effect of the solvent-to-sample ratio, all quadratic effects significantly influenced (p < 0.05) the recovery of sugars and polyols from the coffee beans. RSM suggested the optimized MAE conditions: temperature 52 °C, ethanol concentration in water 18.5%, and solvent-to-sample ratio 17:1. Under the optimum condition, a kinetics study confirmed that 15 min showed high precision and accuracy of the developed method. Ultimately, a real sample application of the developed MAE revealed that the new method successfully described the composition of sugars and polyols in regular and peaberry coffee beans. Additionally, the method also effectively characterized the green and roasted Arabica and Robusta coffee beans.     

基于近红外光谱技术的电子烟油烟碱含量快速检测研究

烟碱是电子烟烟油中的主要成分,其含量决定了电子烟油的风味口感及产品的安全性。为了提高电子烟油烟碱含量的测量效率,该文采用近红外光谱技术和极限学习机回归(ELMR)建立了电子烟油烟碱含量的定量预测模型。实验结果表明：相比于传统的主成分回归(PCR)和偏最小二乘回归(PLSR)模型,所建立的ELMR预测模型的决定系数R2为0.926 2,远高于PCR预测模型的0.859 0和PLSR预测模型的0.860 4;同时,使用ELMR模型的预测均方根误差(RMSEP)为0.026 8,小于PCR预测模型的0.043 1和PLSR预测模型的0.040 9。以上结果说明该文所建立的近红外光谱定量模型能够应用于烟碱含量的快速准确测量,为实现电子烟油烟碱含量的实时在线监测和其它质量参数的快速测量奠定了良好的基础。     

Identification of gasoline adulteration using comprehensive two-dimensional gas chromatography combined to multivariate data processing

A method to detect potential adulteration of commercial gasoline (Type C gasoline, available in Brazil and containing 25% (v/v) ethanol) is presented here. Comprehensive two-dimensional gas chromatography with flame ionization detection (GCxGC-FID) data and multivariate calibration (multi-way partial least squares regression, N-PLS) were combined to obtain regression models correlating the concentration of gasoline on samples from chromatographic data. Blends of gasoline and white spirit, kerosene and paint thinner (adopted as model adulterants) were used for calibration; the regression models were evaluated using samples of Type C gasoline spiked with these solvents, as well as with ethanol. The method was also checked with real samples collected from gas stations and analyzed using the official method. The root mean square error of prediction (RMSEP) for gasoline concentrations on test samples calculated using the regression model ranged from 3.3% (v/v) to 8.2% (v/v), depending on the composition of the blends; in addition, the results for the real samples agree with the official method. These observations suggest that GCxGC-FID and N-PLS can be an alternative for routine monitoring of fuel adulteration, as well as to solve several other similar analytical problems where mixtures should be detected and quantified as single species in complex samples.     

组合偏最小二乘回归方法在近红外光谱定量分析中的应用

针对近红外光谱数据局部效应显著,变量个数多,彼此间常存在严重的复共线性,并多与样品组分含量呈非线性关系,构建一种组合非线性偏最小二乘回归(E-S-QPLSR)方法。它采用无重复采样技术(subag-ging),从训练样本中生成若干子样,然后每个子样通过二次多项式偏最小二乘回归(QPLSR),建立其子模型,并实现对训练样本因变量的定量预测,再将它们交由线性PLS算法用于计算各子模型的组合权系数。将该法应用于80个玉米样品的水组分含量与其近红外光谱的定量关系建模,效果良好,显示出很强的学习能力,所建模型的预报性能也优于其它方法。     

Ochratoxin A and 2′R-Ochratoxin A in Selected Foodstuffs and Dietary Risk Assessment

The aim of this study was to estimate the contamination of grain coffee, roasted coffee, instant coffee, and cocoa purchased in local markets with ochratoxin A (OTA) and its isomerization product 2′R-ochratoxin A (2′R-OTA), and to assess risk of dietary exposure to the mycotoxins. OTA and 2′R-OTA content was determined using the HPLC chromatography with immunoaffinity columns dedicated to OTA. OTA levels found in all the tested samples were below the maximum limits specified in the European Commission Regulation EC 1881/2006. Average OTA concentrations calculated for positive samples of grain coffee/roasted coffee/instant coffee/cocoa were 0.94/0.79/3.00/0.95 µg/kg, with the concentration ranges: 0.57–1.97/0.44–2.29/0.40–5.15/0.48–1.97 µg/kg, respectively. Average 2′R-OTA concentrations calculated for positive samples of roasted coffee/instant coffee were 0.90/1.48 µg/kg, with concentration ranges: 0.40–1.26/1.00–2.12 µg/kg, respectively. In turn, diastereomer was not found in any of the tested cocoa samples. Daily intake of both mycotoxins with coffee/cocoa would be below the TDI value even if the consumed coffee/cocoa were contaminated with OTA/2′R-OTA at the highest levels found in this study. Up to now only a few papers on both OTA and 2′R-OTA in roasted food products are available in the literature, and this is the first study in Poland.     

Prediction of sensory properties of espresso from roasted coffee samples by near-infrared spectroscopy

Thirty-five representative and suitably selected roasted coffee samples were characterised by near-infrared (NIR) spectroscopy and used to prepare the corresponding espresso samples to be subsequently subjected to sensory evaluation by trained panellists. The main purpose was to investigate the relationships between certain crucial sensory attributes of espresso coffees, including perceived acidity, mouthfeel, bitterness and aftertaste, and near-infrared spectra of original roasted coffee samples, in such a way that non-destructive near-infrared reflectance measurements would be used to predict all these sensory properties with a decisive influence from a quality assurance standpoint. Separate calibration models based on partial least squares regression (PLS), correlating NIR spectral data of roasted coffee samples with each sensory attribute of espresso samples studied, were developed. Wavelength selection was also performed applying iterative predictor weighting-PLS (IPW-PLS) in order to take into account only significant and characteristic spectral features, in an attempt to improve the quality of the final regression models constructed. Using IPW-PLS regression, prediction of the four sensory responses modelled was performed with high accuracy, with root mean square errors of the residuals in cross-validation (RMSECV) ranging from 4.7 to 7.0%. Thus, the results provided by the high-quality calibration models proposed in the present study, comparable in terms of accuracy to the evaluations provided by a trained sensory panel, are promising and prove the feasibility of using a similar methodology in on-line or routine applications to predict the sensory quality of unknown espresso coffee samples via their respective NIR roasted coffee spectra.     

High and low frequency unfolded partial least squares regression based on empirical mode decomposition for quantitative analysis of fuel oil samples

Accurate prediction of the model is fundamental to the successful analysis of complex samples. To utilize abundant information embedded over frequency and time domains, a novel regression model is presented for quantitative analysis of hydrocarbon contents in the fuel oil samples. The proposed method named as high and low frequency unfolded PLSR (HLUPLSR), which integrates empirical mode decomposition (EMD) and unfolded strategy with partial least squares regression (PLSR). In the proposed method, the original signals are firstly decomposed into a finite number of intrinsic mode functions (IMFs) and a residue by EMD. Secondly, the former high frequency IMFs are summed as a high frequency matrix and the latter IMFs and residue are summed as a low frequency matrix. Finally, the two matrices are unfolded to an extended matrix in variable dimension, and then the PLSR model is built between the extended matrix and the target values. Coupled with Ultraviolet (UV) spectroscopy, HLUPLSR has been applied to determine hydrocarbon contents of light gas oil and diesel fuels samples. Comparing with single PLSR and other signal processing techniques, the proposed method shows superiority in prediction ability and better model interpretation. Therefore, HLUPLSR method provides a promising tool for quantitative analysis of complex samples.     

Preparation,structure elucidation,and antioxidant activity of new bis(thiosemicarbazone) derivatives

Schiff-base–bearing new bis(thiosemicarbazone) derivatives were prepared from terephthalaldehyde and various thiosemicarbazides. FT–IR, ¹H NMR, ¹³C NMR, and UV–Vis spectroscopic methods and elemental analysis were used to elucidate the identification of the synthesized molecules. The in vitro antioxidant activity of the synthesized compounds was analysed with the 1,1-diphenyl-2-picryl hydrazyl free-radical–trapping process. The synthesized compounds exhibited lower antioxidant activity than the standard ascorbic acid. IC₅₀ values of the synthesized molecules measured from 3.81 ± 0.01 to 29.05 ± 0.11 μM. Among the synthesized compounds, compound 3 had the best antioxidant activity. Moreover, this study explained the structure–activity relationship of the synthesized molecules with different substituents in radical trapping reactions.     

Sequential (step-by-step) detection,identification and quantitation of extra virgin olive oil adulteration by chemometric treatment of chromatographic profiles

An analytical method for the sequential detection, identification and quantitation of extra virgin olive oil adulteration with four edible vegetable oils--sunflower, corn, peanut and coconut oils--is proposed. The only data required for this method are the results obtained from an analysis of the lipid fraction by gas chromatography-mass spectrometry. A total number of 566 samples (pure oils and samples of adulterated olive oil) were used to develop the chemometric models, which were designed to accomplish, step-by-step, the three aims of the method: to detect whether an olive oil sample is adulterated, to identify the type of adulterant used in the fraud, and to determine how much aldulterant is in the sample. Qualitative analysis was carried out via two chemometric approaches--soft independent modelling of class analogy (SIMCA) and K nearest neighbours (KNN)--both approaches exhibited prediction abilities that were always higher than 91% for adulterant detection and 88% for type of adulterant identification. Quantitative analysis was based on partial least squares regression (PLSR), which yielded R2 values of >0.90 for calibration and validation sets and thus made it possible to determine adulteration with excellent precision according to the Shenk criteria.     

Non-Destructive Analysis of Chlorpheniramine Maleate Tablets and Granules by Chemometrics-Assisted Attenuated Total Reflectance Infrared Spectroscopy

Non-destructive analysis of chlorpheniramine maleate (CPM), pharmaceutical tablets, and granules was conducted by chemometrics-assisted attenuated total reflectance infrared spectroscopy (ATR-IR). For tablets, an optimum PLSR model with eight latent factors was obtained from area-normalized and standard normal variate (SNV) pretreated ATR-IR spectral data with correlation coefficients (R²) of calibration and cross-validation of 0.9716 and 0.9602, respectively. The model capability for the 42 test set samples was proven with R² between the reference and model prediction values of 0.9632, and a root-mean-square error of prediction (RMSEP) of 1.7786. The successive PLSR model for granules was constructed from SNV and first derivative pretreated ATR-IR spectral data with two latent factors and correlation coefficients (R²) of calibration and cross-validation of 0.9577 and 0.9450, respectively.     

Quantitative analysis of two adulterants in Cynanchum stauntonii by near-infrared spectroscopy combined with multi-variate calibrations

Authentication of traditional Chinese medicines (TCMs) has become important because they can be adulterated with relatively cheap herbal medicines similar in appearance. Detection of such adulterated samples is needed because their presence is likely to reduce the pharmacological potency of the original TCM and, in the worst cases, the samples may be harmful. The aim of this study was to develop a rapid near-infrared spectroscopy (NIRS) analytical method which was supported by multi-variate calibration, e.g. partial least squares regression (PLSR) and radial basis function artificial neural networks (RBF-ANN), in order to quantify the TCM and the adulterants. In this work, Cynanchum stauntonii (CS), a commonly used TCM, in mixtures with one or two adulterants ?? two morphological types of TCM, Cynanchum atrati (CA) and Cynanchum paniculati (CP), were determined using NIR reflectance spectroscopy. The three sample sets, CS adulterated with CA or CP, and CS with both CA and CP, were measured in the range of 800?C2500 nm. Both PLSR and RBF-ANN calibration models provided satisfactory results, even at an adulteration level of 5 mass %, but the RBF-ANN models with better root mean square error of prediction (RMSEP) values for CS, CA, and CP arguably performed better. Consequently, this work demonstrates that the NIR method of sampling complex mixtures of similar substances such as CS adulterated by CA and/or CP is capable of producing data suitable for the quantitative analysis of mixtures consisting of the original TCM adulterated by one or two similar substances, provided the spectral data are interrogated by multi-variate methods of data analysis such as PLS or RBF-ANN.     

Application of MOS Gas Sensors Coupled with Chemometrics Methods to Predict the Amount of Sugar and Carbohydrates in Potatoes

Five potato varieties were studied using an electronic nose with nine MOS sensors. Parameters measured included carbohydrate content, sugar level, and the toughness of the potatoes. Routine tests were carried out while the signals for each potato were measured, simultaneously, using an electronic nose. The signals obtained indicated the concentration of various chemical components. In addition to support vector machines (SVMs that were used for the classification of the samples, chemometric methods, such as the partial least squares regression (PLSR) method, the principal component regression (PCR) method, and the multiple linear regression (MLR) method, were used to create separate regression models for sugar and carbohydrates. The predictive power of the regression models was characterized by a coefficient of determination (R²), a root-mean-square error of prediction (RMSEP), and offsets. PLSR was able to accurately model the relationship between the smells of different types of potatoes, sugar, and carbohydrates. The highest and lowest accuracy of models for predicting sugar and carbohydrates was related to Marfona potatoes and Sprite cultivar potatoes. In general, in all cultivars, the accuracy in predicting the amount of carbohydrates was somewhat better than the accuracy in predicting the amount of sugar. Moreover, the linear function had 100% accuracy for training and validation in the C-SVM method for classification of five potato groups. The electronic nose could be used as a fast and non-destructive method for detecting different potato varieties. Researchers in the food industry will find this method extremely useful in selecting the desired product and samples.     

Comparison of calibrations for the determination of soluble solids content and pH of rice vinegars using visible and short-wave near infrared spectroscopy

Visible and short-wave near infrared (Vis/SWNIR) spectroscopy combined with chemometrics was investigated for the fast determination of soluble solids content (SSC) and pH values of rice vinegars. Two hundred and twenty-five samples (45 for each variety) were selected randomly for the calibration set, whereas, 75 samples (15 for each variety) for the validation set, and the remaining 25 samples for the independent set. After some preprocessing, partial least squares (PLS) analysis was implemented for calibration models with different wavelength bands including visible, SWNIR and Vis/SWNIR regions. The best PLS models were achieved with Vis/SWNIR (550–1000 nm) region. Furthermore, different latent variables (5, 6, 7, 8 LVs) were used as inputs of least squares-support vector machine (LS-SVM) to develop the LV-LS-SVM models with grid search technique and RBF kernel. The optimal models were obtained with 6 LVs and they outperformed PLS models. Moreover, effective wavelengths (EWs) were selected according to regression coefficients. The EW-LS-SVM models were developed and an excellent prediction precision was achieved, and the effectiveness of the EWs was also validated. The correlation coefficient (r), root mean square error of prediction (RMSEP) and bias for the best prediction by EW-LS-SVM were 0.999, 0.189 and 0.051 for SSC, whereas 0.999, 0.008 and −1.7 × 10⁻³ for pH, respectively. The overall results indicated that Vis/SWNIR spectroscopy combined with LS-SVM could be applied as a high precision and fast way for the determination of SSC and pH values of rice vinegars.     

茶叶浸出液中农药残留表面拉曼光谱检测与初步定量分析

选取甲基对硫磷和水胺硫磷为研究对象,改良了传统的QuEChERS前处理工艺,以自制纳米金溶胶为增强基底,利用表面增强拉曼光谱(SERS)技术,对茶叶浸出液中的农药残留进行检测。通过比对两种有机磷农药的拉曼特征峰进行定性分析。同时,选取570,1034,1107和1202 cm^-1等拉曼位移附近的特征峰光谱数据,利用微分等数学手段,结合偏最小二乘法(PLSR)建立回归方程,预测样品中农药残留含量。所得预测数值与气相色谱-质谱联用(GC-MS)法检测值对比,验证本方法的可行性与可信度。结果表明:基于SERS技术对上述两种有机磷农药的检出限可达0.05 mg/L;通过数学模型分析建立回归方程,其线性相关系数范围为0.9077~0.9824,预测均方根误差(RMSEP)范围为0.77%~2.68%;利用回归方程得到的预测值与GC-MS检测结果基本接近,相对误差范围-5.16%~9.03%,回收率为81.4%~115.1%,说明可以用SERS技术对茶叶浸出液中的有机磷农药残留进行定性和初步定量分析。     

On the Identification and Quantification of Ergothioneine and Lovastatin in Various Mushroom Species: Assets and Challenges of Different Analytical Approaches

In recent years, mushrooms have drawn the attention of agro-industries and food-industries as they were considered to be valuable natural sources of health promoting compounds such as β-glucans, ergothioneine, and lovastatin. The detection and quantification of such compounds by implementing reliable analytical approaches is of the utmost importance in order to adjust mushrooms’ cultivation conditions and maximize the production in different species. Toward this direction, the current study focuses on the comparison of ultraviolet–visible (UV–Vis) spectrometry and liquid chromatography–mass spectrometry (LC–MS) methods (a) by evaluating the content of ergothioneine and lovastatin in mushrooms and (b) by highlighting any possible substrate-based interferences that hinder the accurate determination of these two compounds in order to propose the technique-of-choice for a standardized bioactive compounds monitoring. For this purpose, mushrooms produced by three species (i.e., Agaricus bisporus, Pleurotus ostreatus, and P. citrinopileatus) on various cultivation substrates, namely wheat straw (WS), winery (grape marc (GM)), and olive oil (OL) by-products, were examined. Among the two applied techniques, the developed and validated LC–MS methods, exhibiting relatively short analysis time and higher resolution, emerge as the methods-of-choice for detecting ergothioneine and lovastatin in mushrooms. On the contrary, UV–Vis methods were hindered due to co-absorbance of different constituents, resulting in invalid results. Among the studied mushrooms, P. citrinopileatus contained the highest amount of ergothioneine (822.1 ± 20.6 mg kg⁻¹ dry sample), whereas A. bisporus contained the highest amounts of lovastatin (1.39 ± 0.014 mg kg⁻¹ dry sample). Regarding the effect of different cultivation substrates, mushrooms produced on OL and WS contained the highest amount of ergothioneine, while mushrooms deriving from GM-based substrates contained the highest amount of lovastatin.