# Use of ATR-FTIR Spectroscopy and Chemometrics for the Variation of Active Components in Different Harvesting Periods ofLonicera japonica

Authors: Yun-Chu Zhang; Jing Deng; Xiu-Lian Lin; Ya-Mei Li; Han-Xuan Sheng; Bo-Hou Xia; Li-Mei Lin
Journal: International Journal of Analytical Chemistry
 (2022)

---

## Abstract

Lonicera japonica Thunb is a commonly used Chinese herbal medicine, which belongs to the family Caprifoliaceae. The active components varied greatly during bud development. Research on the variation of the main active components is significant for the timely harvesting and quality control of Lonicera japonica. In this study, the attenuated total reflection Fourier transform infrared spectroscopy (ATR-FTIR) combined with the chemometric method was performed to investigate the variability of different harvesting periods of Lonicera japonica. The preliminary characterization from ATR-FTIR fingerprints showed various characteristic absorption peaks of the main active components from the different harvesting times, such as flavonoids, organic acids, iridoids, and volatile oils. Additionally, principal component analysis (PCA) scatter plots showed that there was a clear clustering trend in the samples of the same harvesting period, and the samples of the different harvesting periods could be well distinguished. Finally, further analysis by the orthogonal partial least-squares discriminant analysis (OPLS-DA) showed that there were regular changes in flavonoids, phenolic acids, iridoids, and volatile oils in different harvesting periods. Therefore, ATR-FTIR, as a novel and convenient analytical method, could be applied to evaluate the quality of Lonicera japonica.


---

## Body

## 1. Introduction

Traditional Chinese medicine has a long history in China. It is a summary of the long-term life experience of the Chinese people. With the development of medical undertakings, Chinese medicines play an increasingly important role in disease treatments, beauty care, regimen, and immunity enhancement. The composition of Chinese medicines is complex and susceptible to the origin, climate, and harvesting periods, making it challenging to control quality. The quality of Chinese medicines is determined mainly by the content of active components, which is closely related to the harvesting period [1–4]. Therefore, the harvesting period is an essential factor affecting the quality of Chinese medicines.Lonicera japonica is the dried flower buds, the initial opening flower, of Lonicera japonica Thunb. It has the efficacy of heat clearing, detoxifying, and evacuating wind heat [5]. The chemical composition of Lonicera japonica is complicated. It contains bioactive components, such as chlorogenic acid and its analogs, flavonoids, iridoid glycosides, and triterpenoid saponins [6]. At present, Lonicera japonica has been widely used in the medicine and food industry. However, the quality of Lonicera japonica on the market is mixed due to its concentrated flowering period, short harvesting time, and the effect of varieties, processing methods, and transportation conditions. Moreover, flower buds of different harvesting periods are easily mixed and distinctly different in the content of the main bioactive components, making it difficult to exert the corresponding clinical effect when used as drugs. Therefore, it is of great significance to study the variation of active components of Lonicera japonica in growth stages for its quality control.Currently, in-depth research on the dissimilarity in different harvesting periods ofLonicera japonica is limited. The research on Lonicera japonica mainly focuses on the relevant aspects of the content changes in active components. Jiang et al. [7] studied the anti-inflammatory activity from different flowering phases of Lonicera japonica by high-performance liquid chromatography (HPLC). They found that the anti-inflammatory effects of Lonicera japonica gradually decrease as the flowering progressed. However, the preprocessing of this method is relatively complicated, and it acquires plenty of time and energy. Fu et al. [8] applied the HPLC to study the content changes in the primary antioxidants at different developmental stages of Lonicera japonica. The study revealed that chlorogenic acid and its derivatives have the strongest antioxidant capacity. Additionally, the highest content of chlorogenic acid was found before flowering. This helps to choose the harvesting period in which Lonicera japonica had the strongest antioxidant capacity. However, the shortcoming is that the flowering period is relatively short and not clear enough. Therefore, there is a need to develop a method and technique that is handy, rapid, economical, strong fingerprinting, and capable of fully characterizing the chemical composition of Lonicera japonica.In addition to the above research methods, the common methods for studying the quality difference in Chinese medicinal materials in recent years include ultraviolet-visible spectrophotometry (UV-Vis), Fourier transform infrared spectroscopy (FTIR), liquid chromatography-mass spectrometer (LC-MS), gas chromatography-mass spectrometer (GC-MS), and nuclear magnetic resonance spectroscopy (NMR). Among them, UV-Vis has low sensitivity to compounds with weak UV absorption. The application scope of GC-MS is relatively limited, and it is often used to determine volatile compounds. The sample pretreatment and detection processes of LC-MS and NMR are complicated, and thus, they can be time-consuming and costly. However, FTIR technology is an effective, economical, fast, and strong fingerprinting detection method. It can be used nondestructively and rapidly to obtain biochemical fingerprints that provide information about molecular structure and composition [9]. Moreover, in recent years, ATR-FTIR technology has been developed based on FTIR technology, which can be equipped with ATR accessories for crystal materials, such as zinc selenide, germanium, and diamond. Only a small number of samples are needed, so sampling is quick and straightforward, with no pretreatment required. This technology is used for qualitative and quantitative analysis [10, 11]. It is a “green analysis method,” which has been extensively used in the research field of medicine and food chemical composition.Therefore, in this study, the ATR-FTIR method was performed to establish the chemical information fingerprint ofLonicera japonica in five different harvesting periods. The sample data from these five periods were further analyzed and compared by applying chemometric methods, such as PCA and OPLS-DA. Moreover, to ensure the sustainable utilization of medicinal resources with high-quality and high-yield, the variation characteristics of the five periods were also discussed.


## 2. Materials and Methods

### 2.1. Instruments

The Fourier transform infrared spectrometer is equipped with ATR accessory (Nicolet iS5, Thermo Scientific, USA), 101 types of electrothermal blowing dry box (Beijing Yongguangming Medical Instrument Factory), DJ-10A pulverizer (Shanghai Longtuo Instrument Equipment Co., Ltd.), and 100 mesh standard sieve tray (Shangyu City, Zhejiang Province Yarn Screening Factory).


### 2.2. Plant Materials and Sample Collection

Lonicera japonica flower buds were collected from the medicine plantation of Hunan University of Chinese Medicine in April in Changsha, Hunan, China (112°54′E; 28°08′N). The well-developed Lonicera japonica samples were harvested at different time points. According to the color of the alabastrum, the flower bud differentiation process was divided into five periods including young alabastrum, green alabastrum, white alabastrum, silvery, and golden flower periods, named ya, ga, wa, sf, and gf, respectively (sample characteristics of each period are shown in Figure 1). The samples of three different positions were mixed in each period as one sample, and five samples were collected in each period.Figure 1
Sample characteristics ofLonicera japonica in five periods.


### 2.3. Sample Preparation

The samples were rinsed with distilled water, were dried in the shade, and then placed in an oven. After drying at a constant temperature of 60°C until achieving constant weight, the samples were crushed and screened through a mesh size of 100. All samples were stored in a dryer until the ATR-FTIR analysis.


### 2.4. ATR-FTIR Spectra Collection

To deduct the interference of the absorption of dihydrogen oxide and carbon dioxide, the ATR-FTIR spectra of the air, accumulating 36 scans per spectra, were performed as background with 8 cm−1 resolution. When concerning sample test, a small amount of Lonicera japonica samples were taken on the dry and clean zinc selenide crystal material for spectral scanning in the range of 4000–500 cm−1 under the room temperature. After collecting every one sample, the ATR crystal was washed with a cotton ball dipped in alcohol until it was non-contaminated, and then, the next sample was collected.


### 2.5. Data Processing

Firstly, the ATR-FTIR fingerprint ofLonicera japonica was processed by Origin 2017 software, and then, SIMCA-P+ 15.0 software was adopted to conduct the first derivative (FD), second derivative (SD), standard normal variate (SNV), multiplicative scatter correction (MSC), Savitzky–Golay smoothing (S-G), exponentially weighted moving average (EWMA), and row center (RC) preprocessing on the spectral data. Each preprocessing method was processed by six scaling methods, namely unit variance (UV), unit variance none (UVN), Pareto (Par), Pareto none (ParN), centering (Ctr), and freeze (their definitions are shown in Table 1 [12, 13]). Then, the classification and discrimination models are constructed by combining PCA, OPLS-DA, and other chemometric methods. These preprocessing methods could improve the stability and accuracy of the model [14]. The optimal preprocessing method was SD. The R2 and Q2 values of the models range between 0 and 1, with higher values (greater than 0.5) indicating higher model fit and predictability [15]. Permutation tests (100 times) and S-line plots were used in the OPLS-DA model for model evaluation and screening of differential components. The R2Y (close to 1) and Q2Y intercepts (less than zero) in the permutation plot were applied to evaluate whether the model was overfitted [16]. S-line plots were used to identify statistically significant and potentially biochemically significant metabolites, based on both contributions to the model (greater than 0.8) and their reliability [17].Table 1
Definition of six scaling methods.
Scaling methodDefinitionUVThe variablej is centered and scaled to “unit variance,” which means that the base weight is computed as the reciprocal of the standard deviation of variable j computed around the mean (sdj)UVNSimilar to UV, the variable is not centered, which means that the standard deviation is computed around zeroParBetween no scaling and UV scaling. The variablej is centered and scaled to Pareto variance, which means that the base weight is computed as 1/sqrt (sdj)ParNSimilar to par, however, the variable is not centeredCtrA column-wise centering to transform values varies around zeroFreezeThe scaling weight of the variable is frozen and will not be recomputed as observations in the work set change or the variable metric is modified


## 2.1. Instruments

The Fourier transform infrared spectrometer is equipped with ATR accessory (Nicolet iS5, Thermo Scientific, USA), 101 types of electrothermal blowing dry box (Beijing Yongguangming Medical Instrument Factory), DJ-10A pulverizer (Shanghai Longtuo Instrument Equipment Co., Ltd.), and 100 mesh standard sieve tray (Shangyu City, Zhejiang Province Yarn Screening Factory).


## 2.2. Plant Materials and Sample Collection

Lonicera japonica flower buds were collected from the medicine plantation of Hunan University of Chinese Medicine in April in Changsha, Hunan, China (112°54′E; 28°08′N). The well-developed Lonicera japonica samples were harvested at different time points. According to the color of the alabastrum, the flower bud differentiation process was divided into five periods including young alabastrum, green alabastrum, white alabastrum, silvery, and golden flower periods, named ya, ga, wa, sf, and gf, respectively (sample characteristics of each period are shown in Figure 1). The samples of three different positions were mixed in each period as one sample, and five samples were collected in each period.Figure 1
Sample characteristics ofLonicera japonica in five periods.


## 2.3. Sample Preparation

The samples were rinsed with distilled water, were dried in the shade, and then placed in an oven. After drying at a constant temperature of 60°C until achieving constant weight, the samples were crushed and screened through a mesh size of 100. All samples were stored in a dryer until the ATR-FTIR analysis.


## 2.4. ATR-FTIR Spectra Collection

To deduct the interference of the absorption of dihydrogen oxide and carbon dioxide, the ATR-FTIR spectra of the air, accumulating 36 scans per spectra, were performed as background with 8 cm−1 resolution. When concerning sample test, a small amount of Lonicera japonica samples were taken on the dry and clean zinc selenide crystal material for spectral scanning in the range of 4000–500 cm−1 under the room temperature. After collecting every one sample, the ATR crystal was washed with a cotton ball dipped in alcohol until it was non-contaminated, and then, the next sample was collected.


## 2.5. Data Processing

Firstly, the ATR-FTIR fingerprint ofLonicera japonica was processed by Origin 2017 software, and then, SIMCA-P+ 15.0 software was adopted to conduct the first derivative (FD), second derivative (SD), standard normal variate (SNV), multiplicative scatter correction (MSC), Savitzky–Golay smoothing (S-G), exponentially weighted moving average (EWMA), and row center (RC) preprocessing on the spectral data. Each preprocessing method was processed by six scaling methods, namely unit variance (UV), unit variance none (UVN), Pareto (Par), Pareto none (ParN), centering (Ctr), and freeze (their definitions are shown in Table 1 [12, 13]). Then, the classification and discrimination models are constructed by combining PCA, OPLS-DA, and other chemometric methods. These preprocessing methods could improve the stability and accuracy of the model [14]. The optimal preprocessing method was SD. The R2 and Q2 values of the models range between 0 and 1, with higher values (greater than 0.5) indicating higher model fit and predictability [15]. Permutation tests (100 times) and S-line plots were used in the OPLS-DA model for model evaluation and screening of differential components. The R2Y (close to 1) and Q2Y intercepts (less than zero) in the permutation plot were applied to evaluate whether the model was overfitted [16]. S-line plots were used to identify statistically significant and potentially biochemically significant metabolites, based on both contributions to the model (greater than 0.8) and their reliability [17].Table 1
Definition of six scaling methods.
Scaling methodDefinitionUVThe variablej is centered and scaled to “unit variance,” which means that the base weight is computed as the reciprocal of the standard deviation of variable j computed around the mean (sdj)UVNSimilar to UV, the variable is not centered, which means that the standard deviation is computed around zeroParBetween no scaling and UV scaling. The variablej is centered and scaled to Pareto variance, which means that the base weight is computed as 1/sqrt (sdj)ParNSimilar to par, however, the variable is not centeredCtrA column-wise centering to transform values varies around zeroFreezeThe scaling weight of the variable is frozen and will not be recomputed as observations in the work set change or the variable metric is modified


## 3. Results and Discussion

### 3.1. ATR-FTIR Fingerprint Analysis

Infrared fingerprints have a large amount of information and strong specificity, which can provide the information of main chemical constituents in medicinal materials. Different chemical components have their own infrared characteristic peaks [18]. The difference in absorption value and peak intensity can reflect the dissimilarity of the main chemical components of Lonicera japonica in different periods. The ATR-FTIR fingerprint of different harvesting periods and the composition and group information of different bands are shown in Figure 2 and Table 2, respectively. The results showed that there was a sharp absorption peak appeared near 1105 cm−1 in the ya stage, which could be distinguished from other periods, corresponding to the overlapped absorption peak of -OH bending vibration and C-O stretching vibration of flavonoids or organic acids. The ga stage has obvious absorption near 1255 cm−1 and 1378 cm−1, corresponding to the characteristic absorption of C-O of organic acid. There was a strong infrared absorption peak near 1055 cm−1 at the wa stage, corresponding to the overlapped absorption peak of -OH bending vibration and C-O stretching vibration of phenolic acid or flavonoids. This indicated that the content of phenolic acids and flavonoids is abundant in this period. The sf stage has obvious absorption near 1656 cm−1, corresponding to the absorption peak of C=C stretching vibration of iridoids. At the gf stage, a sharp absorption peak near 2956 appeared, which could be distinguished from other periods, corresponding to the absorption peak of –CH3 stretching vibration of volatile oils. Additionally, there was an obvious absorption near 1634 cm−1, corresponding to the absorption peak of aromatic ring skeleton vibration of flavonoids or phenolic acid. It can be seen from the above analysis that the dissimilarity of these principal components might be an important reason for the variability of Lonicera japonica in different harvesting periods. However, due to the complexity of infrared fingerprints, it is difficult to distinguish the subtle differences in many bands by manual observation and comparison only; hence chemometric methods are desired to establish classification and discriminant models for further analysis.Figure 2
ATR-FTIR fingerprint ofLonicera japonica in different harvesting periods.Table 2
Infrared spectrum band assignments.
Wavenumber (cm−1)VibrationSuggested biomolecular assignmentReference4000–3500O-H stretchingH2O[14]3800–2600-OH stretchingWater molecules[19]3615νs(OH) freeWater molecules[19]3490νas(OH), HBWater molecules[19]3420, 3400νs(OH), HBWater molecules[19]3400O-H bond vibrationSugar compound[20]3400Free N-H stretchingAmide[19]3392O-H stretching[21]3307-OH deformation mode[22]∼3300ν(NH), HBAmide A, proteins[19]3286–3284N-H stretchingAmide A of proteins[23]3286–3284O-H stretchingPolysaccharides[23]3286N-H stretchingAmide[24]3280νs(OH), HBWater molecules[19]3280H-O-H stretching[25]∼3100ν(NH), HBAmide B, proteins[19]3077νCHAlkane/alkene[24]3030–2800Lipids[23]3020–3010ν(=CH)Unsaturated fatty acids, cholesterol esters[19]3008Olefinic HC=CHUnsaturated lipids[23]3000–2800-CH3 and -CH2 groupsPhospholipids and fatty acids[19]2960–2955νas(CH3)Protein side chains, phospholipids, ceramides, fatty acids[19]2960CH3- stretching (asym)Alkane/alkene[24]2957Asymmetric CH3 stretching[25]2956CH3 antisymmetric stretchingLipids, protein side chains[23]2924–2915νas(CH2)Mainly lipids: phospholipids, ceramides, and fatty acids[19]2923, 2853CH2- stretching[21]2923CH2 antisymmetric stretchingLipids (mainly)[23]2920Asymmetric CH2 stretching[25]2900C-H bondsSugar compounds[20]2875–2872νs(CH3)Protein side chains, phospholipids, ceramides, fatty acids[19]2872CH3- stretching (sym)Alkane/alkene[24]2872CH3 symmetric stretchingProteins (mainly), lipids, carbohydrates, nucleic acids[23]2855–2847νs(CH2)Mainly lipids: phospholipids, ceramides, fatty acids[19]2854CH2 symmetric stretchingLipids (mainly)[26]2844CH2- stretching[22]2700–2330NH+ stretching and overtones or combination bands in Fermi resonanceTertiary amine hydrochloride salts[9]2500–2000Unsaturated hydrocarbons[27]2442–2208O-C-O stretchingCO2[14]2400–2250CO2[26]2360, 2340CO2[26]1900The 2nd overtone of the CO bondCaffeoylquinic acids[28]1800–1650Chlorogenic acid[29]1800–900Fingerprint regionAll molecules[23]1742Carbonyl C=O stretchingCholesterol esters[23]1740Triglycerides[30]1740–1720ν(C=O)Phospholipids, esters, glycerides[19]1734, 1627, 1522, 1440, 1410, 1367, 1315, 1255Chlorogenic acid, flavonoids[22]1734C=O stretchingPolyphenol[22]1720The asymmetric overtone of the C-H bondCaffeoylquinic acids[28]1700–1400Fatty acids, lipids, proteins[30]1690–1610ν(C=O)Amide I (∼70–80% C=O stretch)[19]1682, 1639, 1471, 1284, 1181, 1111, 1032, 982, 950, 822, 789Organic acids, flavonoids[22]1682C=O stretchingChlorogenic acid[22]1650C=O stretchingAmide I (proteins, lipids, and carbohydrates)[24]1650The first C-H overtone[31]1644Triterpenoid saponin[32]1641C=O stretchingß-sheet structureAmide I (protein)[23]1636C=C vibrations (aromatic ring skeletal)[22]1635The crystalline water of sugar compounds[20]1627, 1522C=C vibrations (aromatic ring skeletal)[22]1600–1450Secologanic acid, chlorogenic acid, galuteolin[33]1590Amino acids[30]1560–1500δ(NH), ν(CN)Amide II (∼40–60% N-H in-plane bend, ∼20–40% C-N stretch)[19]1545, 1455, 1450, 1410, 880, 875Carbonate[26]1538N-H bend, C-N stretch a helical structureAmide II (protein)[23]1516Tyrosine[30]1491C=C aromatic ring stretching vibrations[9]1473–1468δ(CH2), δas(CH3)Proteins, lipids[19]1462δ(CH2), δas(CH3)Proteins, lipids[19]1454–1451CH2 bendingLipids[23]1454δ(CH2), δas(CH3)Proteins, lipids[19]1453CH2 scissoring[25]1440, 1410, 1376OH vibrational modesOrganic acid[22]1428, 1407NCH3 bending[9]1420–1300The second overtone of the carbonyl groupSecologanic acid, chlorogenic acid, galuteolin[33]1413C-H deformation vibration (CH3 and -CH2- groups)[28]1412O-C=O symmetric stretchingGlycine[25]1400ν(C=O) + δs(CH3)COO-: proteins, lipids, fatty acids[19]1398CH3 bendingProteins[14]1398COO- stretching (sym)Fatty acids, amino acids[14]1397COO- symmetric stretchingFatty acids[23]1394C=O stretch of COO-[25]1393, 1358νCO/νCOO-, δCHx (x = 1, 2, 3), δNH vibrational modesProteins and lipids[24]1350The second overtone of stretching O-H and C-H bonds[31]1343δ(CH2)C-H wagging[19]1333C-H2 waggingGlycine[25]1315, 1250C-O stretching[22]1243, 1240ν(CN), δ(NH) + νas(PO2−)Amide III (∼40% C-N stretch, 30% N-H in-plane bend, 20% methyl-C stretch)[19]1242AsymmetricPO2− stretch[25]1241νas(PO2−）Nucleic acids, phospholipids[30]1234PO2− antisymmetric stretchingNucleic acids (mainly), phospholipids[23]1200–1000C-O stretchingPolysaccharides, glycosides[21]1200–1000Nucleic acids, hydrocarbons, phosphates[24]1200–900ν(CC), δ(COH)Carbohydrates[19]1190–976, 962, 660–520, 472Phosphate ions[26]1171Ester C-O asymmetric stretch[25]1170–1100The second overtone of C-H (-CH2)Secologanic acid, chlorogenic acid, galuteolin[33]1167C-H stretching (CH3 and -CH2- groups)[28]1166CO-O-C antisymmetric stretchingEster bonds in cholesteryl esters, ribose ring formations: RNA[23]1155C-O stretchingOligosaccharides, triacylglycerols[14]1152, 1080, 1022δ(CHO), ν(C-O)Glycogen, glucose[24]1150–850Phosphate group[26]1150, 1020-OH deformation modes[22]1141, 1079C-O glycosidic bonds[20]1124ν(C-O)RNA[30]1120C-O bondsRibose[30]1118ν(C-O)RNA, carbohydrates[30]1117N-H3 rockingGlycine[25]1091p-Substituted aromatic vibrations[9]1090–1075νs(PO2−)Phospholipids[19]1080C-O stretch[25]1074PO2− symmetric stretchingNucleic acids, phospholipids, glycogen, polysaccharides, glycolipids[23]1051C-O stretchingStarch[14]1039C-N stretchingGlycine[25]1034Glycolytic components and nucleic acids[24]1010 (1225–950)Aromatic C-H in-plane bending[9]996Nucleotides[19]971C-N+-C stretchingRNA[23]950–400H2O[24]919Ribose ring vibrationsRNA/DNA[23]914–600O-C-O bendingCO2[14]1914CH2 rockingGlycine[25]894Theß-configuration and α-configuration[20]877, 744C-H deformation in proteins[19]833, 822Aromatic C-H out-of-plane bending[9]771C-O-C symmetry vibrationD-Glucopyranosyl ring[20]730–718ρ(CH2)Lipids[19]


### 3.2. Overall Variance Analysis by Principal Component Analysis

SIMCA is a supervised classification technique that uses samples with known origin (training samples) to derive a classification rule, which allows classifying new samples (test samples) with unknown origin in one of the classes, based on the values of the features of the new samples [34]. PCA is a mathematical algorithm that reduces the dimensionality of the data and can intuitively reflect the spatial distribution of the sample in the mathematical model while retaining most of the variation in the data set, which helps to understand the overall situation of the data and improve the accuracy of the model [35]. The results of the six scaling methods are shown in Figure 3. UV was chosen as the optimized scaling method for PCA model, and the results are shown in Figure 4 (three-dimensional). The model fitting parameter R2X = 0.784, and the model prediction parameter Q2 = 0.511, both of which are greater than 0.5, indicating that the model has a good fitting degree and strong predictive ability. There was a significant clustering trend of the Lonicera japonica samples in the same harvesting period. The samples in the ya stage and the gf stage could be significantly distinguished from those in other stages, indicating that the chemical components in these two periods are distinct. The spatial distribution of some samples in the ga stage, wa stage, and sf stage was relatively close, indicating that the variation of chemical components in these three periods was comparatively slight.Figure 3
PCA scatter plot of six scaling methods.
(a)(b)(c)(d)(e)(f)Figure 4
PCA scatter plot with UV scaling.


### 3.3. Individual Variance Analysis by Orthogonal Partial Least-Squares Discriminant Analysis

To better distinguish the differences between groups and to specify the dissimilarity between the two adjacent harvesting periods and their differential chemical compositions, the OPLS-DA was performed in this study [36]. The OPLS-DA scatter plot from the ya stage to the ga stage is shown in Figure 5(a). The model fitting parameter R2Y = l, and the model prediction parameter Q2 = 0.96, indicating that the model had a good fitting degree and strong predictive ability. It can be observed that the ya stage and the ga stage could be well distinguished, indicating that there were significant differences between the samples in these two periods. In the permutation test plot (Figure 5(b)), R2 = 0.999, Q2 = −0.0406, and all the points on the left of R2 and Q2 were lower than the rightmost point. It was concluded that the data were not overfitting, and the model could be used for the screening of differential markers. The contribution values greater than 0.8 in the S-line plot (Figure 6(a)) were 2987 cm−1, 2948 cm−1, 2900 cm−1, 1400 cm−1, 1254 cm−1, 1094 cm−1, 1078 cm−1, and 1051 cm−1. According to Table 1, it could be concluded that the main components that caused the change from the ya stage to the ga stage might be flavonoids or phenolic acid compounds, and nutrients, such as lipids, sugars, nucleic acids, and proteins. Research by Cui et al. [37] showed that the flavonoid content increased first and then decreased during the growth and development process of Lonicera japonica using the HPLC method. The increase in flavonoid content was probably caused by the increased activity of the chalcone isomerase, which catalyzes its synthesis [38]. Moreover, Zhang et al. [39] found that the content of phenolic acid compounds increased during the period from the ya stage to the ga stage. These research results were consistent with this study.Figure 5
Scatter plot and permutation test plot of OPLS-DA models.
(a)(b)(c)(d)(e)(f)(g)(h)Figure 6
S-line plot of OPLS-DA models.
(a)(b)(c)(d)The OPLS-DA scatter plot from the ga stage to the wa stage is shown in Figure5(c). The model fitting parameter R2Y = 1, and the model prediction parameter Q2 = 0.871, indicating that there was a good fitting degree and predictability. In the permutation test plot (Figure 5(d)), R2 = 0.988 and Q2 = −0.0889, and when the proportion of the displacing Y variable increases, the Q2 of the stochastic model gradually decreases, indicating that the data were not overfitting. According to the S-line plot (Figure 6(b)) based on the rule of the contribution values greater than 0.8, 3025 cm−1, 2979 cm−1, 2947 cm−1, 1710 cm−1, 1688 cm−1, and 1561 cm−1, 1388 cm−1, 1021 cm−1, and 806 cm−1 were confirmed. Accordingly, it could be concluded that the main components causing the change from the ga stage to the wa stage were similar to the result from the ya stage to the ga stage. However, among them, 1710 cm−1, 1688 cm−1, 1388 cm−1, and 806 cm−1 were the characteristic absorption of caffeoylquinic acid, chlorogenic acid, and flavonoids, which showed more chemical information of phenolic acids and flavonoids than other periods. Furthermore, the phenolic acids and flavonoids showed the most substantial infrared absorption peaks in the wa stage, indicating that the content of phenolic acids and flavonoids (the quality control index of Lonicera japonica) is at the highest level [40]. Meanwhile, Wang et al. [41] studied the content changes in the main active constituents of Lonicera japonica at different stages. The results showed that with the development of flower buds, the content of phenolic acids increased firstly, then decreased, and reached the peak value at the wa stage. Kong et al. [38] discovered that the content of flavonoids increased firstly and then decreased during the growth and development periods of Lonicera japonica. The highest content was also found at the wa stage. From the above study, phenolic acids and flavonoids were the differential components that caused the changes in the ga and wa stages, which were in line with the results of previous researchers.The OPLS-DA scatter plot from the wa stage to the sf stage is shown in Figure5(e). The model fitting parameter R2Y = 1, and the model prediction parameter Q2 = 0.919, indicating an excellent model fit and high predictability. In the permutation test plot (Figure 5(f)), R2 = 0.997 and Q2 = −0.0278, indicating that the model could be used for the screening of differential markers. The contribution values greater than 0.8 in the S-line plot (Figure 6(c)) were 2987 cm−1, 2873 cm−1, 1650 cm−1, 1376 cm−1, 1060 cm−1, and 1039 cm−1. It could be concluded that the main components causing the change from the wa stage to the sf stage might be iridoids, phenolic acids or flavonoids, and nutrients, such as lipids, sugars, nucleic acids, and proteins. Wang et al. [42] found that with the development of Lonicera japonica flower buds, the iridoids showed a trend, which decreased first, then increased, and finally dropped, with the highest content appearing at the wa stage and the sf stage. Fu et al. [8] observed that phenolic acid content was significantly reduced as the flowers opened. Therefore, phenolic acid content showed a downward trend from the wa stage to the sf stage. Cui et al. [37] revealed that the flavonoid content tends to decrease during the period from the wa stage to the sf stage. This might be attributed to the competition of chalcone isomerase for precursors at the sf stage [38]. Therefore, these results in this study were in agreement with previous reports.Finally, the OPLS-DA scatter plot from the sf stage to the gf stage is shown in Figure5(g). Similar to the above-established model, this model also showed a good fitting degree and strong predictive ability (R2Y = 0.999, Q2 = 0.827) and the permutation test plot (Figure 5(h)) showed nonsignificant overfitting (R2 = 0.993, Q2 = −0.187). According to the S-line plot (Figure 6(d)), the contribution values of 3022 cm−1, 2987 cm−1, 1667 cm−1, 1516 cm−1, 1392 cm−1, 1066 cm−1, and 1059 cm−1 were all greater than 0.8. This showed that the main components that caused the change from the sf stage to the gf stage might be volatile oils, phenolic acids or flavonoids, and nutrients, such as lipids, sugars, nucleic acids, and proteins. Wang et al. [41] showed that the content of volatile oil compounds increases in the flower bud stage, peaks in the sf stage, and then decreased. The change in volatile oil content showed certain flowering time dependence, indicating that some chemical components might accumulate during a particular period as the plant adapts to changes in its environment [43, 44]. Kong et al. [38] discovered that after the flowering of Lonicera japonica, the content of phenolic compounds significantly reduces till the gf stage when the content is minimal. The reduction in chlorogenic acid content might result from the increased biosynthesis of lignin using chlorogenic acid as a substrate during the gf stage. Zhang et al. [39] showed that the content of flavonoids decreases during the period from the sf stage to the gf stage of Lonicera japonica. These previous reports indicated that volatile oils, phenolic acids, and flavonoids were the differential components that cause changes in the sf and gf stages, which were in accordance with the results of this study.In summary, ATR-FTIR combined with chemometrics was used to establish chemical information fingerprint and classification and discrimination model ofLonicera japonica in different harvesting periods. It showed that this method could scientifically evaluate the variation of active components among different harvesting periods of Lonicera japonica. Moreover, ATR-FTIR technology is simple, rapid, effective, and applicable to Chinese medicine. However, in the determination of the infrared spectrum, there exist some spectral absorptions of the interfering components and the absorption peak probably generates superposition phenomenon, which has a certain impact on the accuracy of the determination results. Therefore, the determination of the differential components of Lonicera japonica at different harvesting periods still needs to be further studied in conjunction with other analytical techniques.


## 3.1. ATR-FTIR Fingerprint Analysis

Infrared fingerprints have a large amount of information and strong specificity, which can provide the information of main chemical constituents in medicinal materials. Different chemical components have their own infrared characteristic peaks [18]. The difference in absorption value and peak intensity can reflect the dissimilarity of the main chemical components of Lonicera japonica in different periods. The ATR-FTIR fingerprint of different harvesting periods and the composition and group information of different bands are shown in Figure 2 and Table 2, respectively. The results showed that there was a sharp absorption peak appeared near 1105 cm−1 in the ya stage, which could be distinguished from other periods, corresponding to the overlapped absorption peak of -OH bending vibration and C-O stretching vibration of flavonoids or organic acids. The ga stage has obvious absorption near 1255 cm−1 and 1378 cm−1, corresponding to the characteristic absorption of C-O of organic acid. There was a strong infrared absorption peak near 1055 cm−1 at the wa stage, corresponding to the overlapped absorption peak of -OH bending vibration and C-O stretching vibration of phenolic acid or flavonoids. This indicated that the content of phenolic acids and flavonoids is abundant in this period. The sf stage has obvious absorption near 1656 cm−1, corresponding to the absorption peak of C=C stretching vibration of iridoids. At the gf stage, a sharp absorption peak near 2956 appeared, which could be distinguished from other periods, corresponding to the absorption peak of –CH3 stretching vibration of volatile oils. Additionally, there was an obvious absorption near 1634 cm−1, corresponding to the absorption peak of aromatic ring skeleton vibration of flavonoids or phenolic acid. It can be seen from the above analysis that the dissimilarity of these principal components might be an important reason for the variability of Lonicera japonica in different harvesting periods. However, due to the complexity of infrared fingerprints, it is difficult to distinguish the subtle differences in many bands by manual observation and comparison only; hence chemometric methods are desired to establish classification and discriminant models for further analysis.Figure 2
ATR-FTIR fingerprint ofLonicera japonica in different harvesting periods.Table 2
Infrared spectrum band assignments.
Wavenumber (cm−1)VibrationSuggested biomolecular assignmentReference4000–3500O-H stretchingH2O[14]3800–2600-OH stretchingWater molecules[19]3615νs(OH) freeWater molecules[19]3490νas(OH), HBWater molecules[19]3420, 3400νs(OH), HBWater molecules[19]3400O-H bond vibrationSugar compound[20]3400Free N-H stretchingAmide[19]3392O-H stretching[21]3307-OH deformation mode[22]∼3300ν(NH), HBAmide A, proteins[19]3286–3284N-H stretchingAmide A of proteins[23]3286–3284O-H stretchingPolysaccharides[23]3286N-H stretchingAmide[24]3280νs(OH), HBWater molecules[19]3280H-O-H stretching[25]∼3100ν(NH), HBAmide B, proteins[19]3077νCHAlkane/alkene[24]3030–2800Lipids[23]3020–3010ν(=CH)Unsaturated fatty acids, cholesterol esters[19]3008Olefinic HC=CHUnsaturated lipids[23]3000–2800-CH3 and -CH2 groupsPhospholipids and fatty acids[19]2960–2955νas(CH3)Protein side chains, phospholipids, ceramides, fatty acids[19]2960CH3- stretching (asym)Alkane/alkene[24]2957Asymmetric CH3 stretching[25]2956CH3 antisymmetric stretchingLipids, protein side chains[23]2924–2915νas(CH2)Mainly lipids: phospholipids, ceramides, and fatty acids[19]2923, 2853CH2- stretching[21]2923CH2 antisymmetric stretchingLipids (mainly)[23]2920Asymmetric CH2 stretching[25]2900C-H bondsSugar compounds[20]2875–2872νs(CH3)Protein side chains, phospholipids, ceramides, fatty acids[19]2872CH3- stretching (sym)Alkane/alkene[24]2872CH3 symmetric stretchingProteins (mainly), lipids, carbohydrates, nucleic acids[23]2855–2847νs(CH2)Mainly lipids: phospholipids, ceramides, fatty acids[19]2854CH2 symmetric stretchingLipids (mainly)[26]2844CH2- stretching[22]2700–2330NH+ stretching and overtones or combination bands in Fermi resonanceTertiary amine hydrochloride salts[9]2500–2000Unsaturated hydrocarbons[27]2442–2208O-C-O stretchingCO2[14]2400–2250CO2[26]2360, 2340CO2[26]1900The 2nd overtone of the CO bondCaffeoylquinic acids[28]1800–1650Chlorogenic acid[29]1800–900Fingerprint regionAll molecules[23]1742Carbonyl C=O stretchingCholesterol esters[23]1740Triglycerides[30]1740–1720ν(C=O)Phospholipids, esters, glycerides[19]1734, 1627, 1522, 1440, 1410, 1367, 1315, 1255Chlorogenic acid, flavonoids[22]1734C=O stretchingPolyphenol[22]1720The asymmetric overtone of the C-H bondCaffeoylquinic acids[28]1700–1400Fatty acids, lipids, proteins[30]1690–1610ν(C=O)Amide I (∼70–80% C=O stretch)[19]1682, 1639, 1471, 1284, 1181, 1111, 1032, 982, 950, 822, 789Organic acids, flavonoids[22]1682C=O stretchingChlorogenic acid[22]1650C=O stretchingAmide I (proteins, lipids, and carbohydrates)[24]1650The first C-H overtone[31]1644Triterpenoid saponin[32]1641C=O stretchingß-sheet structureAmide I (protein)[23]1636C=C vibrations (aromatic ring skeletal)[22]1635The crystalline water of sugar compounds[20]1627, 1522C=C vibrations (aromatic ring skeletal)[22]1600–1450Secologanic acid, chlorogenic acid, galuteolin[33]1590Amino acids[30]1560–1500δ(NH), ν(CN)Amide II (∼40–60% N-H in-plane bend, ∼20–40% C-N stretch)[19]1545, 1455, 1450, 1410, 880, 875Carbonate[26]1538N-H bend, C-N stretch a helical structureAmide II (protein)[23]1516Tyrosine[30]1491C=C aromatic ring stretching vibrations[9]1473–1468δ(CH2), δas(CH3)Proteins, lipids[19]1462δ(CH2), δas(CH3)Proteins, lipids[19]1454–1451CH2 bendingLipids[23]1454δ(CH2), δas(CH3)Proteins, lipids[19]1453CH2 scissoring[25]1440, 1410, 1376OH vibrational modesOrganic acid[22]1428, 1407NCH3 bending[9]1420–1300The second overtone of the carbonyl groupSecologanic acid, chlorogenic acid, galuteolin[33]1413C-H deformation vibration (CH3 and -CH2- groups)[28]1412O-C=O symmetric stretchingGlycine[25]1400ν(C=O) + δs(CH3)COO-: proteins, lipids, fatty acids[19]1398CH3 bendingProteins[14]1398COO- stretching (sym)Fatty acids, amino acids[14]1397COO- symmetric stretchingFatty acids[23]1394C=O stretch of COO-[25]1393, 1358νCO/νCOO-, δCHx (x = 1, 2, 3), δNH vibrational modesProteins and lipids[24]1350The second overtone of stretching O-H and C-H bonds[31]1343δ(CH2)C-H wagging[19]1333C-H2 waggingGlycine[25]1315, 1250C-O stretching[22]1243, 1240ν(CN), δ(NH) + νas(PO2−)Amide III (∼40% C-N stretch, 30% N-H in-plane bend, 20% methyl-C stretch)[19]1242AsymmetricPO2− stretch[25]1241νas(PO2−）Nucleic acids, phospholipids[30]1234PO2− antisymmetric stretchingNucleic acids (mainly), phospholipids[23]1200–1000C-O stretchingPolysaccharides, glycosides[21]1200–1000Nucleic acids, hydrocarbons, phosphates[24]1200–900ν(CC), δ(COH)Carbohydrates[19]1190–976, 962, 660–520, 472Phosphate ions[26]1171Ester C-O asymmetric stretch[25]1170–1100The second overtone of C-H (-CH2)Secologanic acid, chlorogenic acid, galuteolin[33]1167C-H stretching (CH3 and -CH2- groups)[28]1166CO-O-C antisymmetric stretchingEster bonds in cholesteryl esters, ribose ring formations: RNA[23]1155C-O stretchingOligosaccharides, triacylglycerols[14]1152, 1080, 1022δ(CHO), ν(C-O)Glycogen, glucose[24]1150–850Phosphate group[26]1150, 1020-OH deformation modes[22]1141, 1079C-O glycosidic bonds[20]1124ν(C-O)RNA[30]1120C-O bondsRibose[30]1118ν(C-O)RNA, carbohydrates[30]1117N-H3 rockingGlycine[25]1091p-Substituted aromatic vibrations[9]1090–1075νs(PO2−)Phospholipids[19]1080C-O stretch[25]1074PO2− symmetric stretchingNucleic acids, phospholipids, glycogen, polysaccharides, glycolipids[23]1051C-O stretchingStarch[14]1039C-N stretchingGlycine[25]1034Glycolytic components and nucleic acids[24]1010 (1225–950)Aromatic C-H in-plane bending[9]996Nucleotides[19]971C-N+-C stretchingRNA[23]950–400H2O[24]919Ribose ring vibrationsRNA/DNA[23]914–600O-C-O bendingCO2[14]1914CH2 rockingGlycine[25]894Theß-configuration and α-configuration[20]877, 744C-H deformation in proteins[19]833, 822Aromatic C-H out-of-plane bending[9]771C-O-C symmetry vibrationD-Glucopyranosyl ring[20]730–718ρ(CH2)Lipids[19]


## 3.2. Overall Variance Analysis by Principal Component Analysis

SIMCA is a supervised classification technique that uses samples with known origin (training samples) to derive a classification rule, which allows classifying new samples (test samples) with unknown origin in one of the classes, based on the values of the features of the new samples [34]. PCA is a mathematical algorithm that reduces the dimensionality of the data and can intuitively reflect the spatial distribution of the sample in the mathematical model while retaining most of the variation in the data set, which helps to understand the overall situation of the data and improve the accuracy of the model [35]. The results of the six scaling methods are shown in Figure 3. UV was chosen as the optimized scaling method for PCA model, and the results are shown in Figure 4 (three-dimensional). The model fitting parameter R2X = 0.784, and the model prediction parameter Q2 = 0.511, both of which are greater than 0.5, indicating that the model has a good fitting degree and strong predictive ability. There was a significant clustering trend of the Lonicera japonica samples in the same harvesting period. The samples in the ya stage and the gf stage could be significantly distinguished from those in other stages, indicating that the chemical components in these two periods are distinct. The spatial distribution of some samples in the ga stage, wa stage, and sf stage was relatively close, indicating that the variation of chemical components in these three periods was comparatively slight.Figure 3
PCA scatter plot of six scaling methods.
(a)(b)(c)(d)(e)(f)Figure 4
PCA scatter plot with UV scaling.


## 3.3. Individual Variance Analysis by Orthogonal Partial Least-Squares Discriminant Analysis

To better distinguish the differences between groups and to specify the dissimilarity between the two adjacent harvesting periods and their differential chemical compositions, the OPLS-DA was performed in this study [36]. The OPLS-DA scatter plot from the ya stage to the ga stage is shown in Figure 5(a). The model fitting parameter R2Y = l, and the model prediction parameter Q2 = 0.96, indicating that the model had a good fitting degree and strong predictive ability. It can be observed that the ya stage and the ga stage could be well distinguished, indicating that there were significant differences between the samples in these two periods. In the permutation test plot (Figure 5(b)), R2 = 0.999, Q2 = −0.0406, and all the points on the left of R2 and Q2 were lower than the rightmost point. It was concluded that the data were not overfitting, and the model could be used for the screening of differential markers. The contribution values greater than 0.8 in the S-line plot (Figure 6(a)) were 2987 cm−1, 2948 cm−1, 2900 cm−1, 1400 cm−1, 1254 cm−1, 1094 cm−1, 1078 cm−1, and 1051 cm−1. According to Table 1, it could be concluded that the main components that caused the change from the ya stage to the ga stage might be flavonoids or phenolic acid compounds, and nutrients, such as lipids, sugars, nucleic acids, and proteins. Research by Cui et al. [37] showed that the flavonoid content increased first and then decreased during the growth and development process of Lonicera japonica using the HPLC method. The increase in flavonoid content was probably caused by the increased activity of the chalcone isomerase, which catalyzes its synthesis [38]. Moreover, Zhang et al. [39] found that the content of phenolic acid compounds increased during the period from the ya stage to the ga stage. These research results were consistent with this study.Figure 5
Scatter plot and permutation test plot of OPLS-DA models.
(a)(b)(c)(d)(e)(f)(g)(h)Figure 6
S-line plot of OPLS-DA models.
(a)(b)(c)(d)The OPLS-DA scatter plot from the ga stage to the wa stage is shown in Figure5(c). The model fitting parameter R2Y = 1, and the model prediction parameter Q2 = 0.871, indicating that there was a good fitting degree and predictability. In the permutation test plot (Figure 5(d)), R2 = 0.988 and Q2 = −0.0889, and when the proportion of the displacing Y variable increases, the Q2 of the stochastic model gradually decreases, indicating that the data were not overfitting. According to the S-line plot (Figure 6(b)) based on the rule of the contribution values greater than 0.8, 3025 cm−1, 2979 cm−1, 2947 cm−1, 1710 cm−1, 1688 cm−1, and 1561 cm−1, 1388 cm−1, 1021 cm−1, and 806 cm−1 were confirmed. Accordingly, it could be concluded that the main components causing the change from the ga stage to the wa stage were similar to the result from the ya stage to the ga stage. However, among them, 1710 cm−1, 1688 cm−1, 1388 cm−1, and 806 cm−1 were the characteristic absorption of caffeoylquinic acid, chlorogenic acid, and flavonoids, which showed more chemical information of phenolic acids and flavonoids than other periods. Furthermore, the phenolic acids and flavonoids showed the most substantial infrared absorption peaks in the wa stage, indicating that the content of phenolic acids and flavonoids (the quality control index of Lonicera japonica) is at the highest level [40]. Meanwhile, Wang et al. [41] studied the content changes in the main active constituents of Lonicera japonica at different stages. The results showed that with the development of flower buds, the content of phenolic acids increased firstly, then decreased, and reached the peak value at the wa stage. Kong et al. [38] discovered that the content of flavonoids increased firstly and then decreased during the growth and development periods of Lonicera japonica. The highest content was also found at the wa stage. From the above study, phenolic acids and flavonoids were the differential components that caused the changes in the ga and wa stages, which were in line with the results of previous researchers.The OPLS-DA scatter plot from the wa stage to the sf stage is shown in Figure5(e). The model fitting parameter R2Y = 1, and the model prediction parameter Q2 = 0.919, indicating an excellent model fit and high predictability. In the permutation test plot (Figure 5(f)), R2 = 0.997 and Q2 = −0.0278, indicating that the model could be used for the screening of differential markers. The contribution values greater than 0.8 in the S-line plot (Figure 6(c)) were 2987 cm−1, 2873 cm−1, 1650 cm−1, 1376 cm−1, 1060 cm−1, and 1039 cm−1. It could be concluded that the main components causing the change from the wa stage to the sf stage might be iridoids, phenolic acids or flavonoids, and nutrients, such as lipids, sugars, nucleic acids, and proteins. Wang et al. [42] found that with the development of Lonicera japonica flower buds, the iridoids showed a trend, which decreased first, then increased, and finally dropped, with the highest content appearing at the wa stage and the sf stage. Fu et al. [8] observed that phenolic acid content was significantly reduced as the flowers opened. Therefore, phenolic acid content showed a downward trend from the wa stage to the sf stage. Cui et al. [37] revealed that the flavonoid content tends to decrease during the period from the wa stage to the sf stage. This might be attributed to the competition of chalcone isomerase for precursors at the sf stage [38]. Therefore, these results in this study were in agreement with previous reports.Finally, the OPLS-DA scatter plot from the sf stage to the gf stage is shown in Figure5(g). Similar to the above-established model, this model also showed a good fitting degree and strong predictive ability (R2Y = 0.999, Q2 = 0.827) and the permutation test plot (Figure 5(h)) showed nonsignificant overfitting (R2 = 0.993, Q2 = −0.187). According to the S-line plot (Figure 6(d)), the contribution values of 3022 cm−1, 2987 cm−1, 1667 cm−1, 1516 cm−1, 1392 cm−1, 1066 cm−1, and 1059 cm−1 were all greater than 0.8. This showed that the main components that caused the change from the sf stage to the gf stage might be volatile oils, phenolic acids or flavonoids, and nutrients, such as lipids, sugars, nucleic acids, and proteins. Wang et al. [41] showed that the content of volatile oil compounds increases in the flower bud stage, peaks in the sf stage, and then decreased. The change in volatile oil content showed certain flowering time dependence, indicating that some chemical components might accumulate during a particular period as the plant adapts to changes in its environment [43, 44]. Kong et al. [38] discovered that after the flowering of Lonicera japonica, the content of phenolic compounds significantly reduces till the gf stage when the content is minimal. The reduction in chlorogenic acid content might result from the increased biosynthesis of lignin using chlorogenic acid as a substrate during the gf stage. Zhang et al. [39] showed that the content of flavonoids decreases during the period from the sf stage to the gf stage of Lonicera japonica. These previous reports indicated that volatile oils, phenolic acids, and flavonoids were the differential components that cause changes in the sf and gf stages, which were in accordance with the results of this study.In summary, ATR-FTIR combined with chemometrics was used to establish chemical information fingerprint and classification and discrimination model ofLonicera japonica in different harvesting periods. It showed that this method could scientifically evaluate the variation of active components among different harvesting periods of Lonicera japonica. Moreover, ATR-FTIR technology is simple, rapid, effective, and applicable to Chinese medicine. However, in the determination of the infrared spectrum, there exist some spectral absorptions of the interfering components and the absorption peak probably generates superposition phenomenon, which has a certain impact on the accuracy of the determination results. Therefore, the determination of the differential components of Lonicera japonica at different harvesting periods still needs to be further studied in conjunction with other analytical techniques.



---

*Source: 8850914-2022-03-07.xml*