• <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).

    国产欧美日韩一区二区三 | 伦理电影免费视频| 亚洲国产精品一区三区| 女人爽到高潮嗷嗷叫在线视频| 曰老女人黄片| 丝瓜视频免费看黄片| 一二三四在线观看免费中文在| 国产99久久九九免费精品| 一本大道久久a久久精品| 国产视频一区二区在线看| 国产深夜福利视频在线观看| 亚洲久久久国产精品| 51午夜福利影视在线观看| 亚洲人成电影观看| 日韩中文字幕视频在线看片| 亚洲成人手机| 精品高清国产在线一区| 下体分泌物呈黄色| 狂野欧美激情性bbbbbb| 亚洲精品国产av蜜桃| 在线观看免费视频网站a站| 欧美黑人欧美精品刺激| 最近最新中文字幕大全免费视频| 91av网站免费观看| 婷婷丁香在线五月| 精品国产国语对白av| 精品少妇一区二区三区视频日本电影| 国产亚洲av片在线观看秒播厂| 宅男免费午夜| 国产精品自产拍在线观看55亚洲 | 在线永久观看黄色视频| 中文字幕人妻丝袜制服| 后天国语完整版免费观看| 久久久久久亚洲精品国产蜜桃av| 国产精品一区二区免费欧美 | 满18在线观看网站| 日韩一卡2卡3卡4卡2021年| 伦理电影免费视频| 欧美性长视频在线观看| 男女午夜视频在线观看| 日韩三级视频一区二区三区| 国产一区二区三区av在线| 成人国产av品久久久| 亚洲精品国产色婷婷电影| 91九色精品人成在线观看| 性色av一级| 在线观看一区二区三区激情| 99国产综合亚洲精品| 国产精品成人在线| 国产高清videossex| 亚洲专区字幕在线| 亚洲少妇的诱惑av| 日本黄色日本黄色录像| 色精品久久人妻99蜜桃| 人妻久久中文字幕网| 中文字幕av电影在线播放| 一级a爱视频在线免费观看| 99久久人妻综合| 欧美黑人精品巨大| 男女之事视频高清在线观看| 免费不卡黄色视频| 大片免费播放器 马上看| 一级a爱视频在线免费观看| www日本在线高清视频| 妹子高潮喷水视频| 亚洲av美国av| 嫁个100分男人电影在线观看| 成年女人毛片免费观看观看9 | 欧美变态另类bdsm刘玥| 男女下面插进去视频免费观看| 满18在线观看网站| 日韩一卡2卡3卡4卡2021年| 国产成人av激情在线播放| 成年人午夜在线观看视频| 成人国产一区最新在线观看| 国产一区二区激情短视频 | 欧美性长视频在线观看| 国产在线一区二区三区精| 免费在线观看完整版高清| 天天操日日干夜夜撸| 电影成人av| 免费不卡黄色视频| 一本综合久久免费| 老熟妇仑乱视频hdxx| 人妻 亚洲 视频| 美女脱内裤让男人舔精品视频| 国产精品偷伦视频观看了| 国产成人啪精品午夜网站| 女人爽到高潮嗷嗷叫在线视频| 国产欧美日韩一区二区三区在线| 一二三四社区在线视频社区8| 一级片免费观看大全| 首页视频小说图片口味搜索| 欧美乱码精品一区二区三区| 午夜免费成人在线视频| 久久精品人人爽人人爽视色| 免费高清在线观看视频在线观看| 成人18禁高潮啪啪吃奶动态图| 在线 av 中文字幕| videos熟女内射| 国产三级黄色录像| 久久久精品国产亚洲av高清涩受| 美女福利国产在线| 丝瓜视频免费看黄片| 亚洲国产日韩一区二区| 久9热在线精品视频| 搡老乐熟女国产| 国产精品久久久久久精品电影小说| 国产精品国产三级国产专区5o| 成年动漫av网址| 亚洲七黄色美女视频| 精品国产一区二区久久| 91精品伊人久久大香线蕉| 超碰97精品在线观看| 男女无遮挡免费网站观看| 大陆偷拍与自拍| 日本一区二区免费在线视频| 久久国产亚洲av麻豆专区| 老熟妇仑乱视频hdxx| 一本一本久久a久久精品综合妖精| 性色av乱码一区二区三区2| bbb黄色大片| 国产精品一区二区免费欧美 | 美女脱内裤让男人舔精品视频| 久9热在线精品视频| 亚洲情色 制服丝袜| 777久久人妻少妇嫩草av网站| 丝袜喷水一区| 男女之事视频高清在线观看| 成年女人毛片免费观看观看9 | 深夜精品福利| 麻豆国产av国片精品| 黑人欧美特级aaaaaa片| 超碰97精品在线观看| 久久av网站| 亚洲人成电影观看| 亚洲欧美日韩高清在线视频 | 少妇精品久久久久久久| 搡老岳熟女国产| 欧美激情高清一区二区三区| 午夜激情久久久久久久| 王馨瑶露胸无遮挡在线观看| 久久亚洲精品不卡| 亚洲精品久久久久久婷婷小说| 免费av中文字幕在线| 老司机亚洲免费影院| 99国产极品粉嫩在线观看| 国产一区二区激情短视频 | 在线 av 中文字幕| 在线永久观看黄色视频| 亚洲国产欧美一区二区综合| 一本色道久久久久久精品综合| 亚洲中文日韩欧美视频| 无遮挡黄片免费观看| 交换朋友夫妻互换小说| 午夜福利视频精品| 日韩欧美一区视频在线观看| 超色免费av| 国产精品秋霞免费鲁丝片| 18禁裸乳无遮挡动漫免费视频| 极品人妻少妇av视频| 女人精品久久久久毛片| 在线永久观看黄色视频| 欧美日韩亚洲高清精品| 人妻一区二区av| 人人妻人人爽人人添夜夜欢视频| 亚洲av电影在线进入| 视频在线观看一区二区三区| 这个男人来自地球电影免费观看| 久久精品亚洲av国产电影网| 18禁裸乳无遮挡动漫免费视频| 丰满迷人的少妇在线观看| 精品一区二区三区av网在线观看 | 国产福利在线免费观看视频| tocl精华| 三级毛片av免费| 国产高清videossex| 亚洲精品国产av成人精品| 国产无遮挡羞羞视频在线观看| 国产在线一区二区三区精| 宅男免费午夜| 国产伦人伦偷精品视频| 国产精品亚洲av一区麻豆| 精品卡一卡二卡四卡免费| 自线自在国产av| 免费观看av网站的网址| 后天国语完整版免费观看| 91老司机精品| 女人被躁到高潮嗷嗷叫费观| 岛国毛片在线播放| 亚洲精品自拍成人| 水蜜桃什么品种好| 国产精品一二三区在线看| 丝袜脚勾引网站| 亚洲中文字幕日韩| 极品少妇高潮喷水抽搐| 亚洲,欧美精品.| 99国产综合亚洲精品| 中文欧美无线码| 大片电影免费在线观看免费| 久久久久国产一级毛片高清牌| 男女高潮啪啪啪动态图| 99久久综合免费| 天堂俺去俺来也www色官网| 午夜福利一区二区在线看| 男女床上黄色一级片免费看| 岛国毛片在线播放| 免费高清在线观看视频在线观看| 一本色道久久久久久精品综合| 色视频在线一区二区三区| 亚洲av日韩精品久久久久久密| 亚洲视频免费观看视频| 国产在线免费精品| 天堂中文最新版在线下载| 精品福利永久在线观看| 日本一区二区免费在线视频| svipshipincom国产片| 久久精品国产综合久久久| 老汉色∧v一级毛片| 青春草视频在线免费观看| 亚洲精品粉嫩美女一区| 高清黄色对白视频在线免费看| 麻豆国产av国片精品| 免费不卡黄色视频| 女人精品久久久久毛片| 看免费av毛片| 手机成人av网站| 伊人久久大香线蕉亚洲五| 久久毛片免费看一区二区三区| 国产精品自产拍在线观看55亚洲 | 永久免费av网站大全| 成人国产av品久久久| 午夜福利,免费看| 亚洲人成77777在线视频| 亚洲色图 男人天堂 中文字幕| 亚洲第一欧美日韩一区二区三区 | 女性生殖器流出的白浆| 亚洲激情五月婷婷啪啪| 国产精品偷伦视频观看了| 日韩中文字幕视频在线看片| 久久久久久久精品精品| 97精品久久久久久久久久精品| 亚洲精品第二区| 巨乳人妻的诱惑在线观看| 久久久水蜜桃国产精品网| 又大又爽又粗| 国产精品秋霞免费鲁丝片| 久久久久久久久免费视频了| 中国国产av一级| 欧美精品亚洲一区二区| 亚洲中文字幕日韩| 极品少妇高潮喷水抽搐| 精品少妇黑人巨大在线播放| 永久免费av网站大全| 一区二区av电影网| 夫妻午夜视频| 色婷婷久久久亚洲欧美| 国产xxxxx性猛交| 美女大奶头黄色视频| www.熟女人妻精品国产| 久热爱精品视频在线9| 蜜桃在线观看..| 国产精品成人在线| 搡老乐熟女国产| 欧美日韩黄片免| 19禁男女啪啪无遮挡网站| 视频区图区小说| 美国免费a级毛片| 伊人久久大香线蕉亚洲五| 欧美午夜高清在线| 69av精品久久久久久 | 国产男女超爽视频在线观看| 亚洲天堂av无毛| 欧美中文综合在线视频| 我要看黄色一级片免费的| 午夜精品国产一区二区电影| 欧美午夜高清在线| 欧美xxⅹ黑人| 日韩制服骚丝袜av| 一级毛片电影观看| 夜夜骑夜夜射夜夜干| 久久久久久久大尺度免费视频| 精品人妻在线不人妻| 欧美激情高清一区二区三区| 国产成人av教育| 9色porny在线观看| 亚洲avbb在线观看| 亚洲 欧美一区二区三区| 亚洲欧美日韩另类电影网站| 日韩视频在线欧美| 叶爱在线成人免费视频播放| 亚洲欧美成人综合另类久久久| 我的亚洲天堂| 国产精品99久久99久久久不卡| 乱人伦中国视频| 在线天堂中文资源库| 成人手机av| 午夜福利视频精品| 国产成人系列免费观看| 电影成人av| 人人妻人人添人人爽欧美一区卜| 丝瓜视频免费看黄片| 亚洲国产精品一区二区三区在线| 国产片内射在线| 美女视频免费永久观看网站| 国产一区二区三区在线臀色熟女 | 精品国产国语对白av| 精品国产一区二区三区四区第35| 亚洲av欧美aⅴ国产| 成人18禁高潮啪啪吃奶动态图| 婷婷色av中文字幕| 亚洲,欧美精品.| 亚洲精品自拍成人| 国产精品 欧美亚洲| 国产精品久久久久成人av| 亚洲av成人不卡在线观看播放网 | 一本色道久久久久久精品综合| 免费一级毛片在线播放高清视频 | 狠狠精品人妻久久久久久综合| 国产成人精品在线电影| 高清av免费在线| 欧美av亚洲av综合av国产av| 老司机午夜福利在线观看视频 | av欧美777| 十分钟在线观看高清视频www| 成人国产一区最新在线观看| 99久久国产精品久久久| 伊人久久大香线蕉亚洲五| 免费看十八禁软件| 亚洲中文日韩欧美视频| 国产免费现黄频在线看| 国产高清videossex| 人人妻,人人澡人人爽秒播| 电影成人av| 黄色a级毛片大全视频| 天天影视国产精品| 精品国产一区二区久久| 亚洲三区欧美一区| 啪啪无遮挡十八禁网站| 一本久久精品| 国产欧美亚洲国产| 午夜激情久久久久久久| 精品视频人人做人人爽| 免费不卡黄色视频| 在线亚洲精品国产二区图片欧美| 嫁个100分男人电影在线观看| 国产成人a∨麻豆精品| 精品人妻在线不人妻| 日本撒尿小便嘘嘘汇集6| 国产视频一区二区在线看| 日韩欧美免费精品| 久久国产亚洲av麻豆专区| 在线永久观看黄色视频| 男人操女人黄网站| 91精品伊人久久大香线蕉| 欧美日韩国产mv在线观看视频| 老熟妇仑乱视频hdxx| 国产一区二区激情短视频 | 中文精品一卡2卡3卡4更新| 色婷婷久久久亚洲欧美| 在线天堂中文资源库| 欧美少妇被猛烈插入视频| 国产野战对白在线观看| 极品少妇高潮喷水抽搐| 中文字幕色久视频| 在线亚洲精品国产二区图片欧美| 亚洲精华国产精华精| 久久精品亚洲熟妇少妇任你| 五月开心婷婷网| 可以免费在线观看a视频的电影网站| 99热全是精品| 国产在线视频一区二区| 91麻豆av在线| 波多野结衣一区麻豆| 中文字幕最新亚洲高清| 在线观看人妻少妇| 国产国语露脸激情在线看| 激情视频va一区二区三区| 精品免费久久久久久久清纯 | 他把我摸到了高潮在线观看 | 日本wwww免费看| 美女国产高潮福利片在线看| 国产欧美日韩精品亚洲av| 国产一区二区在线观看av| 一二三四在线观看免费中文在| 久久亚洲国产成人精品v| 成年人黄色毛片网站| 高清在线国产一区| 精品国产一区二区三区四区第35| 极品人妻少妇av视频| 国产成人精品久久二区二区免费| 精品亚洲成a人片在线观看| 日韩熟女老妇一区二区性免费视频| av又黄又爽大尺度在线免费看| 亚洲欧美精品综合一区二区三区| 男人添女人高潮全过程视频| 久久久久久久久久久久大奶| 搡老乐熟女国产| 午夜视频精品福利| 老司机影院毛片| 大片电影免费在线观看免费| 午夜视频精品福利| 成年人黄色毛片网站| 一级a爱视频在线免费观看| 成年美女黄网站色视频大全免费| 黑人操中国人逼视频| 亚洲久久久国产精品| 久久久精品94久久精品| 精品人妻1区二区| 亚洲一码二码三码区别大吗| 午夜久久久在线观看| 亚洲一区中文字幕在线| 久久人人97超碰香蕉20202| 欧美日韩国产mv在线观看视频| 亚洲精品国产av蜜桃| 亚洲国产欧美一区二区综合| 久久影院123| 黄色 视频免费看| 免费高清在线观看日韩| 久久久国产一区二区| 丁香六月天网| 两人在一起打扑克的视频| 999久久久国产精品视频| 欧美日韩视频精品一区| 高清欧美精品videossex| 日韩电影二区| 亚洲色图 男人天堂 中文字幕| 热re99久久国产66热| 国产精品久久久av美女十八| 中文欧美无线码| 一本久久精品| 日韩欧美一区二区三区在线观看 | 亚洲国产精品999| 国产精品久久久久久人妻精品电影 | 欧美 亚洲 国产 日韩一| 国产一区有黄有色的免费视频| 18在线观看网站| 制服人妻中文乱码| 亚洲精品中文字幕在线视频| 国产男女超爽视频在线观看| 人妻久久中文字幕网| 999久久久精品免费观看国产| 最近中文字幕2019免费版| 国产亚洲午夜精品一区二区久久| 午夜两性在线视频| 美女高潮到喷水免费观看| 亚洲第一欧美日韩一区二区三区 | 亚洲国产精品一区二区三区在线| 精品国产一区二区久久| 国产精品一区二区精品视频观看| 日韩,欧美,国产一区二区三区| 人人妻人人澡人人看| 老司机亚洲免费影院| 亚洲成av片中文字幕在线观看| 色婷婷av一区二区三区视频| 一级a爱视频在线免费观看| 黄色视频不卡| 国产成人免费观看mmmm| 久热爱精品视频在线9| 精品一区二区三卡| 免费观看a级毛片全部| 黑人操中国人逼视频| 日韩欧美一区二区三区在线观看 | 久热这里只有精品99| 国产成人av教育| 久久人人爽人人片av| 黄色视频不卡| 一区在线观看完整版| 曰老女人黄片| 亚洲精品美女久久av网站| 欧美一级毛片孕妇| 涩涩av久久男人的天堂| 久久久久精品国产欧美久久久 | 亚洲精品美女久久av网站| 青春草亚洲视频在线观看| 久久毛片免费看一区二区三区| 欧美中文综合在线视频| tube8黄色片| 日本wwww免费看| 老汉色∧v一级毛片| 午夜福利在线观看吧| 不卡av一区二区三区| 国产精品av久久久久免费| 91精品国产国语对白视频| 高清视频免费观看一区二区| 国产精品一区二区在线不卡| 午夜免费观看性视频| 亚洲第一欧美日韩一区二区三区 | 国产精品久久久久久精品古装| a在线观看视频网站| 亚洲五月色婷婷综合| 亚洲精华国产精华精| 国产精品久久久久久精品电影小说| 亚洲九九香蕉| 99九九在线精品视频| 国产在线免费精品| 人人澡人人妻人| kizo精华| 久久ye,这里只有精品| av视频免费观看在线观看| 欧美另类一区| 男人添女人高潮全过程视频| 老司机午夜福利在线观看视频 | 美女国产高潮福利片在线看| 在线看a的网站| 久久精品亚洲熟妇少妇任你| 亚洲欧美一区二区三区黑人| 国产亚洲精品一区二区www | 国产一级毛片在线| 中文字幕色久视频| 欧美黑人欧美精品刺激| 国产成人精品久久二区二区免费| av有码第一页| 免费一级毛片在线播放高清视频 | 一二三四社区在线视频社区8| 国产在视频线精品| 亚洲五月婷婷丁香| 精品人妻一区二区三区麻豆| 久久久久精品国产欧美久久久 | 高潮久久久久久久久久久不卡| 日韩欧美免费精品| 激情视频va一区二区三区| 亚洲av男天堂| a级片在线免费高清观看视频| 美女福利国产在线| 最近最新免费中文字幕在线| 亚洲av成人不卡在线观看播放网 | 亚洲中文av在线| 9热在线视频观看99| 在线观看免费午夜福利视频| 精品一区二区三卡| 日日爽夜夜爽网站| 亚洲国产欧美一区二区综合| 亚洲专区中文字幕在线| 国产精品欧美亚洲77777| 欧美日韩中文字幕国产精品一区二区三区 | 精品一区二区三区av网在线观看 | 女人高潮潮喷娇喘18禁视频| 亚洲免费av在线视频| 91成人精品电影| 精品欧美一区二区三区在线| 男女无遮挡免费网站观看| 午夜成年电影在线免费观看| 80岁老熟妇乱子伦牲交| 中国国产av一级| 久久毛片免费看一区二区三区| 青春草亚洲视频在线观看| 国产不卡av网站在线观看| 青春草亚洲视频在线观看| 在线精品无人区一区二区三| 日日夜夜操网爽| www.av在线官网国产| 亚洲av片天天在线观看| 黄色视频,在线免费观看| 欧美日韩国产mv在线观看视频| 中文字幕人妻丝袜制服| 老司机亚洲免费影院| 高清黄色对白视频在线免费看| 欧美日韩av久久| 汤姆久久久久久久影院中文字幕| 黄色 视频免费看| 三级毛片av免费| 午夜福利视频在线观看免费| 午夜激情久久久久久久| 欧美午夜高清在线| 精品国产乱子伦一区二区三区 | 一级毛片精品| 亚洲人成77777在线视频| 老司机深夜福利视频在线观看 | www日本在线高清视频| 男女边摸边吃奶| 建设人人有责人人尽责人人享有的| 国产精品 国内视频| 国产精品国产三级国产专区5o| 欧美日韩视频精品一区| 成人国产一区最新在线观看| 一本—道久久a久久精品蜜桃钙片| 国产一区有黄有色的免费视频| 亚洲国产中文字幕在线视频| 宅男免费午夜| 丁香六月欧美| 亚洲九九香蕉| 精品福利观看| 欧美亚洲 丝袜 人妻 在线| 每晚都被弄得嗷嗷叫到高潮| 老司机午夜十八禁免费视频| 日本91视频免费播放| 亚洲国产中文字幕在线视频| 高清视频免费观看一区二区| 91精品国产国语对白视频| 香蕉丝袜av| 亚洲欧美精品综合一区二区三区| 久久性视频一级片| 韩国精品一区二区三区| 大码成人一级视频| 久久久欧美国产精品| av国产精品久久久久影院| 俄罗斯特黄特色一大片| 欧美日韩av久久| 国产av又大| 十八禁高潮呻吟视频| 美女扒开内裤让男人捅视频| 久久中文看片网| 午夜福利视频精品| 狠狠狠狠99中文字幕| 亚洲欧美激情在线| 国产国语露脸激情在线看| 亚洲第一av免费看| 久久ye,这里只有精品| 黄色视频不卡| 老司机深夜福利视频在线观看 | 国产精品国产三级国产专区5o| 国产欧美日韩综合在线一区二区| tube8黄色片| 青春草亚洲视频在线观看|