Research Article | Open Access
Quality Control of Mutton by Using Volatile Compound Fingerprinting Techniques and Chemometric Methods
A method for chromatographic fingerprinting of flavor was established for the quality control of mutton. Twenty-five mutton samples that were chosen from twelve batches were investigated by gas chromatography-mass spectroscopy (GC-MS) and gas chromatography-olfactometry (GC-O). Spectral correlative chromatograms combined with GC-O assessment were employed, and 32 common odor-active compounds that characterize mutton flavor fingerprint were obtained. Based on the flavor chromatographic fingerprint data, principal component analysis (PCA) and partial least squares-discriminant analysis (PLS-DA) were designed and employed as chromatographic fingerprint methods. Defined categories were perfectly discriminated after PLS-DA was conducted on the fused matrix, demonstrating a 100% accurate classification. Fourteen constituents were further screened with PLS-DA to be the main chemical markers, and they were used to develop similar approaches for the determination of mutton quality and traceability. The flavor fingerprint of mutton established using SPME-GC-MS/O coupled with PLS-DA is appropriate for differentiating and identifying samples, and the procedure would be used in quality control.
Mutton is very popular in Western countries and is widely consumed because of its unique flavor and its nutritional qualities . Meat flavor is mainly produced by the interactions of water-soluble compounds and lipids . Water-soluble compounds, which are present in lean meat, generate the basic meat flavor. Lipids (meat fat) impart the characteristic flavor of distinct meat species . Studies have shown that if fat was removed from meat, meat species are difficult to distinguish based only on their sensory attributes . As such, mutton is often adulterated by low-priced meats (e.g., duck or pork) combined with mutton fat to maximize profit. A continuous supply of adulterated mutton products is mixed with real mutton, resulting in the gradual deterioration of consumer confidence on mutton industry. Therefore, reliable and accurate methods should be established for the quality control of mutton.
Mutton is a complex system with multiple components . Thus, its authenticity and quality cannot be effectively evaluated by simply measuring one or few compounds or indices. Prior to this study, no report has been published on quality control of mutton. Aroma is one of the most important aspects and composition indicators of meat flavor . The sensory attributes of adulterated mutton samples are similar to those of real mutton, but their aromatic volatiles are possibly distinct in terms of category, quantity, and composition. Thus, volatile compound evaluation is prioritized in specific flavor studies, and the unique composition and content of various aromatic compounds in mutton can be used to identify mutton samples and serve as the standard in quality evaluation. Therefore, unlike conventional methods that focus mainly on determining certain whole compounds, flavor fingerprinting can characterize a complex system with quantitative reliability .
Chromatographic fingerprinting technique is a unique method used in distinguishing authentic from counterfeit products and in evaluating the quality of fruits, vegetables , tea , beverages , flavors, and fragrances , and so on. This technique was also used to evaluate the quality of meat products. For instance, Tian et al.  confirmed that the moderate oxidized suet may be a desirable precursor for mutton process flavors by using a PLSR model. Considered effective for accessing the quality of raw materials and/or their processing products, this technique is a system model that is established on the basis of data acquired using chromatographic techniques, thin-layer chromatography, X-ray, gas chromatography (GC), high-performance liquid chromatography, and so on to establish the specific patterns of recognition for multiple compounds in the samples . This technique is basically characterized by integrity and ambiguity . Ambiguity emphasizes fingerprint similarity rather than identity between control and test samples, whereas integrity emphasizes the integral composition and features of sample components, that is, the integral chromatographic features. However, minor differences between the mutton and fraud meat products may not be observed because of the complexity of the analyzed samples. Nonetheless, pattern recognition methods, such as linear discriminant analysis (LDA), principal component analysis (PCA), partial least squares-discriminant analysis (PLS-DA), and soft-independent modeling of class analogy, facilitate the visualization of the information in the fingerprints. In this research, a method was designed in which flavor fingerprints were used in PLS-DA and PCA to build a target pattern recognition system for the quality control of mutton.
This research was primarily aimed at determining the characteristic flavor fingerprint of mutton through GC-MS and GC-O. Moreover, PCA and PLS-DA were used to distinguish mutton from other kinds of meat or adulterated products for the quality control of mutton. Considering the above-mentioned analysis, the new GC-MS/GC-O analysis described here coupled with chemometric methods is a promising methodology to ensure mutton quality and achieve comprehensive quality control for mutton.
2. Materials and Methods
2.1. Materials and Chemicals
Different batches of mutton samples obtained from different sources, as well as horsemeat samples, were provided by Beitun EHe steppe Food Co., Ltd. (Xinjiang, China). The sheep were reared indoors and were fed a maize silage-based diet. They had free access to hay and straw; they were slaughtered at 280 d of age. Pork, chicken, duck, and beef were purchased from the Wuxi Wal-Mart store (Wuxi, China). C6–C26, a mixture of serial -alkanes, was obtained from Sinopharm Chemical Reagent Co., Ltd. (Shanghai, China). Reference components were purchased from Sigma Aldrich Co., Ltd. (Shanghai, China).
2.2. Sample Preparation
Bones were immediately removed from the mutton samples, which were divided into five parts for analysis (neck, forelegs, sternal rib, abdomen, and draw legs with buttocks). Each part weighed 500 g and was wrapped in aluminum foil. Parts were then stored at −20°C for further analysis. Prior to analysis, the meat samples were defrosted for 24 h at 4°C and cut into 0.5 cm × 0.5 cm × 0.5 cm cubes. Each part of the meat samples was accurately weighed (100.0 g; total weight, 500.0 g), mixed with 500 mL of tap water, and pressure-cooked (121 Pa) for 30 min. The samples were then immediately cooled in ice water, sealed, and then stored at −18°C for subsequent analysis.
2.3. GC-MS Analysis
The volatiles were sampled by using a 75 μm CAR/PDMS SPME fiber (Supelco, USA). Three grams of the meat samples was distributed evenly in a 15 mL vial and then mixed by using a magnetic stirring bar. Prior to trapping, 1,2-dichlorobenzene (4 μL, 0.555 μg/μL) used as an internal standard was added to each sample. The vial was closed and sealed by a PTFE/BYTL septum, placed in a 55°C water bath. Then the fiber was exposed in the sample headspace, maintained for 30 min, and then transferred immediately into the GC injection port. Then the fiber was desorbed for 5 min at 250°C.
Chromatographic separation was performed using a Finnigan trace GC-MS (Finnigan, USA) equipped with a FID detector. The volatile compounds were separated on a DB-WAX column (J & W Scientific Inc., USA; 30 m length, 0.25 mm i.d., 0.25 μm ). The carrier gas was helium (99.999%) and had a flow rate of 1.8 mL/min. The oven temperature was programmed to run for 3 min at 40°C, subsequently increased to 80°C at 6°C/min, then increased to 230°C at 10°C/min, and finally held for 10 min. The detection temperature was 250°C. MS spectra were acquired under electron impact mode at an emission current of 35 mA and ionization voltage of 70 eV. The transfer line temperature was 250°C, whereas the quadrupole mass filter was employed at 150°C. A mass range of to 450 was noted at the rate of 4.45 scan/s.
Comparing the Kovats retention indices (KI) and the mass spectra of individual components with those of authentic standard compounds allowed identification of the volatile components. They were tentatively identified using WILEY and NIST and KI. The relative concentration of volatile components was calculated on the basis of their peak area with those of 1,2-dichlorobenzene. The expression of KI value is defined as follows:where is the unknown compound and and are the number of carbon atoms of the n-alkanes eluted before and after the unknown compound , respectively. , , and are the retention times of the unknown compound , n-alkanes of carbon number , and , respectively.
2.4. GC-O Analysis
GC-O was conducted with Finnigan trace gas chromatograph equipped with an OP275 sniffing port (GL Sciences Inc., Japan) and an FID. The effluent from the DB-WAX column was divided at 1 : 1 (v/v) between the sniffing port and FID. Nitrogen was the carrier gas and it had a flow velocity of 2.1 mL/min. The samples were collected, and the column and analyzing conditions were all similar to those in GC-MS. A detection frequency method involving 8 well-trained panelists (4 males and 4 females, aged range 26–48) was used to acquire the odor profile of mutton samples. The method was used as previously described .
2.5. Statistical Data Processing
ANOVA was carried out using SPSS18.0. PCA and PLS-DA were conducted using Unscrambler X version 9.6. The chemometric methods were performed using the average of 3 gas chromatographic peaks that were acquired from triplicates of each sample.
Repeatability, precision, and stability were evaluated by calculating their relative standard deviations (RSDs) to validate analytical method for the flavor fingerprinting. Repeatability was assessed by analyzing 5 different preparations of the same sample. Precision was evaluated by reduplicating continuous injections of the same sample for 5 times. Regarding sample stability, it was performed by analyzing the same treated sample that was stored for 36 h at RT.
The common volatile compounds in the fingerprints acquired through GC-MS were analyzed using PCA to obtain a visual representation of the similarity among samples. First, thirty-two characteristic compounds of mutton flavor fingerprint were extracted from the mutton samples. The data for the compound were recorded as 0 if this compound was not detected in one sample. Then, the dataset for the analytical samples was performed to be the input data through mean centering of the 32 common characteristic components in place of the full spectrum of the chromatogram. Then, this dataset was integrated into the map matrix for further analysis. In this part of the method, PCA was carried out to distinguish and identify different meat or adulterated samples. The samples included 25 mutton samples, 5 beef samples, five chicken samples, 5 pork samples, 5 horse meat samples, 5 duck samples, and 5 fraud mutton samples. Each sample was analyzed in triplicate under similar conditions. The sample data matrices were denoted as M1–25 (mutton), B1–5 (beef), C1–5 (chicken), P1–5 (pork), H1–5 (horse meat), D1–5 (duck), and F1–5 (fraud samples).
PLS-DA was conducted to predict the identity of the mutton samples by using the common volatile data. This technique combines PLSR and DA. Thus, this process demonstrates the advantages of PLS, that is, feature extraction and noise reduction, and is applicable particularly in situations involving few samples but numerous variables . The samples analyzed and the methods for data processing in PLS-DA were similar to those in PCA. First, six of the above-mentioned known meat samples and fraud samples were selected for PLS modeling. Based on the fingerprint information, the classified variables of 25 random mutton samples, considered as training samples, were set as 1, and those of the other samples were labeled as 0. The samples in the training set were subsequently subjected to regression analysis to establish the PLS models of the chromatographic data and classified variables. The -axis represents the 32 characteristic compounds, and the -axis represents the classified variables (1 or 0). All variables were scaled and centered to 1/Sedv. The PLS model was then validated using full cross-validation to set the unit variance and a zero mean for each variable before performing PLS analysis. This ensured an unbiased contribution of each variable to a criterion. The other 25 external mutton samples and 5 fraud mutton samples that were not used in constructing the model were predicted on the basis of the threshold value (ultimate result with ±0.5 deviation) calculated by the model. The mutton samples were identified according to a cut-off value of 0.5. In other words, samples with prediction value of >0.5 and deviation of <0.5 were identified as mutton, and samples with prediction value of <0.5 and deviation of <0.5 were identified as other meats or fraud samples. When the deviation exceeded 0.5, the established identification model was regarded as unstable.
3. Results and Discussion
3.1. Standardization of Chromatographic Fingerprint
A chromatographic fingerprint is a mutual pattern of several common compounds with some characteristics in all the analyzed samples. It should feature the fundamental attributes of fuzziness and integrity or differences and sameness. The mutton could be accurately authenticated and identified (integrity) with these fingerprints even if the concentrations of the characteristic chemical components were not exactly similar among samples (fuzziness). Thus, the most important aspect of building a reliable fingerprint and standardizing its information feature is the extraction and quantification of the common characteristic compounds. Twenty-five mutton samples (M1–M25) were analyzed by GC-MS using the same experimental conditions.
Generally, the average chromatograms are considered standardized characteristic fingerprint information, and this was also partly adopted in this study. In this study, 162 peaks were detected but only 73 existed in all mutton samples. These 73 coexistent compounds are considered the “common components” used to establish a fingerprint. For flavor fingerprints, however, the amount of a compound did not necessarily reflect the intensity of the perceived aroma due to the different odor thresholds for distinct compounds. Therefore, contrary to traditional fingerprints, the common characteristic compounds in flavor fingerprint should also be odor-active, which can be tested through GC-O assessment with a detection frequency of >50%. Based on the detection frequency results, only 32 compounds possessed odor activities. Thus, these compounds can be assigned as the common characteristic compounds. Then, the methods used in quantifying the above 32 common compounds to establish a standard fingerprint require further investigation. However, several phenomena, especially for the changes in the response values of the peaks during data processing, may reduce the accuracy of the fingerprints. This problem can usually be resolved by selecting a stable peak as reference peak or adding an internal standard substance. In this study, 1,2-dichlorobenzene was chosen to be an internal standard when calculating the relative concentrations. The average chromatographic information of selected common odor-active components was regarded as the standardized characteristic flavor fingerprint of mutton, according to this method (Table 1).
|Code representing the 32 active-odor compounds observed in GC-MS-O analysis; KI: Kovats index in column DB-Wax; values reported are means ± S.D. (µg/kg); identification method: (A) mass spectrum and KI agree with that of the authentic compound run under similar GC-MS conditions; (B) tentative identification based on interpretation of mass spectrum and comparison with similar compound. (C) Mass spectrum and KI agree with literature data or NIST Chemistry WebBook (http://webbook.nist.gov/chemistry/name-ser/).|
Among the compounds mentioned in Table 1, three have been determined as the decisive characteristic flavor compounds of mutton: 4-ethyloctanoic acid (cheese, boiled mutton), 4-methylnonanoic acid (tallow, roasted, mutton), and hexanoic acid (lamb, oily) . Moreover, some aliphatic aldehydes, especially some short-chain aliphatic aldehydes (C5–C9), were also reported in sheep meat; these aldehydes are as follows: hexanal (fruity, fatty, and green), heptanal (nutty, oily, and fruity), octanal (fruity, floral, and fatty), nonanal (oily, fatty), and 2-octenal (fruity) . Furthermore, 6-methyl-5-hepten-2-one, 3-methylbutanal, 2-methyl-5-(methylthio) furan, 2-undecenal, and 2,3,5,6-tetramethyl pyrazine were related to emitting meaty odors, and all of these components influence the overall mutton flavor to a certain extent. 4-Methyloctanoic acid is described as the odor of mutton , but it was not detected in this study.
The method was then validated regarding its repeatability, precision, and stability. The instrument/injection precision of the method, as represented by RSD, was lower than 1.02% for relative retention time (RRT) and within 0.28% to 2.96% for relative peak area (RPA). The RSD of the repeatability of RRT and RPA were lower than 0.64% and 3.78%, respectively. In addition, sample stability, which was also evaluated by RSD, ranged from 0.21% to 1.32% for RRT and from 0.79% to 5.01% for RPA. These results showed that the method is applicable and reliable for the analysis and establishment of the mutton flavor fingerprint and subsequent analysis.
3.2. Application of the Fingerprint in Quality Control of Mutton
The developed mutton flavor fingerprint contained information that discriminates mutton samples from other meat samples and from adulterated samples as verified by the PCA and PLS-DA loadings.
The derived PCA showed three significant PCs, which explained 61% of the total variance. Thus, PC1-PC2 and PC1-PC3 were mapped (Figure 1). According to the fundamentals of PCA, proximity to score plots is associated with high similarity with respect to the two series of components concerned; otherwise, they differ from each other . As shown in Figure 1, it is observed that 25 mutton samples and 25 other meat samples were located in different positions, which indicated that the samples significantly differed in two-dimensional space projection. For adulterated mutton samples, (pork mixed with artificial mutton flavor), (duck mixed with synthetic simulated mutton flavor), (duck mixed with 20% lean mutton), and (pork mixed with 20% lean mutton) were easily distinguished from mutton samples. However, (lean pork mixed with mutton suet) are nearly similar to mutton samples (Figure 1(b)). This result demonstrated that sample pork mixed with mutton fat can produce a flavor similar to that of normal mutton. This phenomenon may be attributed to the fact that the characteristic flavor of mutton is mainly attributed to mutton suet, and lean meat only provides the basic meat flavor. This result is similar to those of other studies. Wasserman and Talley (1968) found that when beef tallow was added to veal, evaluators could distinguish veal from beef; however, when lamb fat was added to veal, veal is more easily mistaken as lamb. The above results indicated that PCA performed based on fingerprint information is generally useful in the quality control of mutton, but it is limited for adulterated samples, particularly for samples mixed with mutton suet.
Considering the above results, an alternative strategy can be performed, which is using a supervised classification technique, such as PLS-DA. The data matrix in the PLS-DA model was comprised of 55 rows (samples in the training set) and 32 columns (data on common characteristic volatiles). Figure 2 showed the plot correlating predicted versus measured classified variables is a reasonably good predictor of mutton samples in most cases. The root-mean-square error for the prediction was 0.209 (lower than the error standard <0.5), and the explained variance was 0.880 (83.9% by cross-validation). These results revealed that the PLS-DA model achieved a good fitting without overfitting and is thus applicable in mutton sample classification and identification.
The observations in the different groups of the prediction set, namely, mutton, other meats, and fraud mutton samples, were accurately clustered close to their corresponding groups in the PLS-DA model (Figure 3(a)). The model achieved was validated through an external test set of samples to test the stability of PLS-DA. The identification and prediction results are shown in Figure 4, whereas the specific analytical data are shown in Table 2. The predicted scores of the 25 external mutton samples achieved prediction values of >0.640 (latent variable 1) with deviation of <0.22, whereas the adulterated samples, which were analyzed via a method similar to PCA, showed prediction values of 0.189–0.473 with deviation of <0.31. According to the PLS-DA fundamentals, mutton samples can be identified correctly by using the PLS-DA model, indicating that the proposed model can sufficiently discriminate mutton and other samples.
|The labels of the samples correspond to the code in Figure 4.|
Furthermore, PLS-DA loadings were graphed to identify the possible volatiles that contributed to sample discrimination. Loadings (E)-2-nonenal (tallowy, green), (E)-2-decenal (rosin, green), (E,E)-2,4-octadienal (fatty, boiled meat), tetradecanal (buttery, fruity, green), octadecanal (fatty, candle), acetic acid (pungency, sour), 4-ethyloctanoic acid (cheese, boiled mutton), decanoic acid (rancid, oily), nonanoic acid (fatty, cheese), 2-undecanone (oily, fruity), 2-pentadecanone (oily), ethyl dodecanoate (wax, fruity, fatty), 2-pentylfuran (earthy, meaty), and 2-methyl-5-(methylthio)furan (roasted, peanut) were the variables that contributed the most to sample discrimination (Figure 3(b)). These 14 compounds should be used as chemical markers for discriminating the internal quality mutton samples from that of the other meat samples, which would indicate that the quality of mutton samples could be assessed quantitatively by determining the amounts of these chemical markers. Thus, the above-mentioned compounds possibly influence the discrimination of mutton samples.
GC-MS combined with GC-O was established for flavor fingerprint analysis to monitor the quality of mutton. Information on the common flavor fingerprint was acquired from 32 aroma-active volatile compounds present in 25 mutton samples. The results were then validated for the sample stability, repeatability, and reproducibility. The developed flavor fingerprint is appropriate for differentiating and identifying mutton samples by geographical origin when used in combination with the appropriate chemometric method (PLS-DA). This method can thus be used for quality control.
First author is Ping Zhan.
Conflicts of Interest
The authors have declared no conflicts of interest.
This research was conducted with support from National Key Research and Development Program of China (Grant no. 2016YFD0400705), National Natural Science Foundation of China (Grant no. 31460408), and China Postdoctoral Science Foundation Funded Project (Grant no. 2016M591029). The authors would like to thank Professor Liu Jingsheng for their critical help in the design of research content.
- V. Banskalieva, T. Sahlu, and A. L. Goetsch, “Fatty acid composition of goat muscles and fat depots: a review,” Small Ruminant Research, vol. 37, no. 3, pp. 255–268, 2000.
- Y.-C. Huang, H.-J. Li, Z.-F. He, T. Wang, and G. Qin, “Study on the flavor contribution of phospholipids and triglycerides to pork,” Food Science and Biotechnology, vol. 19, no. 5, pp. 1267–1276, 2010.
- J. D. Wood, M. Enser, A. V. Fisher et al., “Fat deposition, fatty acid composition and meat quality: a review,” Meat Science, vol. 78, no. 4, pp. 343–358, 2008.
- F. B. Whitfield and D. S. Mottram, “Volatiles from interactions of maillard reactions and lipids,” Critical Reviews in Food Science and Nutrition, vol. 31, no. 1-2, pp. 1–58, 1992.
- M. M. Sutherland and J. M. Ames, “The effect of castration on the headspace aroma components of cooked lamb,” Journal of the Science of Food and Agriculture, vol. 69, no. 4, pp. 403–413, 1995.
- V. C. Resconi, M. M. Campo, F. Montossi, V. Ferreira, C. Sañudo, and A. Escudero, “Relationship between odour-active compounds and flavour perception in meat from lambs fed different diets,” Meat Science, vol. 85, no. 4, pp. 700–706, 2010.
- S. Song, X. Zhang, K. Hayat, Z. Xiao, Y. Niu, and K. Eric, “Coordinating fingerprint determination of solid-phase microextraction/gas chromatography-mass spectrometry and chemometric methods for quality control of oxidized tallow,” Journal of Chromatography A, vol. 1278, pp. 145–152, 2013.
- A. Z. Berna, S. Buysens, C. Di Natale, I. U. Grün, J. Lammertyn, and B. M. Nicolaï, “Relating sensory analysis with electronic nose and headspace fingerprint MS for tomato aroma profiling,” Postharvest Biology and Technology, vol. 36, no. 2, pp. 143–155, 2005.
- M. Dumarey, I. Smets, and Y. Vander Heyden, “Prediction and interpretation of the antioxidant capacity of green tea from dissimilar chromatographic fingerprints,” Journal of Chromatography B: Analytical Technologies in the Biomedical and Life Sciences, vol. 878, no. 28, pp. 2733–2740, 2010.
- B. J. Savary, P. Vasu, A. Nunez, and R. G. Cameron, “Identification of thermolabile pectin methylesterases from sweet orange fruit by peptide mass fingerprinting,” Journal of Agricultural and Food Chemistry, vol. 58, no. 23, pp. 12462–12468, 2010.
- X. Zhu, Y. Gao, Z. Chen, and Q. Su, “Development of a chromatographic fingerprint of tobacco flavor by use of GC and GC-MS,” Chromatographia, vol. 69, no. 7-8, pp. 735–742, 2009.
- H. Tian, P. Zhan, W. Li et al., “Contribution to the aroma characteristics of mutton process flavor from oxidized suet evaluated by descriptive sensory analysis, gas chromatography, and electronic nose through partial least squares regression,” European Journal of Lipid Science and Technology, vol. 116, no. 11, pp. 1522–1533, 2014.
- T. Cajka, K. Riddellova, M. Tomaniova, and J. Hajslova, “Recognition of beer brand based on multivariate analysis of volatile fingerprint,” Journal of Chromatography A, vol. 1217, no. 25, pp. 4195–4203, 2010.
- B. Bertrand, R. Boulanger, S. Dussert et al., “Climatic factors directly impact the volatile organic compound fingerprint in green Arabica coffee bean as well as coffee beverage quality,” Food Chemistry, vol. 135, no. 4, pp. 2575–2583, 2012.
- P. Zhan, H. Tian, X. Zhang, and L. Wang, “Contribution to aroma characteristics of mutton process flavor from the enzymatic hydrolysate of sheep bone protein assessed by descriptive sensory analysis and gas chromatography olfactometry,” Journal of Chromatography B: Analytical Technologies in the Biomedical and Life Sciences, vol. 921-922, pp. 1–8, 2013.
- S. Mildner-Szkudlarz, H. H. Jeleń, and R. Zawirska-Wojtasiak, “The use of electronic and human nose for monitoring rapeseed oil autoxidation,” European Journal of Lipid Science and Technology, vol. 110, no. 1, pp. 61–72, 2008.
- F. Caporaso, J. D. Sink, P. S. Dimick, C. J. Mussinan, and A. Sanderson, “Volatile flavor constituents of ovine adipose tissue,” Journal of Agricultural and Food Chemistry, vol. 25, no. 6, pp. 1230–1234, 1977.
- E. Wong, C. B. Johnson, and L. N. Nixon, “The contribution of 4-methyloctanoic (hircinoic) acid to mutton and goat meat flavour,” New Zealand Journal of Agricultural Research, vol. 18, no. 3, pp. 261–266, 1975.
- H. Tian, P. Zhan, Z. Deng, H. Yan, and X. Zhu, “Development of a flavour fingerprint by GC-MS and GC-O combined with chemometric methods for the quality control of Korla pear (Pyrus serotina Reld),” International Journal of Food Science and Technology, vol. 49, no. 12, pp. 2546–2552, 2014.
Copyright © 2017 Ping Zhan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.