Simultaneous Quantiﬁcation of Diarylheptanoids and Phenolic Compounds in Juglans mandshurica Maxim. by UPLC–TQ-MS

: The immature epicarps of Juglans mandshurica and Juglans regia have been used as folk medicine for the treatment of cancer in China. Other parts of the J. mandshurica plant, including leaves, branches, barks, and stems, have reported antitumor activities. We previously found that various diarylheptanoids and phenolic compounds isolated from J. mandshurica epicarps show signiﬁcant antitumor activities. However, there are no reports of quantitative analysis of diarylheptanoids and phenolic compounds of J. mandshurica . In this study, a validated quantitative method, based on ultraperformance liquid chromatography coupled with triple-quadrupole tandem mass spectrometry, was employed to determine the contents of eight diarylheptanoids and seven phenolic compounds in the epicarps of J. mandshurica during different growth periods, in different parts of the plant, and in the epicarps of two Juglans species. The most successful J. mandshurica epicarp harvesting time fell between Day 12 and Day 27. The leaves of J. mandshurica showed potential for medical use as they had the highest content of the 15 compounds (3.399 ± 0.013 mg/g). We showed for the ﬁrst time that the total content of diarylheptanoids in J. mandshurica is higher than that in J. regia , though, conversely, J. regia has higher contents of phenolic compounds. The method developed in this study is practical and simple and can be applied for quantitative analysis for evaluating the intrinsic quality of J. mandshurica .

We recently reported on the structures and cytotoxic activities of various diarylheptanoids and phenolic compounds that are constituents of green walnut (J. mandshurica and J. regia) [18,19]. Diarylheptanoids are distinctive biological components responsible for

Reagents and Materials
The epicarps from different growth periods of J. mandshurica were collected every 2 weeks from Harbin Medical University, Harbin, Heilongjiang, China, from 18 June 2020 until the fruits fell. The samples from 18 June, 30 June, 15 July, 28 July, 11 August, and 30 August were marked as Day 0, Day 12, Day 27, Day 40, Day 54 and Day 73, respectively. The sample collected on 30 August 2020 was denoted as HLJ-1 in the comparison to samples of J. regia. The epicarps, branches, leaves, stems, and stem bark of J. mandshurica were picked from Harbin Medical University, Harbin, Heilongjiang, China on 30 August 2020. The epicarps of J. mandshurica, including JiLin-1, JiLin-2, and JiLin-3, were collected from the regions of Shulan, Changchun, Songyuan, respectively, in Jilin, China, on 30 August 2020. Additionally, the epicarps of J. regia, including Luguo1hao, Zhipihetao, Xinzaofeng, Zha343, Xinjufeng, Lipin1hao, and Xilin3hao, which are the most common local varieties, cultivated at a high rate in Xinjiang, were collected from Aksu, Xinjiang, China, on 1 September 2020. After collection, the fresh samples were dried at room temperature in the shade, and then the dry raw materials were stored in sealed plastic bag at room temperature before analysis. All the plant materials were identified by Professor Jinhui Wang (College of Pharmacy, Harbin Medical University, Harbin, China). Voucher specimens (202009XJLKY) were deposited in the Xinjiang Academy of Forestry.

Preparation of Standard and Sample Solutions
Fifteen reference substances were separately dissolved in methanol to prepare the individual standard stock solutions with a concentration of 1.0 mg/mL. Then, the 15 mixed standards were progressively diluted to 10 μg/mL. Finally, the 10 μg/mL mixed standard solution was serially diluted with methanol to obtain nine standard working solutions with different concentrations for the construction of calibration curves. All solutions were stored at −20 °C and filtered through 0.22 μm membrane filter before analysis.
Each sample was accurately weighed to 250 mg and mixed with 5 mL 70% MeOHwater (v/v). The samples were extracted ultrasonically for 30 min and cooled to room temperature; we compensated for weightlessness and then centrifuged the samples at 12,000 rpm for 5 min. Then, 1 mL of the supernatant was filtered through an HLB Cartridge Plus Short and eluted with 4 mL methanol in a 5 mL volumetric flask. Methanol was added to bring the volume to scale and mixed thoroughly. Three replicates were prepared for every sample. LC-MS-grade acetonitrile and methanol were purchased from Fisher Scientific Co. (Loughborough, UK). HPLC-grade methanol was obtained from Sigma-Aldrich (St. Louis, MO, USA). Deionized water (18.2 Ω) was further purified using a Milli-Q system (Millipore, Billerica, MA, USA). Oasis ® PRiME HLB Cartridge Plus Short (335 mg) was provided by Waters Technology (Milford, MA, USA).

Preparation of Standard and Sample Solutions
Fifteen reference substances were separately dissolved in methanol to prepare the individual standard stock solutions with a concentration of 1.0 mg/mL. Then, the 15 mixed standards were progressively diluted to 10 µg/mL. Finally, the 10 µg/mL mixed standard solution was serially diluted with methanol to obtain nine standard working solutions with different concentrations for the construction of calibration curves. All solutions were stored at −20 • C and filtered through 0.22 µm membrane filter before analysis.
Each sample was accurately weighed to 250 mg and mixed with 5 mL 70% MeOHwater (v/v). The samples were extracted ultrasonically for 30 min and cooled to room temperature; we compensated for weightlessness and then centrifuged the samples at 12,000 rpm for 5 min. Then, 1 mL of the supernatant was filtered through an HLB Cartridge Plus Short and eluted with 4 mL methanol in a 5 mL volumetric flask. Methanol was added to bring the volume to scale and mixed thoroughly. Three replicates were prepared for every sample.

Quantitative Analysis Conditions
The quantitation method was performed using an Acquity UPLC I-class system (Waters, Milford, MA, USA) with an Acquity UPLC BEH C18 column (100 mm × 2.1 mm, 1.7 µm, Waters, Milford, MA, USA) used at 40 • C. A binary solvent system was used, consisting of solvent A (water containing 0.1% formic acid) and solvent B (acetonitrile). The A Xevo TQ (Waters MS Technologies, Manchester, UK), equipped with an electrospray ionization ion (ESI) source, was applied for mass spectrometry data acquisition in negative ionization mode. The ion scan mode adopted was multiple reaction monitoring (MRM) mode and the MS experimental conditions were set as follows: electrospray capillary voltage 2.5 kV, cone voltage 40 V, source temperature 110 • C, desolvation temperature 450 • C, cone gas flow 50 L/h, and desolvation gas 800 L/h. Nitrogen and argon were used as the nebulizer and collision gas, respectively. The UPLC-TQ-MS system was controlled by the MassLynx Table S3).

Validation of the Quantitative Method
The calibration curve for each diarylheptanoid was constructed by plotting the peak area (y) against corresponding concentrations (x). The limit of detection (LOD) and limit of quantification (LOQ) were determined by serial dilution of the mixed standard solution at the lowest concentrations of each calibration curves at signal-to-noise ratios (SNRs) of 3 and 10.
The precision was determined by six repetitive injections. To determine the repeatability, the same samples were measured in six replicates. The sample solution was stored in sample manager and analyzed at 0, 2, 4, 8, 12, and 24 h to investigate the stability. The accuracy of the method was expressed according to the recovery. We placed the known mixed standard solutions into the same sample. The results were evaluated by the relative standard deviation (RSD) values, and six replicates were performed for the analysis.

Data Analysis
Statistical analyses of the quantitative data were conducted using Microsoft Office Excel 2016 (Microsoft, Redmond, WA, USA), and the quantitative contents of the representative components were imported into the SIMCA 14.1 software (Sartorius, Goettingen, Germany) to perform chemometric analysis, including construction of the principal component analysis (PCA) and orthogonal projections to latent structures discriminant analysis (OPLS-DA) models. R2X and Q2 (cumulative) are usually used to evaluate the quality of the PCA model; the closer both R2X and Q2 (cumulative) are to 1, the higher the quality of the PCA model. R2X, R2Y, and Q2 were comparably used in the OPLS-DA models. The closer the obtained values are to 1, the higher the model quality.

Method Validations
The chemical structures of the 15 reference substances are shown in Figure 1. The extracted ion chromatogram of the reference substances (ESI−) in J. mandshurica is depicted in Figure 2.

Method Validations
The chemical structures of the 15 reference substances are shown in Figure 1. The extracted ion chromatogram of the reference substances (ESI−) in J. mandshurica is depicted in Figure 2. According to the peak area of each standard solution (y) and the concentration of each analyte (x), a series of the linear calibration curves of the reference standard solutions with seven concentrations was plotted. The correlation coefficient (R 2 ) of linearity was greater than 0.9990 for all investigated compounds. The method showed a good precision in the range of 0.47-4.11%. The recoveries ranged from 99.93% to 108.81% with an RSD% According to the peak area of each standard solution (y) and the concentration of each analyte (x), a series of the linear calibration curves of the reference standard solutions with seven concentrations was plotted. The correlation coefficient (R 2 ) of linearity was greater than 0.9990 for all investigated compounds. The method showed a good precision in the range of 0.47-4.11%. The recoveries ranged from 99.93% to 108.81% with an RSD% less than 3.38%. The RSD% values of the stability and repeatability were less than 2.72% and 2.28%, respectively. The calibration curves, correlation coefficient, linear range, LOD, and LOQ of each compound are listed in Table 1, and they confirm the feasibility of the method for the quantitative assessment of J. mandshurica. * LOD-limit of detection (S/N = 3); LOQ-limit of quantification (S/N = 10).

Quantification of the 15 Compounds in the Epicarps of J. mandshurica
In this study, six J. mandshurica epicarp samples were collected during the period of fruit growth to maturity (18 June-30 August, Day 0-Day 73) at two-week intervals. Then, eight diarylheptanoids and seven phenols in epicarps were quantitatively analyzed using UPLC-TQ-MS. The changing trends and compositions in compound contents for each sample are described below.
The detailed contents of the 15 compounds in J. mandshurica epicarps are listed in Supplementary Tables S1 and S2. By comparing the content of diarylheptanoids in the epicarps collected from June to August ( Figure 3A), we found that with the increase in growth duration, the total content of diarylheptanoids in the epicarps gradually decreased. We also found that the total content of diarylheptanoids in the epicarps of J. mandshurica in June (2.86 ± 0.007 mg/g on Day 0 and 2.761 ± 0 mg/g on Day 12, Table S1) was remarkably higher than in July and August. Furthermore, MM, 11-OL, and 3-OL were the main compounds in epicarps. According to the line graphs of the content of each compound in Figure 3B,C, on Day 12, we found that except for 11-OL and AK, the contents of the other diarylheptanoid compounds increased slightly, and the contents of all compounds dropped sharply by 15 July (Day 27) and levels were stable in August. Most of the time, the MM content was the highest (from 0.043 ± 0 to 1.011 ± 0.001 mg/g, Table S1), except for Day 0, when the content of 11-OL was the highest (1.165 ± 0.001 mg/g, Table S1), possibly due to conversion of 11-OL to MM.
Similarly, Figure 3D shows that the total content of phenolic compounds in the epicarps gradually decreased over time and increased slightly at the end of August (Day 73). Except for compounds SE and CA, all other compounds were the lowest on Day 54, and their levels increased to a certain extent by Day 73 (Figure 3E,F). Among these compounds, GA was the main constituent and had the highest content in J. mandshurica epicarps (from 0.044 ± 0 to 0.997 ± 0 mg/g, Table S2). In addition, the content of these phenolic compounds was relatively stable from Day 12 to Day 40.
To further explore the relationship between the contents of the studied compounds in J. mandshurica epicarps and various growth periods, PCA was conducted and the OPLS-DA models of the 15 compounds were constructed (Figure 4). In a PCA score plot, the closer the samples, the more similar they are in terms of their contents of the studied compounds [33]. On the contrary, a sample distant from others indicates it has a significantly different profile of compound contents. The R 2 X and Q 2 of the PCA models were 0.996 and 0.985, respectively, indicating the model fit the data well. In Figure 4A, samples from Day 0 and Day 12 are notably separated from others, which indicates that the epicarps from June were remarkably different from those of July and August. Additionally, OPLS-DA ( Figure 4B) was performed to explore the difference between the epicarps in the different growth periods based on hierarchical clustering analysis (HCA) ( Figure S1). R 2 X, R 2 Y, and Q 2 were 0.994, 1, and 1 for the OPLS-DA model, respectively, which revealed that the model was well-matched. In the OPLS-DA model, the latter four growth periods (Day 27-73) were significantly separated from the first two growth periods (Day 0-12) with respect to the x-axis. The scatter plot (S plot) and variable importance plot (VIP) from the OPLS-DA model was used to identify the most relevant variables that were distinguished among the first two growth periods and the latter four growth periods. The S plot is a graphical description of the covariance and the correlation structure between the x-loading variables and the predictive score. Therefore, variables situated in the upper-right and lower-left of the S plot (VIP > 1) are potential markers with both high model influence and high reliability [34]. By filtering the compounds with VIP > 1.0, we found that 11-OL, MM, 3-OL, and GA contributed more to distinguishing the different growth periods (Supplementary Figure S2A,B). These compounds are good candidates for further exploring how their profiling can be used to identify growth stages J. mandshurica epicarps. Similarly, Figure 3D shows that the total content of phenolic compounds in the epicarps gradually decreased over time and increased slightly at the end of August (Day 73). Except for compounds SE and CA, all other compounds were the lowest on Day 54, and their levels increased to a certain extent by Day 73 (Figure 3E,F). Among these compounds, GA was the main constituent and had the highest content in J. mandshurica epicarps (from 0.044 ± 0 to 0.997 ± 0 mg/g, Table S2). In addition, the content of these phenolic compounds was relatively stable from Day 12 to Day 40.
To further explore the relationship between the contents of the studied compounds in J. mandshurica epicarps and various growth periods, PCA was conducted and the OPLS-DA models of the 15 compounds were constructed (Figure 4). In a PCA score plot, the closer the samples, the more similar they are in terms of their contents of the studied compounds [33]. On the contrary, a sample distant from others indicates it has a significantly different profile of compound contents. The R 2 X and Q 2 of the PCA models were 0.996 and 0.985, respectively, indicating the model fit the data well. In Figure 4A, samples from Day 0 and Day 12 are notably separated from others, which indicates that the epicarps from June were remarkably different from those of July and August. Additionally, OPLS-DA ( Figure 4B) was performed to explore the difference between the epicarps in the different growth periods based on hierarchical clustering analysis (HCA) ( Figure S1). R 2 X, R 2 Y, and Q 2 were 0.994, 1, and 1 for the OPLS-DA model, respectively, which revealed that the model was well-matched. In the OPLS-DA model, the latter four growth periods MM, 3-OL, and GA contributed more to distinguishing the different growth periods (Supplementary Figure S2A,B). These compounds are good candidates for further exploring how their profiling can be used to identify growth stages J. mandshurica epicarps. From the above results, we conclude that the overall trend in the total content of diarylheptanoids and phenolic compounds indicates their content was higher before Day 27 and Day 40, respectively. However, the walnut fruits only start to form in June, so their total amount was insufficient [35]. As such, we preliminarily determined that the harvesting time of epicarps in Harbin is between Day 12 and Day 27. This study of diarylheptanoids and phenolic compounds lays a basic foundation for the establishment of a rea- From the above results, we conclude that the overall trend in the total content of diarylheptanoids and phenolic compounds indicates their content was higher before Day  27 and Day 40, respectively. However, the walnut fruits only start to form in June, so their total amount was insufficient [35]. As such, we preliminarily determined that the harvesting time of epicarps in Harbin is between Day 12 and Day 27. This study of diarylheptanoids and phenolic compounds lays a basic foundation for the establishment of a reasonable period for harvesting J. mandshurica epicarps.

Quantification of the 15 Compounds in the Different Parts of J. mandshurica
The detailed contents of the 15 compounds in the five plant parts of J. mandshurica are shown in Supplementary Tables S1 and S2. Comparing the total contents of diarylheptanoids in the five parts (epicarps, branches, leaves, stem barks, and stems) of J. mandshurica ( Figure 5A), we determined that the total content in leaves was the highest (0.557 ± 0.006 mg/g, Table S1). The trend from high to low content was leaves, branches, stem bark, stems, and epicarps. Unfortunately, JG, JB, MC, and AK, which are all diphenyl diarylheptanoids, were not detected in the stem bark or stems of J. mandshurica, whereas 3-OL, which is a chain diarylheptanoid, was highest in stem bark (0.314 ± 0.002 mg/g, Table S1) and stems (0.183 mg/g, Table S1). In addition, JG and MC were not detected in the branches or epicarps of J. mandshurica, respectively, but MM, which is a diphenyl ether diarylheptanoid, was highest in the epicarps (0.043 ± 0 mg/g, Table S1) and leaves (0.326 ± 0.003 mg/g, Table S1). Furthermore, as the compound percentage diagram in Figure 5B shows, the percentage compositions of epicarps and leaves were similar, and MM and 3-OL accounted for a higher proportion in the leaves, epicarps, and branches.
Comparison of the total content of phenolic compounds ( Figure 5C), the highest content was found in the leaves (2.842 ± 0.007 mg/g, Table S2), followed by branches, stem bark, epicarps, and stems. IQ, QN, and CA were not detected in the stems. However, IQ (1.868 ± 0.005 mg/g) and HE (0.758 ± 0.007 mg/g) were highest in the leaves and GA was highest in the stem bark (0.396 ± 0 mg/g). In terms of percentage content ( Figure 5D), GA accounted for a higher proportion in the stems and stem bark, IQ in the leaves, and HE in the branches and epicarps. To further explore the similarities and differences in compound profiles of the five parts of J. mandshurica, PCA and OPLS-DA models were built ( Figure 6). In the score plot of the PCA model (R 2 X and Q 2 were 1 and 0.999, respectively; Figure 6A), epicarps, branches, stem bark, and stems grouped into one category, with leaves being far from the other parts based on HCA ( Figure S3). The score plot of the OPLS-DA models (R 2 X, R 2 Y, and Q 2 were 0.853, 1, and 0.999, respectively; Figure 6B) shows a good separation between the leaves and other parts, which suggests the content in the leaves is different from the To further explore the similarities and differences in compound profiles of the five parts of J. mandshurica, PCA and OPLS-DA models were built ( Figure 6). In the score plot of the PCA model (R 2 X and Q 2 were 1 and 0.999, respectively; Figure 6A), epicarps, branches, stem bark, and stems grouped into one category, with leaves being far from the other parts based on HCA ( Figure S3). The score plot of the OPLS-DA models (R 2 X, R 2 Y, and Q 2 were 0.853, 1, and 0.999, respectively; Figure 6B) shows a good separation between the leaves and other parts, which suggests the content in the leaves is different from the others. According to the S plot and filtering the compounds with VIP > 1, IQ, HE, and MM were the main contribution to the differences between the leaves and other parts of J. mandshurica.
As such, we found that the contents of both diarylheptanoids and phenolic compounds are higher in leaves. As a renewable resource, leaves have considerable potential for utilization due to their variety of chemical components and biological activities [7]. Hence, the leaves of J. mandshurica should be further explored to determine their therapeutic potential and application in medicine. To further explore the similarities and differences in compound profiles of the five parts of J. mandshurica, PCA and OPLS-DA models were built ( Figure 6). In the score plot of the PCA model (R 2 X and Q 2 were 1 and 0.999, respectively; Figure 6A), epicarps, branches, stem bark, and stems grouped into one category, with leaves being far from the other parts based on HCA ( Figure S3). The score plot of the OPLS-DA models (R 2 X, R 2 Y, and Q 2 were 0.853, 1, and 0.999, respectively; Figure 6B) shows a good separation between the leaves and other parts, which suggests the content in the leaves is different from the others. According to the S plot and filtering the compounds with VIP > 1, IQ, HE, and MM were the main contribution to the differences between the leaves and other parts of J. mandshurica.
As such, we found that the contents of both diarylheptanoids and phenolic compounds are higher in leaves. As a renewable resource, leaves have considerable potential for utilization due to their variety of chemical components and biological activities [7]. Hence, the leaves of J. mandshurica should be further explored to determine their therapeutic potential and application in medicine.  To investigate the contents and composition of the 15 studied compounds in the epicarps of J. mandshurica (JMM) and J. regia (JRL), which belong to the same genus Juglans and have similar properties [25], we collected the samples of HLJ-1 and JiLin-1~3, belonging to J. mandshurica, and selected seven cultivars of J. regia, common walnuts widely cultivated in Xinjiang, including Luguo1hao, Zhipihetao, Xinzaofeng, Zha343, Xinjufeng, Lipin1ha, and Xilin3hao. All the samples were collected before the fruits had ripened and fallen.
The detailed contents of the 15 compounds in the epicarps of J. mandshurica and J. regia are provided in Supplementary Tables S1 and S2. Figure 7A shows that the total content of diarylheptanoids of JiLin-1 (1.005 ± 0.002 mg/g, Table S1) was the highest, followed closely by JiLin-2 (0.632 ± 0.004 mg/g) and JiLin-3 (0.304 ± 0 mg/g). However, despite being the same species, their total content of HLJ-1 (0.129 ± 0 mg/g) was much lower. We did not detect the presence of AK in any of the seven cultivars of J. regia and 3-ONE was not detected in six cultivars of J. regia, only in Xinjufeng. Interestingly, in terms of compound composition (Figure 7B), the HLJ-1 in J. mandshurica was similar to the samples of JiLin-1~3. Particularly, JB, MM, and 3-OL were the main components in the epicarps of J. mandshurica, whereas in the epicarps of J. regia, the main components were JB and MC. Figure 7C shows that the total content of phenolic compounds in the epicarps of the different local varieties of J. regia was generally higher than that of J. mandshurica; Xinjufeng had the highest content of phenolic compounds (3.662 ± 0.006 mg/g, Table S2). In the compound percentage graph ( Figure 7D), SE and GA are shown as the main constituents in the epicarps of J. regia, except for Zhipihetao, whereas HE and IQ accounted for a higher proportion in the epicarps of JiLin 1-3. HE, SE, and GA were the main contributing compounds in HLJ-1 and Zhipihetao. lowed closely by JiLin-2 (0.632 ± 0.004 mg/g) and JiLin-3 (0.304 ± 0 mg/g). However, de-spite being the same species, their total content of HLJ-1 (0.129 ± 0 mg/g) was much lower. We did not detect the presence of AK in any of the seven cultivars of J. regia and 3-ONE was not detected in six cultivars of J. regia, only in Xinjufeng. Interestingly, in terms of compound composition (Figure 7B), the HLJ-1 in J. mandshurica was similar to the samples of JiLin-1~3. Particularly, JB, MM, and 3-OL were the main components in the epicarps of J. mandshurica, whereas in the epicarps of J. regia, the main components were JB and MC.  Figure 7C shows that the total content of phenolic compounds in the epicarps of the different local varieties of J. regia was generally higher than that of J. mandshurica; Xinjufeng had the highest content of phenolic compounds (3.662 ± 0.006 mg/g, Table S2). In the compound percentage graph ( Figure 7D), SE and GA are shown as the main constituents in the epicarps of J. regia, except for Zhipihetao, whereas HE and IQ accounted To determine the similarities and differences of the compounds in the epicarps of the different varieties, PCA and OPLS-DA models ( Figure 8A,B) were built. The R 2 X and Q 2 of the PCA model of the two species were 0.991 and 0.910, respectively, and R 2 X, R 2 Y, and Q 2 were 0.891, 0.987, and 0.980, respectively, for J. mandshurica vs. J. regia, which showed that the models fit the data well. Based on HCA ( Figure S5), the Euclidean distance and Ward clustering algorithm were selected; the score plot of the PCA model of the two Juglans species revealed four groups ( Figure 8A). JL-1 and JL-2 (green plots) were separated into group 1; Xinzaofeng and Xinjufeng were in group 2 (blue plots); HLJ-1 and JL-3 constituted group 3 (red plots); and Zhipihetao, Lipin1hao, Zha343, Lugo1hao, and Xilin3hao were designated as group 4 (yellow plots). To further explore the differences between J. mandshurica and J. regia, a score plot of the OPLS-DA model was constructed ( Figure 8B). We found that the samples of J. mandshurica and J. regia had good separation with respect to the x-axis, which indicated remarkable differences in JMM and JRL. With the combined S plot and VIP (Table S6), by filtering the compounds with VIP > 1, we discovered that SE, GA, MM, and 3-OL were the major contributions to the variance between the two species.
In the above results, we observed a significant difference in the contents and the compositions of the 15 compounds between J. mandshurica and J. regia. Notably, the total content of diarylheptanoids was much higher in J. mandshurica, whereas phenolic compounds content was higher in J. regia.
to the x-axis, which indicated remarkable differences in JMM and JRL. With the combined S plot and VIP (Table S6), by filtering the compounds with VIP > 1, we discovered that SE, GA, MM, and 3-OL were the major contributions to the variance between the two species.
In the above results, we observed a significant difference in the contents and the compositions of the 15 compounds between J. mandshurica and J. regia. Notably, the total content of diarylheptanoids was much higher in J. mandshurica, whereas phenolic compounds content was higher in J. regia.

Conclusions
In this study, we successfully developed a highly efficient and reliable UPLC-TQ-MS method based on chemical profiling to facilitate simultaneous determination of eight diarylheptanoids and seven phenolic compounds within 8 min in five plant parts of J. mandshurica and the epicarps of J. regia. By the analysis of the contents of the 15 compounds in the epicarps of J. mandshurica for different growth periods, we found that the total content of diarylheptanoids and phenolic compounds is the most abundant in June, gradually decreases in July, and then levels off in August; the most suitable harvest time for epicarps in Harbin ranges from Day 12 to Day 27. In addition, leaves are a valuable part, having a variety of chemical constituents and bioactivities that should be studied further. To the best of our knowledge, the differences in the diarylheptanoids and phenolic compounds between different Juglans species have not been previously investigated. The results showed that the profiles of compound contents in the epicarps of J. mandshurica are different from those of J. regia. This paper provides a practical and simple quantitative method and references for further research on the intrinsic qualities of J. mandshurica.
Supplementary Materials: The following are available online at www.mdpi.com/xxx/s1. Table S1: Contents (μg/g) of eight diarylheptanoid compounds from different growth periods, medicinal

Conclusions
In this study, we successfully developed a highly efficient and reliable UPLC-TQ-MS method based on chemical profiling to facilitate simultaneous determination of eight diarylheptanoids and seven phenolic compounds within 8 min in five plant parts of J. mandshurica and the epicarps of J. regia. By the analysis of the contents of the 15 compounds in the epicarps of J. mandshurica for different growth periods, we found that the total content of diarylheptanoids and phenolic compounds is the most abundant in June, gradually decreases in July, and then levels off in August; the most suitable harvest time for epicarps in Harbin ranges from Day 12 to Day 27. In addition, leaves are a valuable part, having a variety of chemical constituents and bioactivities that should be studied further. To the best of our knowledge, the differences in the diarylheptanoids and phenolic compounds between different Juglans species have not been previously investigated. The results showed that the profiles of compound contents in the epicarps of J. mandshurica are different from those of J. regia. This paper provides a practical and simple quantitative method and references for further research on the intrinsic qualities of J. mandshurica.