• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Rapid qualitative and quantitative analysis of strong aroma base liquor based on SPME-MS combined with chemometrics

    2021-05-20 08:53:28ZongboSunJunkuiLiJinfengWuXioboZouChiTngHoLimingLingXiojingYnXunZhou

    Zongbo Sun, Junkui Li, Jinfeng Wu, Xiobo Zou, Chi-Tng Ho,Liming Ling, Xiojing Yn, Xun Zhou,

    a School of Food and Biological Engineering, Jiangsu University, Zhenjiang 212013, China

    b Department of Food Science, Rutgers University, New Brunswick, New Jersey 08903, USA

    c Jiangsu King’s Luck Brewery Co. Ltd., Lianshui 223411, China

    Keywords:

    SPME-MS

    Strong aroma base liquor (SABL)

    Chemometrics

    Grade identification

    Ester compounds

    ABSTRACT

    To objectively classify and evaluate the strong aroma base liquors (SABLs) of different grades, solid-phase microextraction-mass spectrometry (SPME-MS) combined with chemometrics were used. Results showed that SPME-MS combined with a back-propagation artificial neural network (BPANN) method yielded almost the same recognition performance compared to linear discriminant analysis (LDA) in distinguishing different grades of SABL, with 84% recognition rate for the test set. Partial least squares (PLS), successive projection algorithm partial least squares (SPA-PLS) model, and competitive adaptive reweighed samplingpartial least squares (CARS-PLS) were established for the prediction of the four esters in the SABL.CARS-PLS model showed a greater advantage in the quantitative analysis of ethyl acetate, ethyl butyrate,ethyl caproate, and ethyl lactate. These results corroborated the hypothesis that SPME-MS combined with chemometrics can effectively achieve an accurate determination of different grades of SABL and prediction performance of esters.

    1. Introduction

    There are many brands of liquors globally such as whisky, vodka and gin among others. Liquors in China are very different from others in terms of flavor and taste. Chinese liquors are generally made with grains as main raw materials using four processing method, namely,cooking, saccharification, fermentation and distillation. These liquors make up the sixth-largest distilled spirits, as it is embedded in the Chinese long history of ‘liquor culture’. Almost all the liquors are made by blending different base liquors. The base liquors are unblended and made directly from the distillation of fermented grains with 60%-70% alcohol concentration.

    Strong aroma liquor (SAL) is one of the most important liquors cultivated in China, with the largest market share and the most consumed volume in China’s current liquor market. SAL belongs to a distilled liquor with their aroma compounds mainly composed of ethyl caproate. Moreover, this liquor is produced by using grains such as sorghum as the raw material, and processed by solid-state fermentation, storage and blending [1]. The grades of the SABL are mainly determined by following process. Firstly, the grade of the SABL preliminarily determined by experienced and skilled liquor workers during the distilling process as well as an evaluation based on preliminary grade identification. Secondly, different grades of SABL initially determined by workers are subjected to sensory analyses by professional panelists. During sensory evaluation,some SABLs which have a poor coordination or a peculiar smell are downgraded. The result of sensory evaluation, which is easily affected by subjective factors (experience, emotion, physiological conditions, and environmental factors) are full of uncertainty and limitation. Additionally, researchers have shown that the main aroma components in Chinese liquor are esters with a fruity aroma [2] that account for 60% of the total flavor components with a strong aroma.Ethyl caproate, ethyl lactate, ethyl acetate, and ethyl butyrate are the 4 main esters which play an important role in forming the characteristic flavor of strong aroma liquor. The content of the 4 esters need to meet relevant standards in different grades of SABL [3,4].

    In addition to the sensory evaluation being used for liquor classification of different grades and quality assessment, fluorescence spectra combined with simulated annealing algorithm are used for the prediction of mild aroma Chinese liquors. Yin et al. [5] have applied an electronic nose in the spirit classification. Near-infrared spectroscopy and ensemble classification have also been exploited for discrimination between authentic and adulterated liquors [6]. Although the above classification methods are considered to be faster, their poor stability, weak reproducibility and cross interference have limited their accuracy and broad application. Ester content can be generally detected by gas chromatography (GC) method. Notwithstanding,the GC method is tedious, time-consuming, and cannot achieve realtime and rapid detection. Therefore, it is necessary to establish a fast,accurate and stable classification and ester content determination method for the grade evaluation of SABL.

    A non-separative mass spectrometry (MS) method, based on a direct combination of the mass spectrometer and solid-phase microextraction (SPME) has been developed [7]. In the non-separative MS method, the chromatographic column of the gas chromatographymass spectrometry (GC-MS) instrument is replaced by a fused silica tube. The volatile portion of the sample is immediately desorbed at the injection port at a high temperature and subsequently injected into the ionization chamber without a chromatographic separation[8]. In spite of the limit from the poor volatility and high polarity of target analytes due to a poor partition and reduced recovery into the gas phase [9], SPME has unique advantages of a solvent-free, easy to operate and high reproducibility [10,11]. The use of solid-phase microextraction-mass spectrometry (SPME-MS) has become an effective and rapid detection technique [12-14]. In particular, SPMEMS coupled with one or more chemometrics is widely used in the food industry and has tremendous potential for classification and identification according to quality grades [15-18], geographical origin[19,20], pesticide residue [21], raw materials species [22] and others[23-25]. It has been applied for the quantification of benzoic acid in beverages for many years [27].

    Selecting the most appropriate method from a plethora of chemometrics techniques is very important. The method of principal component analysis (PCA), can eliminate unimportant data resulting in a dimensionality reduction. PCA is applied in the first stage of data processing to extract principal components (PCs), which are used for sample analyses and group visualization. Nevertheless, precise clustering in PCA cannot be provided [28,29]. In qualitative analysis,linear discriminant analysis (LDA), a supervised pattern recognition method is used to provide clustering information. The back propagation artificial neural network (BPANN) method is a multilayer forward feedback neural network used for prediction. In quantitative analysis, partial least squares (PLS), successive projection algorithmpartial least squares (SPA-PLS), and competitive adaptive reweighed sampling-partial least squares (CARS-PLS) can reduce time and effort in model calculation thus effectively improve the precision of the model.

    This study aimed to utilize the HS-SPME-MS method combined with multivariate analysis to discriminate grade quality and predict the ester content of SABL. The research goals were focused on 5 key points: (1) to obtain the ion intensity of SABL samples via HSSPME-MS technology; (2) to screen out irrelevant information using PCA, develop LDA and BPANN models, compare the two models and find out the optimal model for discriminating grade quality of SABL; (3) to obtain the true value of esters including ethyl caproate,ethyl lactate, ethyl acetate and ethyl butyrate by GC, combine the ion intensity of SABL samples with the true value of esters to build PLS model, analyze and predict performance; (4) to select optimal feature ions according to SPA and CARS, and used to develop esters prediction model; and (5) to compare the prediction performance of PLS, SPA-PLS, CARS-PLS models built by different ion selection methods and find out the optimal ions election method and model for predicting ester profiles in the SABL.

    2. Materials and methods

    2.1 Samples and chemicals

    Seventy-five SABL samples were collected from different production workshops of King’s Luck in the Huaian City of Jiangsu Province with about 70% ethanol content. All samples were stored at room temperature and the analysis conducted within 15 days after sample collection. Standards of ethyl acetate, ethyl butyrate, ethyl caproate, ethyl lactate and amyl acetate (internal standard) with 99.5%purity were purchased from Tianjin Guangfu Fine Chemical Research Institute (Tianjin, China).

    2.2 Grade identification of SABL

    The employees with experience in the process separated the SABL into 3 grades as it was produced: G1 (best), G2 (middle), G3(worse). Then, the expert assessors evaluated G1 and G2 according to its coordination, aroma, and peculiar smell. In G1 and G2, some unqualified SABL samples with poor coordination or a peculiar smell were downgraded. Each round of 6 samples was evaluated by 3 expert assessors that were selected randomly from 10 expert assessors.

    2.3 GC analysis

    The content of esters was measured in a GC (7890A, Agilent,Santa Clara, CA, USA) with a flame ionization detector (FID). The carrier gas was nitrogen with a flow rate of 2.0 mL/min. A total of 1 μL sample was injected in the split mode (1:40) into a CP-Wax 57CB column (50 m × 0.25 mm i.d., 0.2 μm film thickness, Varian, CA,USA). The programmed temperature was 35 °C for 3 min, increased to 90 °C at 4 °C/min, then amplified to 130 °C at 10 °C/min, finally ramped to 210 °C at 15 °C/min and held for 10 min. Injector and detector temperatures were 200 °C and 210 °C, respectively.

    2.4 SPME-MS analysis

    The SPME fiber (2 cm, 50/30 μm DVB/CAR/PDMS) was purchased from Sigma (St. Louis, MO, USA). Before using, the fiber was placed in GC injector port at a temperature of 270 °C to eliminate any probable contaminants.

    The SPME was carried out according to the following procedures.SABL (1 mL), 4 mL of distilled water and 2 g NaCl were placed in a 15 mL vial and tightly capped with a silica gel pad. Then the vial was placed on a magnetic stirrer (Corning, NY, USA) at 50 °C for 5 min and extracted for 20 min by the extraction fiber at the same temperature. After the extraction, the fiber was inserted into the inlet of GC and the components were desorbed for 5 min at 250 °C in the splitless mode.

    The SPME-MS was carried out using an GC (6890N, Agilent,USA) equipped with an MS detector (5973, Agilent, USA). A fused silica capillary column without a coating (1 m × 0.15 mm i.d, Agilent Technologies Inc., Shanghai, China) that replaced the capillary chromatographic column was inserted directly to MS. Helium (He) was chosen as the carrier gas at a constant flow rate (0.4 mL/min). The oven temperature was set to 200 °C to prevent the volatile components from recondensation and maintained for 5 min. The spectrometer operated in electron impact (EI) mode (70 eV). The ion source temperature was 230 °C, and the interface temperature was 250 °C. Detection was carried out in scan mode in a range ofm/z33-400.

    2.5 Chemometric analysis

    All the data processing and analysis were performed using MATLAB R2009a software. PCA was used as the first step of data analysis to visualize information and to detect patterns in data whereas LDA was used to calculate classification rules for sample discrimination. Also, BPANN was used to predict the SABL classification. Besides, PLS, SPA-PLS, and CARS-PLS were used to reduce the amount of calculation and improve the precision of the model in the quantitative analysis. Each data set was randomly divided into a training set (used to calculate the classification rules for two-thirds of the samples) and a test set (used to evaluate the prediction ability of rules and models for one-third of the samples).

    3. Results and discussion

    The SPME-MS method employed only MS to collect data within a few minutes to obtain the ion abundance. Then, an average mass spectrum (Fig. 1) including the abundance of allm/zof 33-400 ions were generated by HP-ChemStation system. Almost all the average mass spectrum had similarm/zions. For example,m/z45 ion was a base peak which is a characteristic of ethanol, andm/z88 ion can be ascribed to ethyl esters. Ethyl esters have a fruity aroma and it is of great importance in the quality assessment of liquors. Meanwhile,ethyl caproate, ethyl lactate, ethyl acetate and ethyl butyrate are the main esters for strong aroma liquor and were analyzed quantitatively by GC with n-amyl acetate as the internal standard. Firstly, these esters were identified by their respective retention times (RT) based on a mixture of their standards, and the concentration ranges of these four esters’ are shown in Table 1. The standard curve was calculated by concentration (X) of esters as theX-axis and the peak area ratio (Y)of the esters to the internal standard asY-axis. Table 1 showed thatXandYhad a linear relationship with a correlation coefficient greater than 0.999 6. Finally, the ester content was calculated according to the standard curve.

    Fig. 1 The average mass spectrum in the range of m/z 33-400 of the SABL generated by SPME-MS.

    Table 1Concentration range, retention times, standard curves and correlation coefficients of esters.

    3.1 Qualitative analysis of different grades of SABL

    The samples were divided into G1 (32 samples), G2 (24 samples),and G3 (19 samples). In G1 samples, 10 good samples (GG1) not only had a pure rich aroma and no peculiar smell but also tasted sweet and fragrant. However, 22 bad samples of G1 (BG1) had a bad smell(bitter, astringent, burnt, and other unacceptable smells).

    For G2 samples, 10 good samples (GG2) had a rich aroma,elegant and lasting and no off-smell and 14 bad samples (BG2) had also recorded an off-smell. Therefore, a total of 5 grades, GG1, BG1,GG2, BG2, and G3 were evaluated.

    It was difficult to intuitively differentiate various grades of the SABL according to the average mass spectrum although the ion abundance of every grade from the average mass spectrum was subtly different. To classify and evaluate the SABL of different grades objectively, the PCA method was applied in the first stage of data processing. This was followed by the LDA and BPANN model to identify the grade of SABL. From the average mass spectrum, the 75 samples were listed as rows and 368 ion abundance as columns,thus a 75 × 368 data matrix was generated while 50 samples for the training set and 25 samples for the test set were submitted to establish the model.

    PCA is one of the most frequently used chemometric methods and it is applied in the first stage of data processing. It works by replacing a large number of potentially correlated variables with a new linear set of fewer uncorrelated variables, namely PCs, which represent the majority of the original data. These new variables permit the construction of a multivariate model where it is possible to extract useful information from the original mass spectrum data by eliminating overlapping information [30]. Hence, PCA was implemented to determine the particular grade within the 5 grades of the SABL in this study.

    The PCA results (Fig. 2) could broadly reflect the distribution of the 5 grades of the SABL samples in the original variable space.Fig. 2 showed clearly a three-dimensional (3D) space of the SABL samples represented by PC1, PC2, and PC3. PC1, PC2 and PC3 showed 50.43%,18.29% and 10.96% of the variance, respectively.Segregation was observed only in the G3 samples but overlapping appeared in other grades. This might be duo to the G3 samples having some different qualities compared with the other grades. Therefore,PCA cannot de fine the boundaries of the other grades or used directly as a tool for determining the grade of the SABL sample unless the grade of the sample is in G3.

    Fig. 2 3-Dimension score plot of five grades of the SABL after PCA.

    To improve the classification of investigated samples, LDA which is probably the most frequently used technique was carried out to obtain a suitable classification of the SABL of various grades. LDA is a supervised classification technique based on the determination of linear discriminant functions which maximize the ratio of betweenclass variance and minimize the ratio of within-class variance [31].

    It was essential to use a PCA to determine a reduced-dimension of the mass spectrum data to extract optimal PCs as the input variable of the LDA identification model. The reduced-dimension data not only held enough of the original information but also achieved the compression of a certain amount of the original data which was available for subsequent processing. In Fig. 3A, the result of the training set and test set of different grades of the SABL based on SPME-MS analysis showed that the number of PCs affected the results of the training set and the test set in LDA. When the number of PCs was less than 7, the recognition rate of the model was not stable.As the number of PCs reached 7, the established LDA model obtained the best effect on the classification of SABL with discrimination rate of 86% for the training set and 84% rate for the test set. Besides, the minimal gap of the rate between training set and test set demonstrated good versatility of the LDA model with the principal component of 7.Although the prediction effect of the model was improved when the principal component is 15 compared to when the principal component is 7, it will increase the amount of data calculation and have the great possibility to cause over fitting. Additionally, from the 2-dimension score plot of 5 grades of the SABL after Fisher- LDA in Fig. 3B, the contribution rate of the first discriminant factor (LD1) and the second discriminant factor (LD2) were 88.98% and 9.82% respectively, and the cumulative contribution rate reached 98.80%. Samples of G3 were separated from all other samples. Additionally, discrimination between good samples and bad samples in the same grade was also satisfactory. Yet, LDA could not de fine the boundaries of the samples of GG1 and GG2 due to the existence of overlapping portions. In general, the proposed method improved the recognition performance of SABLs’ recognition performance.

    Fig. 3 (A) Discriminant result of LDA of the training and test set in the different PCs; (B) 2-Dimension score plot offive grades of the SABL after Fisher- LDA.

    BPANN, a familiar neural network model, is capable of dealing with many complex problems including nonlinear problems. BP is a commonly used learning algorithm in ANN applications, which uses the back-propagation algorithm as the gradient descent technique to minimize network error [32]. A neural network composes of an input layer, one or more hidden layers, and an output layer. Each layer in the BPANN has several neurons and each neuron transmits input values and processes to the next layer [33]. The number of neurons in the hidden layer could affect the error and convergence rate of the network to some extent. Similarly, the number of PCs also had a direct impact on the prediction performance of the model. Moreover, BPANN parameters such as the learning rate factor, the momentum factor, the target error and iteration times also have some influence on the performance of BPANN model. Therefore, it is crucial to optimize factors.

    The optimal BPANN model was achieved when the learning rate factor was set to 0.1, momentum factor at 0.7, target error from 10-8,and iteration times selected as 1 000. Table 2 shows the result of the training set and test set performed by BPANN using different numbers of PCs. As observed, when the number of PCs was increasing, the recognition rate of the training set and the test set always fluctuated.Compared to other BPANN results, better recognition performance was witnessed when the number of PCs and neurons were 4. Under those conditions, the recognition rate of the training set and the test set obtained 90% and 84% respectively. Hence, SPME-MS combined with the BPANN method yielded almost the same recognition performance compared to LDA in distinguishing different grades of SABL, with 84% recognition rate for the test set.

    Table 2Discriminant result of BPANN in the training set and test set under different PCs.

    3.2 Quantitative analysis of ester content

    The content of 4 esters (ethyl acetate, ethyl butyrate, ethyl caproate, and ethyl lactate) needs to meet corresponding GRADE criteria in different grades SABL. Therefore, a quick and accurate determination of the 4 esters plays a critical role in the grading of SABL. Table 3 showed the reference value of the ester content.The concentration range of these ester compounds was wide.This corroborates its high representativeness and could be used to establish a model.

    Table 3Reference values of esters of SABL samples in the training and test set.

    PLS is a well-known multivariate calibration method based on factor analysis and mostly used for multivariate analysis methods in quantitative analysis. It not only maps the original dependent variable data into a few latent variables whose information is very focused but also selects the latent variables which are greatly correlated with the dependent variables as the principal component variables to establish the model [34,35]. SPA can minimize the collinearity among variables and largely reduce the number of variables. By comparing the magnitude of the projected vector between different wavelengths,SPA takes the wavelength with the largest projection vector as the candidate wavelength and uses the robustness of the training model to determine the final characteristic wavelength. CARS, a feature wavelength selection method is based on Monte Carlo (MC) sampling and PLS regression coefficient [36,37]. CARS first establishes the corresponding PLS model through the sample of calibration set selected by MC sampling, establishes the absolute weight of the wavelength regression coefficient in this sampling, and removes the wavelength variable with smaller absolute value, The number of removed wavelength variable is then determined by exponential decrease function (EDF). Based on the remaining wavelength variables, the PLS model is established by selecting the wavelength using adaptive reweighted sampling (ARS). The corresponding wavelength of the PLS model with the smallest RMSECV is the selected characteristic wavelength of the PLS model [38].

    The best models are chosen based on the lowest root mean square error of cross-validation (RMSECV), the lowest root mean square error of prediction (RMSEP), and the highest correlation coefficients for the calibration data set (Rc) as well as the prediction data set(Rp). RMSECV is mainly used for the evaluation of feasibility and the predictive ability of the PLS, SPA-PLS, and CARS-PLS model.RMSEP is mainly used to evaluate the prediction ability of the model to the external sample. The size ofRshows the correlation between the predicted value and the measured value [39].

    The PLS, SPA-PLS, and CARS-PLS models were established for 4 ester compounds of different grades of the SABL. This study takes ethyl acetate as an example to analyze the established model.

    Fig. 4 Prediction results of the PLS model for ethyl acetate content. (A)Selection of optimum PCs. (B) Results of training set of PLS model based on full variables. (C) Results of test set of PLS model based on full variables.

    In the process of establishing the PLS model using full variables,two-thirds of the SABL samples were randomly selected as a training set and the other one-third of the trials were used as a test set. Fig. 4 demonstrated that the established PLS model yielded the best predictions, under the condition of utilizing the first 6 principal components by obtainingRc(0.828 8),Rp(0.777 8), RMSECV(102.00), and RMSEP (122.00). TheRcandRpgenerated by the PLS model established by full variables were less than 0.9, which demonstrated the ineffectiveness of the model in the prediction of the ethyl acetate in the SABL. This may be due to the existence of many ionic variables that are not related to ethyl acetate content which resulted in the degradation of prediction performance and screening of characteristic variables appeared to be a feasible strategy to solve this problem.

    The range of the number of characteristic variables was set to 1-18 when using SPA to select characteristic variables that were determined by the root mean square error (RMSE). The selection result is shown in Fig. 5A. When the number of variables was 13, the obtained RMSE value was the smallest. The 13 variables selected by SPA were used to develop the PLS model. Fig. 5B-5D illustrated that the established SPA-PLS model yielded the best predictions under the condition of utilizing 9 principal components by achieving Rc(0.848 4),Rp(0.803 2), RMSECV (96.90), and RMSEP (110.00). Compared to the PLS model, SPA-PLS obtained a higher Rcand Rp.

    Fig. 5 Prediction results of SPA-PLS model for ethyl acetate content. (A)Process of screening characteristic variables by SPA. (B) Selection of optimum PCs. (C) Results of training set of PLS model based on SPA. (D) Results of test set of PLS model based on SPA.

    Fig. 6 Prediction results of CARS-PLS model for ethyl acetate content.(A) Process of screening characteristic variables by CARS. (B) Selection of optimum PCs. (C) Results of training set of PLS model based on CARS. (D)Results of test set of PLS model based on CARS.

    As shown in Fig. 6A1, A2, A3, after sampling the MC (19 times), 35 optimal variables were screened out to establish the PLS model. As shown in Fig. 6B, with the growth of number of PLS components, RMSECV changed in a ‘V’ shape and finally achieved the minimum value when PLS components were 7 and the established CARS-PLS model yielded the best predictions with the RMSECV (51.50) of the training set,Rcof 0.956 3,RMSEP of 52.90 for the test set, andRpof 0.937 0 in the test set under this condition, which showed in Fig. 6C, 6D. Compared to PLS and SPA-PLS model, CARS-PLS shown greater probability of models to predict ethyl acetate in SABL. In addition, 30 PLS models developed by randomly selected sample set which have the same characteristic variables asm/z36, 41, 50, 53, 54, 57, 61, 76, 81, 85, 90, 100, 103,112, 114, 115, 116, 118, 129, 133, 134, 150, 151, 173, 186, 189, 191,197, 251, 255, 264, 283, 285, 341 and 342. Statistical analysis showed that the averageRcof 0.955 8, andRpof 0.928 6 and RMSECV of 50.75, RMSEP of 63.38, which may convey an important message that the PLS model developed by characteristic variables selected from CARS possess excellent stability and high fitting performance.

    The prediction results from the model developed by PLS, SPAPLS, and CARS-PLS in ethyl butyrate, ethyl caproate, and ethyl lactate are shown in Table 4. CARS-PLS possessed the best prediction result amongst these 4 esters followed by SPA-PLS except ethyl lactate. After variable selection, the CARS-PLS model displayed elevation prediction performance while SPA-PLS showed degradation compared to the PLS model built with full variables in ethyl lactate.This may indicate that the SPA removed some variables related to the content of ethyl lactate in the process of variable selection. This inference was also attested by the phenomenon that, the characteristic variables selected by CARS are significantly higher than by SPA in screening the results of the characteristic variables of the other three esters. Besides, it is worth noting that PLS model combined with CARS showed the lowest RMSECV and RMSEP compared with PLS and SPA-PLS, which demonstrated that the model established by the variables extracted by CARS is closer to the true value of the predicted value of the ester content in the base liquor. Therefore, CARS as a variable screening method can be well applied to the establishment of a prediction model for the characteristic ester content of SABL.

    Table 4The optimal model of the four ester compounds.

    4. Conclusion

    The SPME-MS technique combined with chemometrics resulted in high classification accuracy for the SABL from 5 different sensory grades. The results indicated that SPME-MS combined with the BPANN method yielded almost the same recognition performance juxtaposed to linear discriminant analysis (LDA) in distinguishing different grades of SABL with 84% of the recognition rate of the test set.Additionally, it was found that the CARS-PLS model showed a greater advantage in the quantitative analysis of ethyl acetate, ethyl butyrate,ethyl caproate, and ethyl lactate compared to PLS and SPA-PLS model.

    Therefore, by overcoming the complex and time-consuming problems of the traditional detection technique, SPME-MS technology combined with chemometrics successfully demonstrated as a tool for distinguishing the quality of different grades of the SABL and quantitative analysis of four main ester compounds.

    conflict of interest

    The authors declare no conflict of interest.

    Acknowledgments

    The study was supported by the Key Research and Development Program of Jiangsu Province (BE2020312), National Natural Science Foundation of China (31671844), Open Project of National Engineering Laboratory for Agri-product Quality Traceability (AQT-2019-YB7), Science Foundation for Postdoctoral in Jiangsu Province(1501100C), Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

    国产一区二区在线观看av| 亚洲av成人精品一区久久| 九九久久精品国产亚洲av麻豆| 日本黄色日本黄色录像| 亚洲人成网站在线播| 国产69精品久久久久777片| 91久久精品国产一区二区三区| 观看av在线不卡| 夜夜看夜夜爽夜夜摸| 国产成人午夜福利电影在线观看| tube8黄色片| 啦啦啦啦在线视频资源| 啦啦啦中文免费视频观看日本| 22中文网久久字幕| 美女大奶头黄色视频| av播播在线观看一区| 日本免费在线观看一区| 国产精品国产三级专区第一集| 哪个播放器可以免费观看大片| 色婷婷av一区二区三区视频| 欧美性感艳星| av女优亚洲男人天堂| 中文字幕人妻熟人妻熟丝袜美| 在线观看免费视频网站a站| 少妇的逼水好多| 国产成人精品无人区| 欧美日韩在线观看h| 我的女老师完整版在线观看| 亚洲色图综合在线观看| 视频中文字幕在线观看| 啦啦啦啦在线视频资源| 80岁老熟妇乱子伦牲交| 国产有黄有色有爽视频| 日韩中字成人| 丰满乱子伦码专区| 一本久久精品| 亚洲美女搞黄在线观看| 亚洲丝袜综合中文字幕| 亚洲精品国产av蜜桃| 国产精品久久久久久精品电影小说| 成人二区视频| tube8黄色片| 久久av网站| .国产精品久久| 亚洲精品视频女| 美女国产视频在线观看| 国产av国产精品国产| 在线 av 中文字幕| 国产日韩欧美亚洲二区| 伦精品一区二区三区| 美女脱内裤让男人舔精品视频| 18在线观看网站| av福利片在线| 青春草亚洲视频在线观看| 亚洲高清免费不卡视频| 性色av一级| 日日啪夜夜爽| 一个人看视频在线观看www免费| 久久人妻熟女aⅴ| 久久久亚洲精品成人影院| 26uuu在线亚洲综合色| 欧美精品国产亚洲| 国语对白做爰xxxⅹ性视频网站| 在线免费观看不下载黄p国产| 丝瓜视频免费看黄片| 久久毛片免费看一区二区三区| 久久热精品热| 18禁裸乳无遮挡动漫免费视频| 国产高清国产精品国产三级| 久久国内精品自在自线图片| 亚洲精品456在线播放app| 成人综合一区亚洲| 狠狠婷婷综合久久久久久88av| 91aial.com中文字幕在线观看| 亚洲精品色激情综合| 91在线精品国自产拍蜜月| 国产精品嫩草影院av在线观看| 久久精品久久久久久久性| 丁香六月天网| .国产精品久久| 免费人妻精品一区二区三区视频| 精品久久久久久久久av| 亚洲精品国产色婷婷电影| 校园人妻丝袜中文字幕| 激情五月婷婷亚洲| 国产高清三级在线| 久久人人爽人人片av| 亚洲不卡免费看| 亚洲av免费高清在线观看| 久久青草综合色| 热99国产精品久久久久久7| 女的被弄到高潮叫床怎么办| 日日啪夜夜爽| 成人综合一区亚洲| 人人妻人人澡人人爽人人夜夜| 母亲3免费完整高清在线观看 | 久久久久人妻精品一区果冻| 成年人午夜在线观看视频| 特大巨黑吊av在线直播| 亚洲性久久影院| 中文天堂在线官网| 久久久久久久久大av| 国产成人av激情在线播放 | 国产日韩欧美在线精品| 女性生殖器流出的白浆| 男人操女人黄网站| 夫妻午夜视频| av不卡在线播放| 午夜精品国产一区二区电影| 熟女av电影| 制服丝袜香蕉在线| 18禁裸乳无遮挡动漫免费视频| av国产精品久久久久影院| 少妇的逼好多水| 男男h啪啪无遮挡| 久久99热6这里只有精品| av在线老鸭窝| 视频在线观看一区二区三区| 狂野欧美激情性xxxx在线观看| 老司机亚洲免费影院| 黑人巨大精品欧美一区二区蜜桃 | 性高湖久久久久久久久免费观看| 亚洲精品日韩在线中文字幕| av国产精品久久久久影院| 99热这里只有是精品在线观看| 欧美一级a爱片免费观看看| 国产精品久久久久久久电影| 日韩精品有码人妻一区| 成人漫画全彩无遮挡| 少妇人妻精品综合一区二区| 国产午夜精品一二区理论片| 国产成人精品久久久久久| 亚洲av成人精品一二三区| 国产成人av激情在线播放 | 中文字幕最新亚洲高清| 麻豆乱淫一区二区| 日韩一区二区视频免费看| 国产av码专区亚洲av| 精品国产一区二区久久| 一区二区三区四区激情视频| 蜜桃久久精品国产亚洲av| 欧美日韩综合久久久久久| av免费观看日本| 久久精品久久久久久噜噜老黄| 黄色怎么调成土黄色| 国产乱来视频区| 亚洲精品久久成人aⅴ小说 | 一区二区三区免费毛片| 999精品在线视频| 亚洲丝袜综合中文字幕| 欧美97在线视频| 在线看a的网站| 精品午夜福利在线看| 免费av中文字幕在线| 亚洲国产精品成人久久小说| 午夜免费观看性视频| 午夜福利在线观看免费完整高清在| 桃花免费在线播放| 大片免费播放器 马上看| 飞空精品影院首页| 大香蕉久久网| 99热网站在线观看| 最后的刺客免费高清国语| 内地一区二区视频在线| 国产男女超爽视频在线观看| 欧美bdsm另类| 亚洲,一卡二卡三卡| 亚洲国产欧美日韩在线播放| 国产精品一区二区三区四区免费观看| 亚洲欧洲国产日韩| 少妇猛男粗大的猛烈进出视频| 18禁动态无遮挡网站| 亚洲精品美女久久av网站| 成人黄色视频免费在线看| 国产日韩一区二区三区精品不卡 | 乱人伦中国视频| 国产在线视频一区二区| 日韩制服骚丝袜av| 99久国产av精品国产电影| 91成人精品电影| 免费大片黄手机在线观看| 国产午夜精品久久久久久一区二区三区| 夫妻午夜视频| 狂野欧美白嫩少妇大欣赏| 考比视频在线观看| 蜜臀久久99精品久久宅男| 搡老乐熟女国产| 国产成人aa在线观看| 精品一品国产午夜福利视频| 日韩人妻高清精品专区| av天堂久久9| 热re99久久国产66热| 一二三四中文在线观看免费高清| 欧美日韩视频高清一区二区三区二| av卡一久久| 欧美成人精品欧美一级黄| 欧美日韩亚洲高清精品| 国产午夜精品一二区理论片| 亚洲av国产av综合av卡| 久久人人爽av亚洲精品天堂| 久久99热6这里只有精品| 欧美老熟妇乱子伦牲交| 美女主播在线视频| 在线观看美女被高潮喷水网站| 日本黄色片子视频| 另类亚洲欧美激情| 免费观看a级毛片全部| 午夜免费鲁丝| 国产爽快片一区二区三区| 五月开心婷婷网| 观看av在线不卡| av国产精品久久久久影院| 欧美成人精品欧美一级黄| 国产在线视频一区二区| 少妇被粗大猛烈的视频| 啦啦啦视频在线资源免费观看| 老司机亚洲免费影院| 99久久人妻综合| 两个人免费观看高清视频| 少妇高潮的动态图| 七月丁香在线播放| 亚洲av免费高清在线观看| 精品一品国产午夜福利视频| 汤姆久久久久久久影院中文字幕| 黄色怎么调成土黄色| 三上悠亚av全集在线观看| 国产在线免费精品| 日本wwww免费看| 中文欧美无线码| 91在线精品国自产拍蜜月| 免费黄网站久久成人精品| 亚洲久久久国产精品| 欧美xxⅹ黑人| 亚洲高清免费不卡视频| 成年女人在线观看亚洲视频| 男人操女人黄网站| 亚洲国产最新在线播放| 欧美日本中文国产一区发布| 搡女人真爽免费视频火全软件| 美女内射精品一级片tv| 秋霞伦理黄片| 18在线观看网站| 高清欧美精品videossex| 免费少妇av软件| 日韩av在线免费看完整版不卡| 激情五月婷婷亚洲| 晚上一个人看的免费电影| 免费大片黄手机在线观看| 人人妻人人澡人人看| 最近手机中文字幕大全| 天堂俺去俺来也www色官网| 尾随美女入室| 国产伦精品一区二区三区视频9| 王馨瑶露胸无遮挡在线观看| 免费黄频网站在线观看国产| 人妻夜夜爽99麻豆av| 九色成人免费人妻av| 精品人妻熟女av久视频| 精品人妻熟女毛片av久久网站| 日韩伦理黄色片| 久久午夜福利片| 女的被弄到高潮叫床怎么办| 精品酒店卫生间| 国产午夜精品久久久久久一区二区三区| 亚洲综合色网址| 在线 av 中文字幕| 伊人久久精品亚洲午夜| 亚洲欧美成人综合另类久久久| 亚洲国产最新在线播放| videosex国产| 免费久久久久久久精品成人欧美视频 | 国产免费福利视频在线观看| 色哟哟·www| 超碰97精品在线观看| 久久女婷五月综合色啪小说| 亚洲国产av新网站| 97精品久久久久久久久久精品| 欧美人与善性xxx| 母亲3免费完整高清在线观看 | 国产精品一区二区在线不卡| 性色av一级| 在线精品无人区一区二区三| 色吧在线观看| 九色亚洲精品在线播放| 免费看光身美女| 国产国语露脸激情在线看| 精品卡一卡二卡四卡免费| 亚洲第一av免费看| 欧美人与性动交α欧美精品济南到 | 久久99热这里只频精品6学生| 麻豆精品久久久久久蜜桃| 岛国毛片在线播放| 亚洲综合色网址| 最后的刺客免费高清国语| 亚洲综合精品二区| 免费看光身美女| 久久午夜综合久久蜜桃| 久久久欧美国产精品| 一本色道久久久久久精品综合| 特大巨黑吊av在线直播| av免费观看日本| 麻豆乱淫一区二区| freevideosex欧美| 亚洲欧美日韩另类电影网站| 日韩av在线免费看完整版不卡| 在现免费观看毛片| 狠狠婷婷综合久久久久久88av| 高清不卡的av网站| 不卡视频在线观看欧美| 中文字幕最新亚洲高清| 97超碰精品成人国产| 2018国产大陆天天弄谢| 欧美一级a爱片免费观看看| www.色视频.com| 国产欧美日韩一区二区三区在线 | 国产片特级美女逼逼视频| 国产男女超爽视频在线观看| 午夜91福利影院| 性高湖久久久久久久久免费观看| 卡戴珊不雅视频在线播放| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品久久久久久久电影| 精品少妇内射三级| 欧美最新免费一区二区三区| 啦啦啦视频在线资源免费观看| 国产色婷婷99| 日本色播在线视频| 大香蕉97超碰在线| 亚洲无线观看免费| 最新中文字幕久久久久| a级毛片黄视频| 精品人妻熟女毛片av久久网站| 少妇的逼好多水| 国产精品久久久久久久久免| 男女啪啪激烈高潮av片| 亚洲欧美一区二区三区国产| 国产永久视频网站| av女优亚洲男人天堂| 亚洲第一区二区三区不卡| 性色avwww在线观看| 免费av不卡在线播放| 欧美少妇被猛烈插入视频| 亚洲精品中文字幕在线视频| 国产熟女欧美一区二区| 狂野欧美激情性xxxx在线观看| 免费黄网站久久成人精品| 永久网站在线| 日本猛色少妇xxxxx猛交久久| 嘟嘟电影网在线观看| 汤姆久久久久久久影院中文字幕| 亚洲欧美日韩另类电影网站| 免费人妻精品一区二区三区视频| 婷婷色av中文字幕| 男人操女人黄网站| 18禁动态无遮挡网站| 亚洲精品日韩av片在线观看| 中文字幕制服av| 中文欧美无线码| 国产视频内射| 一区在线观看完整版| 韩国高清视频一区二区三区| 日本色播在线视频| 久久久久久人妻| 免费播放大片免费观看视频在线观看| 乱人伦中国视频| 亚洲国产av影院在线观看| 精品99又大又爽又粗少妇毛片| 国产精品一国产av| 妹子高潮喷水视频| 国产精品国产三级专区第一集| 亚洲欧美一区二区三区国产| 欧美亚洲 丝袜 人妻 在线| 一级黄片播放器| 亚洲一区二区三区欧美精品| 一本一本综合久久| 欧美精品一区二区大全| 亚洲欧美一区二区三区国产| 高清黄色对白视频在线免费看| 2021少妇久久久久久久久久久| 亚洲性久久影院| 亚洲综合色网址| 亚洲欧美色中文字幕在线| 国产极品粉嫩免费观看在线 | 中文字幕人妻熟人妻熟丝袜美| 久久女婷五月综合色啪小说| 天堂8中文在线网| 欧美97在线视频| 美女主播在线视频| 亚洲精品第二区| 插逼视频在线观看| 精品久久久久久久久av| 午夜视频国产福利| 18禁裸乳无遮挡动漫免费视频| 欧美成人精品欧美一级黄| 一级,二级,三级黄色视频| 2021少妇久久久久久久久久久| 看十八女毛片水多多多| 亚洲精品国产色婷婷电影| 91成人精品电影| 精品久久久久久久久亚洲| 美女国产高潮福利片在线看| 国产成人精品无人区| 能在线免费看毛片的网站| 亚洲国产日韩一区二区| 久久影院123| 美女中出高潮动态图| 国产伦精品一区二区三区视频9| 在线观看www视频免费| 日本色播在线视频| 热re99久久精品国产66热6| 日韩一本色道免费dvd| 亚洲四区av| 国产男女超爽视频在线观看| 少妇猛男粗大的猛烈进出视频| 麻豆成人av视频| 一级片'在线观看视频| 婷婷色av中文字幕| 国产伦精品一区二区三区视频9| 亚洲精品中文字幕在线视频| 这个男人来自地球电影免费观看 | av视频免费观看在线观看| 亚洲av.av天堂| 成年av动漫网址| 欧美精品一区二区免费开放| 99久国产av精品国产电影| 精品熟女少妇av免费看| av国产久精品久网站免费入址| 免费大片黄手机在线观看| 黑人欧美特级aaaaaa片| 午夜免费男女啪啪视频观看| 18禁裸乳无遮挡动漫免费视频| 高清毛片免费看| 国产午夜精品久久久久久一区二区三区| 久久精品夜色国产| 日韩一本色道免费dvd| 91久久精品国产一区二区成人| 国产成人午夜福利电影在线观看| 91在线精品国自产拍蜜月| 伦精品一区二区三区| 久久ye,这里只有精品| 国产精品成人在线| 免费观看性生交大片5| av在线app专区| 亚洲国产精品一区二区三区在线| 日韩av免费高清视频| 建设人人有责人人尽责人人享有的| 国产色婷婷99| 国产极品天堂在线| 91久久精品国产一区二区成人| 91精品国产九色| 欧美激情国产日韩精品一区| 亚洲国产精品专区欧美| 中文乱码字字幕精品一区二区三区| 80岁老熟妇乱子伦牲交| 26uuu在线亚洲综合色| 婷婷色综合www| 80岁老熟妇乱子伦牲交| 美女主播在线视频| 精品久久久精品久久久| 久热这里只有精品99| 欧美人与性动交α欧美精品济南到 | 亚洲美女黄色视频免费看| 亚洲精品成人av观看孕妇| 亚洲经典国产精华液单| 哪个播放器可以免费观看大片| 3wmmmm亚洲av在线观看| 亚洲综合色网址| 美女福利国产在线| av在线老鸭窝| 极品人妻少妇av视频| 亚洲不卡免费看| 国产 一区精品| 成人二区视频| 男人操女人黄网站| 免费黄色在线免费观看| 亚洲av欧美aⅴ国产| 丰满少妇做爰视频| 久久久久久久久久久免费av| 日本av手机在线免费观看| 欧美激情国产日韩精品一区| 久久综合国产亚洲精品| 欧美人与善性xxx| 国产综合精华液| 不卡视频在线观看欧美| xxxhd国产人妻xxx| 人人妻人人爽人人添夜夜欢视频| 精品久久久久久久久av| 亚洲丝袜综合中文字幕| 国产成人91sexporn| 国产一级毛片在线| 欧美三级亚洲精品| av电影中文网址| 精品人妻熟女毛片av久久网站| 伦理电影免费视频| 亚洲三级黄色毛片| 制服诱惑二区| 日韩 亚洲 欧美在线| 国产精品一区二区在线观看99| 天堂8中文在线网| 中文字幕人妻熟人妻熟丝袜美| 国产黄色免费在线视频| 国产精品嫩草影院av在线观看| 久久国产精品大桥未久av| 成人漫画全彩无遮挡| 精品一品国产午夜福利视频| 色吧在线观看| 久久久久久久国产电影| 久久久午夜欧美精品| 99久国产av精品国产电影| 看免费成人av毛片| 一级a做视频免费观看| 久久久久久久久久久丰满| 欧美人与性动交α欧美精品济南到 | 香蕉精品网在线| av免费观看日本| 国产精品蜜桃在线观看| 精品久久久久久久久亚洲| 中文字幕制服av| 精品国产一区二区久久| 久久久国产欧美日韩av| 亚洲欧美成人综合另类久久久| 高清毛片免费看| 免费不卡的大黄色大毛片视频在线观看| 亚洲欧美色中文字幕在线| 又粗又硬又长又爽又黄的视频| h视频一区二区三区| 亚洲欧洲日产国产| 尾随美女入室| 久久97久久精品| 国产av一区二区精品久久| 亚洲色图 男人天堂 中文字幕 | 欧美性感艳星| 国产精品99久久久久久久久| 一区二区三区乱码不卡18| 成年美女黄网站色视频大全免费 | 亚洲精品aⅴ在线观看| 老司机亚洲免费影院| 美女大奶头黄色视频| 一边摸一边做爽爽视频免费| 久久99精品国语久久久| 另类精品久久| 又大又黄又爽视频免费| 极品少妇高潮喷水抽搐| 黑人高潮一二区| 精品一区二区三卡| 99久久中文字幕三级久久日本| 久久久国产欧美日韩av| 午夜激情久久久久久久| 精品视频人人做人人爽| 国产成人a∨麻豆精品| 欧美日韩综合久久久久久| 亚洲四区av| 22中文网久久字幕| 国产亚洲一区二区精品| 精品国产乱码久久久久久小说| av黄色大香蕉| 五月玫瑰六月丁香| 另类亚洲欧美激情| 国产黄片视频在线免费观看| 五月开心婷婷网| 国产在视频线精品| 伊人亚洲综合成人网| 国产视频内射| 亚洲欧美成人综合另类久久久| 久久精品国产鲁丝片午夜精品| 久久久久久久久大av| 国产不卡av网站在线观看| 黄片播放在线免费| 色94色欧美一区二区| 精品久久久精品久久久| av网站免费在线观看视频| 99久久人妻综合| 18禁在线播放成人免费| 久久久国产精品麻豆| 人妻人人澡人人爽人人| 国产女主播在线喷水免费视频网站| 成年av动漫网址| 亚州av有码| 少妇被粗大的猛进出69影院 | 国产成人精品婷婷| 午夜视频国产福利| xxxhd国产人妻xxx| 色网站视频免费| 亚洲国产精品成人久久小说| 婷婷色麻豆天堂久久| 在线亚洲精品国产二区图片欧美 | 在线亚洲精品国产二区图片欧美 | 婷婷成人精品国产| 69精品国产乱码久久久| 日韩视频在线欧美| 成人漫画全彩无遮挡| 久久久国产精品麻豆| 国产精品嫩草影院av在线观看| 亚洲精品自拍成人| 超碰97精品在线观看| 国产成人91sexporn| 久久国产精品男人的天堂亚洲 | 国产亚洲最大av| kizo精华| 少妇猛男粗大的猛烈进出视频| 天堂中文最新版在线下载| 亚洲精品成人av观看孕妇| 亚洲av成人精品一二三区| 青春草视频在线免费观看| 亚洲欧洲日产国产| 日韩视频在线欧美| 免费不卡的大黄色大毛片视频在线观看| 黑人猛操日本美女一级片| av电影中文网址| 国产男女超爽视频在线观看| 亚洲,一卡二卡三卡| 精品人妻熟女毛片av久久网站| 国模一区二区三区四区视频| 国产成人精品无人区| a级毛色黄片|